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Preface 



The constant scaling of CMOS VLSI technologies has produced orders of 
magnitude of increase in packing density and operation frequency of the 
circuits, moving more of the signal processing in to the digital domain. This 
put increasingly difficult demands on the design of analog-to-digital/digital- 
to-analog converters, as the partition between analog and digital sub-circuits 
change. The robustness against circuit non-idealities has made Delta-Sigma 
modulators one of the best choices for high-resolution data conversion. 
Despite the fact that Delta-Sigma modulation is a mature subject, there are 
still many unanswered questions. While a detailed coverage for all these 
questions in a single text is impossible, an honest effort has been made in 
this research monograph to address the exact analysis and speed limitations 
of Delta-Sigma modulators. 

Due to the incorporation of a highly non-linear element (the quantizer) in 
the feedback loop, the exact analysis of Delta-Sigma modulators is a very 
challenging task. In this text, rigorous analysis of Delta-Sigma modulators 
(both multi-stage and single-stage) with rational dc inputs and non-zero 
initial conditions are presented. It is concluded that for a first-order 
modulator the binary quantizer error exhibits a purely discrete power 
spectrum and the conventional white noise assumption is strictly violated. 
Second-order modulators produce the same first-order statistics as those of 
the white noise, but they still fail to generate an uncorrelated quantizer error 
sequence. Third-order modulators are capable of complying with the signal- 
independent white noise assumption provided that an irrational initial 
condition is invoked on the first accumulator. 

It has been well known that Delta-Sigma modulators are non-linear 
chaotic systems, and two infinitesimally small different initial conditions 
may result in completely different and divergent outputs. This book presents 
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the first exact analysis results (to the best knowledge of the authors) that 
rigorously prove the effectiveness of the initial condition in terms of the 
modulator tonality performance. The irrational initial condition principle 
rigorously analyzed here eliminates the need for dithering or chaotic 
operation in which reduction in the limit cycle oscillations and tones comes 
with some performance penalties and hardware overhead. 

During the development of this research, it has become clear that the 
elimination of limit cycle oscillations and tones in some applications are of 
paramount importance. Such a need, for instance, occurs in Fractional-N 
PLL frequency synthesis applications, which may be used to implement a 
local oscillator in a wireless transceiver. In Fractional-N PLL applications, 
the dc input to the modulator poses a significant problem and is responsible 
for the generation of significant amount of spur signals. In this book, an 
important application of our novel mathematical analysis in a Fractional-N 
synthesizer is presented as a solution to cancel spurious signals. The solution 
is simply to set “1” LSB initial condition on the first accumulator, so that the 
initial condition resembles an irrational number eliminating unwanted 
spurious components. 

Increased integration capabilities provided by current-state-of-the-art 
CMOS technologies as well as constant scaling trends encourage analog 
designers to employ more parallelism in Delta-Sigma based A/D converter 
designs. In the architecture domain, the research monograph introduces a 
new method to obtain efficient architectures for time-interleaved Delta- 
Sigma modulators. Additionally, a novel parallel conversion technique, 
which we have dubbed as the “Zero-Insertion Time-Interleaving” concept is 
also proposed. In this approach, the input only needs to be sampled at the 
operating frequency of the parallel channels. Thus, the high sampling rate 
input multiplexer involved in the regular Time-Interleaved approach is 
completely eliminated resulting in further significant reductions in the 
hardware complexity. Such a multiple channel Delta-Sigma modulator may 
be very useful for implementing high-resolution converters for wide 
bandwidth input signals, at the expense of moderate increase in the hardware 
cost. 

This research monograph resulted from a doctoral study completed by the 
first author at the University of Westminster, London, UK, while working 
under the supervision of Prof. Izzet Kale. 



Miicahit Kozak, Izzet Kale 
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INTRODUCTION 



1. INTRODUCTION 

The introduction of oversampling noise shaping technique in [1] by Inose 
and Yasuda has resulted in a very popular method for converting signals 
between the analog and the digital domains. This principle has led to the 
rapid development of robust converters known as delta-sigma^ (AX) 
modulator based converters [2]. In fact, an alternative description of this 
concept was disclosed even earlier in a patent by Cutler [3], which is often 
regarded as the first-order AX modulator with the error-feedback topology. 
Although the AX concept has existed for a long time, there are still many 
problems associated with these modulators. Among these problems, this 
research monograph is mainly concerned with the exact analysis of AX 
modulators, and the exploitation of these analysis results in Fractional-N 
frequency synthesis applications, as well as novel topologies for time- 
interleaved AX modulators. 

1.1 Aims and Motivations 

In AX modulator based data conversion, negative feedback combined 
with a coarse quantization at a high sampling rate shapes the spectrum of the 
quantization noise away from the baseband frequencies [4]. The input signal 
is sampled at a frequency that far exceeds the Nyquist rate to spread the 



‘ In the literature, they are also referred as sigma-delta (LA) modulators. Throughout this 
research monograph, however, we will use the original name, delta-sigma, as it was first 
proposed in [1]. 
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quantization noise over a bandwidth, which is much larger than the signal 
band [5]. The quantized signal is the integrated (sigma) version of the 
difference (delta) between input signal and the analog representation of the 
binary coded output. As a result of noise shaping, only a small portion of the 
noise power lies within the signal band. 

One significant advantage of A£ modulators is that they can achieve 
high- resolution conversion without having component match precision as 
required by the conventional Nyquist rate converters such as the flash 
Analog-to-Digital (A/D) converters. Moreover, oversampling A£ converters 
are well suited for on-chip Very Large Scale Integration (VLSI) 
implementation since they comprise very few components. 

In AZ data conversion, the overall resolution is mainly determined by 
three factors: the oversampling ratio, number of internal quantization bits, 
and order of the noise shaping function [2], [5] and [6]. In practical 
implementations these three degree of freedoms cannot be chosen arbitrarily 
due to requirements imposed by the desired resolution, signal bandwidth, 
available hardware budget (i.e., chip area, power consumption, computation 
speed, and building block count) and other practical limitations. For 
example, it is very common in the art to use a single-bit internal quantizer, 
which consists of a comparator and its associated 1-bit Digital-to- Analog 
(D/A) converter. Therefore, the sampling rate usually needs to be several 
orders of magnitude higher than the Nyquist rate for high-resolution 
conversion using a single-bit internal quantizer. Consequently, this technique 
has been extensively used in relatively low bandwidth applications such as 
speech and audio processing with excellent results [5]-[8]. On the other 
hand, high-resolution conversion of wide band signals using AZ modulators 
may require very high-speed circuitry, which may be very difficult to 
achieve with the current circuit technologies. 

One direct solution to overcoming this problem may be to use a higher 
order modulator with a reduced oversampling ratio so that the circuit speed 
providing the desired signal-to-noise-ratio at the output can be reduced to an 
attainable level. However, this can negate another essential advantage of AZ 
modulators since the anti-aliasing filter does now need a much sharper cut- 
off. In addition, the decimation filter becomes more complicated with the 
increase in noise shaping order [5], [7]. Another possible solution is the use 
of multi-bit internal quantization, which introduces some important pitfalls 
such as requirement of a highly linear feedback D/A converter [7]. The 
multi-bit output also increases the hardware complexity of the digital 
decimation filter. As an alternative solution to this, an attempt has been 
made in this research monograph to provide efficient architectures for time- 
interleaved AZ modulators, which exploit parallelism to mitigate the speed 
limitation of these modulators. It should also be mentioned that this research 
monograph is aimed at providing a method for architectural development of 
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time-interleaved structures. Although the limitations of our method are 
explored in a comparative manner to the existing time-interleaving approach, 
there are still some practical issues that need further research. 

Another significant problem associated with AZ modulators is their 
rigorous analysis, which is made difficult by the presence of a coarse 
quantizer in negative feedback. The quantization error^ spectrum from a AZ 
modulator is usually non-white, unlike what the linear assumption predicts. 
Especially, in an all-digital implementation it is not random noise like, but 
purely deterministic, while in analog implementations it is partially 
randomized by the existing circuit noise. The resulting tones and oscillations 
in the output quantization noise^ spectrum are very undesirable in almost 
every application. Needless to say, the need exists to determine the quantizer 
error distribution through an exact analysis and compare it against the 
traditionally made white noise assumption. It is therefore a subject of this 
research monograph to provide exact analysis of AZ modulators in their 
digital implementations. 

Apart from data conversion, AZ modulators are also used in the 
Fractional-N Phase-Locked Loop (PLL) frequency synthesis applications. 
The simulation of such systems is very difficult because of its mixed-signal 
nature. Circuit simulations (for example using SPICE [9], [10]) is not run- 
time efficient. This is because the difference in time resolutions between the 
output and input of the Fractional-N PLL is very large, and therefore very 
small time steps should be taken to obtain realistic simulation results. 
Although circuit simulation is more accurate and is a vital step towards the 
real implementation, it is not a good starting point to test the spurious 
suppression performance of a particular modulator type. Higher-order 
modulators produce low level of tones and oscillations when an input of 
sufficient activity is employed. In Fractional-N applications, however, the 
input to the modulator is the worst it could ever be, i.e. a constant, resulting 
in spurious components. These periodic contents can be eliminated by using 
our theoretical analysis results — appropriately setting the initial conditions 
in a AZ modulator. However, in order to establish the link between our 
theoretical findings and the Fractional-N spurs a time domain difference 
equation simulation model is needed that can process enormous amounts of 
data. This research monograph is partly aimed at developing a behavioral 
level simulation environment for a complete Fractional-N system, which 
could be used for these purposes. 



^ In this research monograph, the quantizer (or quantization) error is referred as the difference 
between the quantizer output and input. 

^ In this research monograph, we define the output quantization noise as the difference 
between the modulator output and input (i.e., it is essentially the quantizer error sequence 
shaped by the noise transfer function of the modulator). 
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The exact analysis of modulators performed here is confined to only dc 
inputs. There are three main reasons for this. Firstly, rigorous analysis of 
modulators for general time varying inputs is very challenging (at least with 
our method), and there is no reported results in the art concerning this issue 
so far. Secondly, our particular application area, namely Fractional-N PLL 
frequency synthesis, requires that the input be a constant. Lastly and most 
importantly, it could be well argued that if a modulator does not produce any 
tones and oscillations for dc inputs (worst case inputs), then it will be tone 
free for all other input types too. A logical basis for this is that for more 
active inputs the quantizer error becomes more randomized, as expected. 
Note that this conclusion was arrived to through thorough experimental 
verification of our findings. 

1.2 Original Contributions 

This research monograph contains a number of significant contributions 
to the state-of-the-art, which are as follows: 

- A novel rigorous analysis has been performed for the first-order AZ 
modulator with rational dc inputs and non-zero initial condition. It was 
shown that it is impossible to have a tone free spectrum in a first-order 
loop regardless of the initial condition. 

- The rigorous analysis for rational dc inputs and non-zero initial 
conditions has been extended to the higher-order modulator case (both 
multi-stage and single-stage). It has been proven that provided that the 
modulator order is greater than two an irrational initial condition on the 
first accumulator provides a tone free output spectrum. It is also worth 
mentioning that the irrational initial condition principle, analyzed and 
proposed in this research monograph, is a dramatic departure from the 
traditional known methods. There have been mainly two approaches to 
reduce the idle tones at the output of a modulator: dithering and chaotic 
operation. Both of these methods have disadvantages such as requirement 
of additional hardware, jeopardizing stability and increasing the in-band 
noise power. In these methods, the reduction in the tonal performance 
does not come without a cost. The complete elimination of the tones 
achieved with the irrational initial condition principle, however, does not 
come with a performance penalty nor does it require any additional 
hardware. To the author’s best knowledge, this research monograph 
presents the first rigorous analysis of AZ modulators with non-zero initial 
conditions. 

- One question that immediately comes to the mind is how to represent an 
irrational initial condition in the digital domain; otherwise the proposed 
approach would be impractical. Implementation aspects of the irrational 
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initial condition in digital implementations are thoroughly discussed. To 
this effect, we have introduced a novel measure, which is referred to as 
the Noise Shaping Quality factor. Using this measure, for a given 
modulator order, the optimum range of the accumulator size is 
investigated where the irrationality assumption is still acceptable. 

- The Noise Shaping Quality factor can vary between 0 and 1 . It measures 
the difference between the theoretical (assuming an ideal signal- 
independent white noise model for the quantizer) and the actual noise 
shaping performance of the modulator. For a modulator that produces 
near ideal white noise quantizer error, the Noise Shaping Quality is close 
to 1, where for a tonal modular it reduces toward zero. 

- A novel pipelined implementation for a digital multi-stage modulator is 
proposed, which is to be used in a Fractional-N PLL frequency 
synthesizer. The operation frequency is almost doubled through using 
pipelining but without appreciable hardware overhead. 

- A behavioral level computer simulation is developed for Fractional-N 
PLL synthesizers. Besides forming an initial design step toward 
appropriately setting the system parameters, another main advantage of 
this simulator is that it enables rapid assessment of the spurious 
suppression performance of a AZ modulator in a complete Fractional-N 
system. Specifically, by means of this simulator it is shown how the 
spurious components are prevented with the aid of our proposed 
irrational initial condition principle. 

- Efficient architectures for time-interleaved AZ modulators are also 
proposed and explored, for the sake of overcoming the speed limitations 
of oversampled converters. Additionally, a new time-interleaving 
approach is proposed, referred to as the zero-insertion interpolation time- 
interleaving approach, which provides considerable reduction in the 
implementation. 

1.3 Outline of the Research Monograph 

The remainder of this research monograph is organized as follows: 
Chapter 2 presents background material and a review of the state-of-the-art 
achievements in the field of AZ modulation required to follow the work 
presented in the subsequent chapters. 

Chapter 3 deals with the exact analysis of the multi-stage modulators 
with dc inputs and non-zero initial conditions. Chapter 4 extends this exact 
analysis to the higher-order single-stage case. 

Chapter 5 is devoted to the behavioral level simulation of Fractional-N 
PLL frequency synthesizers. A digital implementation for a pipelined fourth- 
order multi-stage modulator is also presented and explained in this chapter. 
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Chapter 6 introduces a new method for designing hardware efficient 
time-interleaved AS modulators. Several simulations are presented to 
support our theoretical developments. 




Chapter 2 

BASIC PRINCIPLES OF DELTA-SIGMA 
MODULATION 



2. BASIC PRINCIPLES OF DELTA-SIGMA 
MODULATION 

Despite the fact that real world signals in electronics are analog in nature, 
it is often desirable to perform the required processing in the digital domain. 
The reasons for this shift from analog to digital processing are the improved 
accuracy, reliability and flexibility of digital circuits, as well as the low-cost, 
high-speed and low-power production offered by modern digital VLSI 
technologies. It is therefore very important to design and develop high- 
resolution A/D and D/A converters that are capable of exploiting the 
advantages of digital processing. 

Recently, oversampled AL modulators have received a lot attention for 
high-resolution data conversion applications, such as high-quality digital 
audio, instrumentation and measurement, and integrated transducer and 
sensor applications. [6]. In this chapter, the basic principles of these 
oversampled modulators are briefly explained along with a survey on the 
state-of-the-art architectures and analysis. We start first by discussing the 
classical Nyquist rate converters and their limitations [11]. The 
oversampling principle is introduced next, which can be used to improve the 
resolution of the conversion [5]. The fundamental basics of AZ modulators, 
which combine oversampling and noise shaping techniques, are described 
later. We shall look at the performance modeling of simple first and second- 
order systems followed by discussions on typical higher-order architectures. 
Oversampling noise shaping techniques have been extensively investigated 
and reported in the literature [6], [7]. The review material presented in this 
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chapter is solely intended to provide a brief survey of the state-of-the-art, 
rather than a full description, so as enable the reader to follow the author’s 
contributions, which are described in the later chapters. 

2.1 Nyquist Rate Converters 



A Nyquist rate converter generates a series of output levels where each of 
these levels has a one-to-one correspondence to a particular input value. 
Shown in Figure 2.1 are the primary functions of a Nyquist rate A/D 
converter"^ that converts an analog input signal into a series of digital codes. 
The continuous time input signal is sampled at a uniform rate /; , and 

then quantized to a finite set of output levels. The effect of the sampling 
process is to create repeated versions of the signal spectrum at multiples of 
the sampling frequency 112]. According to the well-known Nyquist sampling 
theorem, the signal must be band limited to half the sampling frequency, so 
that there is no interference between these repeated versions of the signal 
spectrum. In other words, the sampling frequency /, must be greater than 
twice the signal bandwidth (i.e., >2/^) to avoid aliasing. An anti- 
aliasing filter preceding the sample and hold (S/H) is therefore required to 
ensure that the signal is indeed band limited. 



analog input 







S/H 



Figure 2-1. Block diagram of a Nyquist rate A/D converter 

Nyquist rate converters are usually classified according to the number of 
clock cycles, which are required to complete a single conversion cycle. For 
instance, flash type converters need only one clock cycle to perform the 
quantization and are the standard approach of achieving very high-speed 
A/D converters. The input signal in an A^-bit flash A/D is compared 
simultaneously with all possible quantization levels by using 2^ 
comparators connected in parallel. The limitation of this type of converter, 
however, is in their resolution, because of the matching and accuracy 

Most of the discussion here is referred to A/D conversion, but much of it directly applies to 
D/A conversion too without loss of generality. 
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requirements put on some analog components (resistors or capacitors). For 
an additional bit of resolution, the matching requirements and the circuit 
complexity will double. Matching of components to greater than 10 bits or 
equivalently to one part in one-thousand (2“’) is almost impossible to 
achieve in the current VLSI technologies, unless very expensive laser 
trimming or self-calibration techniques are used [8]. Furthermore, the 
complexity of the implementation exponentially increases with the number 
of bits of resolution. Consequently, flash type A/D converters have been 
used in high-speed applications with at most 10-bit resolution. 

To mitigate these limitations of flash converters, a number of alternative 
architectures have been proposed in the literature and used in the industry 
[13]. For example, in a subranging converter there are two steps to complete 
a conversion cycle. In the first step, N/2 most significant bits (MSBs) are 
determined by using an N/2 -bit flash A/D sub-converter. To determine the 
remaining N/2 least significant bits (LSBs), the digital signal obtained in the 
first step is converted back to the analog domain using an N/2 -bit D/A sub- 
converter, and subtracting the equivalent analog MSBs from the input signal. 
This difference is then subject to a fine quantization by means of another 
N/2 -bit A/D sub-converter. With this method, rather than requiring 2^' 
comparators for an N -bit flash converter, only comparators are 

needed in subranging converters. However, note that the latency is now 
increased by two, which results in the decreased conversion bandwidth for a 
given clock frequency. Additionally, the component matching requirements 
are not relaxed. Although, this approach makes use of N/2 -bit flash A/D 
sub-converters, the accuracy of these converters would still need to match to 
the overall resolution [5], [13]. 

In some other methods, where one-bit is resolved per clock cycle, only 
one comparator is required. Converters employing the successive 
approximation method fall into this category [13]. The successive 
approximation converter applies a binary search algorithm to determine the 
closest digital word to the input signal amplitude. In the first clock cycle, the 
MSB is determined by comparing the input signal with a feedback voltage. 
In the following cycles, the feedback voltage is successively changed in 
accordance with the comparison result obtained in the previous cycle, until 
all bits are determined. Thus, a successive approximation converter needs N 
clock cycles to complete a single conversion with A-bit resolution. 

All of these Nyquist rate converters put stringent requirements on the 
anti-aliasing filter design. This is because the input signal is sampled very 
close to the Nyquist rate and the anti-aliasing filter must have a very sharp 
cut-off transition, which significantly increases the hardware complexity [5], 
[8]. On the other hand, in AS converters the transition bandwidth is much 
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wider due to the oversampling, and consequently the implementation of the 
anti-aliasing filter is significantly simplified [14]. Another important 
advantage of NL converters over Nyquist rate ones is that the component 
matching requirements are also significantly relaxed by use of a coarse 
quantizer (usually single-bit) in a negative feedback loop combined with 
linear filtering. Table 2.1 summarizes the advantages and disadvantages of 
traditional Nyquist rate and oversampled AZ data conversion techniques. 

Table 2-1. Comparison of oversampling AZ and Nyquist rate converters 





Oversampling AZ 


Nyquist rate 


Bandwidth 


low to medium 


high to medium 


Resolution 


high (>16 bits) 


medium (<12 bits) 


Hardware complexity 


low 


high 


Sensitivity to component 


low 


very high 


accuracy 

Anti-aliasing filter 


very simple 


complicated 



2.2 Quantization Noise 

The quantization process involved in A/D and D/A conversion is an 
inherently non-linear operation and introduces errors to the conversion [15]. 
A typical input-output transfer characteristic of a uniform quantizer is a 
staircase function as shown in Figure 2.2, where v[n] and y[n] are the input 
and output sequences, respectively. The difference between adjacent output 
levels is called the quantizer bin width A . Basically, the quantizer maps any 
input amplitude to the closest output level. Clearly, an infinite number of 
input amplitudes can be mapped to a unique output level resulting in loss of 
exact information of input amplitude (and hence the quantization error). One 
of the primary challenges in data conversion applications is to understand 
and correctly model the behavior of these errors [16]. 

The quantizer characteristic in Figure 2.2 is a mid-riser type [5]. 
Throughout this research monograph we will make use of the mid-riser 
quantizer characteristic, however, it should be noted that the mid-thread 
quantizer characteristic is also equally applicable. The quantization error is 
defined as the difference between its output and input (i.e., e[n] =y[n]-v[n]). 
This allows one to write the output as the sum of the input and the 
quantization error. One should immediately recognize that this definition 
does not depend on any assumption and it is not a linear model in the sense 
that the quantization error is not an independent signal but may be strongly 
dependent on the input [11]. This mathematical representation, however, 
becomes an approximate linear model when some assumptions are made on 
the statistics of the quantization error. Although conceptually simple, the 
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quantizer is very difficult to analyze, because of its severe non-linearity. In 
order to make the analysis tractable, the quantizer is often linearized by 
using an input-independent additive white noise model, as shown in Figure 
2.3. 



y 




Figure 2-2. Input-output transfer characteristic of a quantizer 



quantizer 




I 

linear model 



signal-independent white noise 




Figure 2-3. Quantizer and its linear model 

In the white noise model, a number of assumptions are traditionally made 
on the error process and its statistics, which follow [17], [18]: 

- the quantization error sequence e[n] is statistically independent of the 
input sequence, (i.e., e[n] is uncorrelated with v[^] for all n,k ) 
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— the Probability Density Function (PDF) of e[n\ is uniform over the range 
of half the bin width, (i.e., e[n] is uniformly distributed in [- A/2, A/2]) 

— the random variables of the error process are uncorrelated (i.e., e[n] is a 
white noise) 

Now assume a Nyquist rate converter that samples the input signal at the 
Nyquist rate /, =2/^ (/^ is the signal bandwidth) and employs A^-bit 
quantization. Let the maximum output levels of the quantizer be ±b as in 
Figure 2.2. In this case, the bin width is A = - 1). The magnitude of the 

quantizer error e{n\ does not exceed half a bin width provided |v[«]|<2'^"'A 

for all n . Under these circumstances, the quantizer is said to operate in the 
no-overload region. Figure 2.4 and Figure 2.5 show the PDF and the Power 
Spectral Density (PSD) for the quantization error, according to the linear 
model described above. The mean value of e[n] is then 



1 

U=— \ede = 0 

A J 

-A/l 



The average power of the quantizer error is given by^ 



<7, 





( 2 . 1 ) 



( 2 . 2 ) 



Pe(e) 




Figure 2-4. PDF of the quantizer error 



^ Recall that for a zero-mean random process, its variance <7 ‘ is equal to its average power 
[19]. 
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As v[n] and e[n] are assumed statistically independent by the linear 
model, the noise power at the data converter output is equal to that of the 
quantization error power crj . If we treat the input signal as a zero-mean 
random process with average power , then the Signal-to-Noise Ratio 
(SNR) at the output is [5]: 



SNR = mog, 




= lOlogi, 






+ 4.77 + 6.02A^ {dB) 



(2.3) 



SAf) 




Figure 2-5. PSD of the quantizer error 

Notice that for each additional bit in the quantizer resolution, the SNR 
increases approximately by 6 dB. It is therefore very common in the art to 
equate the differences of SNR performances in dB to the equivalent bits of 
resolution. For instance, a data converter with a 3dB better SNR is said to 
have a Vi bit better resolution. 

2.3 Oversampling Advantage 

Oversampling is a technique that improves the resolution of a Nyquist 
rate data converter. This improvement is achieved by sampling the input 
signal at a significantly faster rate than the Nyquist rate. The ratio between 
the sampling rate and two times the signal bandwidth is defined as the 
Oversampling ratio (OSR): 




(2.4) 
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After oversampling, the signal is subject to quantization and then low- 
pass filtered and finally down-sampled to or near to the Nyquist rate^ . The 
low-pass filter eliminates the high frequency quantization noise. After 
filtering the sampling rate can be reduced to the Nyquist rate without 
affecting the SNR performance. The low-pass filtering and down sampling 
operations are collectively known as decimation [20], [21]. 

Assume that the oversampled signal is quantized with A-bit resolution, 
then the average quantization error power is given by (2.2). Obviously, this 
is exactly the same power produced by a Nyquist rate converter. The total 
quantization error power is the same, but its frequency distribution is 
different due to oversampling. The linear model assumes that the quantizer 
error is white, implying that its power is uniformly distributed over -fJ2 
and /, /2 . Figure 2.6 reveals the power spectral density of the quantizer error 
for a Nyquist rate and oversampling converter. In the oversampling case, the 
same amount of error power is spread over a much wider frequency 
bandwidth, leaving a relatively small amount of error power within the 
signal bandwidth. 



SAf) 



J 



Nyquist rate converter 



Oversampling converter 



- f ,/2 



A 



A/2 



■>/ 



Figure 2-6. PSD of the quantizer error in oversampled conversion 

For an oversampled converter, the power spectral density of the quantizer 
error is SAf) = <y] j f, for - f,j2 < f < fj2 , where o] equals to aYI 2 . If one 
filters the quantizer output with an ideal brick-wall low-pass filter with cut- 
off frequency A » the in-band error power at the output of the data converter 
is then 



® In D/A applications, the signal is filtered by an analog low-pass filter, and there is no down 
sampler required. Moreover, an interpolation filter is used to oversample the digital input. 
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a:. 





OSR 



(2.5) 



Since the signal is unaltered after low-pass filtering, its power remains 
the same and the SNR at the output of the oversampled converter is 



V 




+ l01og,„(OSfl)(dB) 



(2.6) 



For each doubling of OSR (increasing the sampling rate by two), the SNR 
increases by 3 dB or equivalently a Vi bit of resolution improvement is 
obtained. In the oversampling scheme, we trade-off speed for better 
resolution. In other words, higher resolutions are obtained by operating the 
circuits at a faster rate. This trade-off becomes more advantageous through 
the use of AX modulation that employs both oversampling and noise shaping 
techniques simultaneously, to good effect. 



2.4 AZ Modulation 



As noted in the previous section, oversampling reduces the in-band 
quantization noise by redistributing a fixed amount of noise power over a 
bandwidth much larger than the signal bandwidth. The in-band noise is 
further attenuated by pushing much of the noise power out of the signal 
band. This is process is known as the “noise shaping technique”, and is 
employed by AX modulators. The output of the AX modulator can then be 
low-pass filtered and finally down sampled to or near to the Nyquist rate, 
like conventional oversampled converters. 

Shown in Figure 2.7(a) is a block diagram for a general AX modulator 
that incorporates a quantizer (preferably single-bit) along with a filter H{z) 
in negative feedback. If the white noise approximation for the quantizer is 
used as described in Section 2.2, the linear model of the modulator is 
illustrated in Figure 2.7(b). Performing straight forward z-domain analysis on 
the linear model, one obtains 



Y{z) = X{z) 



H{z) 
//(z) + l 



+ £(z) 



1 

//(2) + l 



(2.7) 



where A(z), Y{z) and E{z) are the z-transforms of the input, the output 
and the quantization error, respectively. 

In AX modulation, the signal that is subject to quantization is not the 
input signal itself, but a filtered version of the difference between the input 
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signal and the encoded output. The filter H{z), usually called the feed- 
forward filter, in a first-order AS modulator is a discrete-time integrator with 
transfer function //(z) = 2 ''/(l -z''). Using this, for a first-order AS modulator 
we can write 

Y{z) = X{z)z'^-E{z){\-z-'] (2.8) 



e\ti\ 




(b) 

Figure 2-7. General AS modulator: (a) block diagram, (b) linear model 

It follows that in a first-order modulator, the output is given by the 
superposition of the one-sample-delayed version of the input signal, plus the 
quantization noise that is shaped by a first-order differentiation (high-pass 
filtered). The input signal is subject to filtering by the signal transfer 
function (STF) S{z) = z'\ while the quantization error undergoes filtering 
through the noise transfer function (NTF) A^(z) = (l-z'’). Observe that the 
NTF is a high-pass filter with a zero at dc, and thereby reduces the 
quantization noise around dc. To evaluate the SNR performance, one needs 
to find the signal and noise components at the output of the modulator. We 
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make use of the fact that if a stationary random process with power spectral 
density S{f) is input to a linear filter //(/), then the output is another 

stationary random process with power spectral density 5(/)|//(/)|^ [19]. 
Realizing that the magnitude of the STF is 1 (i.e., |5(z)|^ =|z“'P = l) and 
further assuming that y[n] will be filtered by an ideal low-pass filter, the 
signal power at the output remains intact as c ] . In order to find the in-band 
noise power, one should use the following relation [5] 




To find the magnitude of the NTF we let z = and write the 

following 



A^(/) = 1 _ = 2 



2J 



(2.10) 



Taking the magnitude of both sides 
|A'(/)| = 2sin(;5f//,) 

Applying (2.11) into (2.9), we get 



A ( ) 



( 2 . 11 ) 



( 2 . 12 ) 



For large OSRs (i.e., (2/i,//J«l), we can apply a Taylor series 
expansion to the sinusoid term as [22] 



sin(x) = x 

6 



(2.13) 



This reduces (2.12) to 






(2,14) 



The SNR at the output of a first-order modulator is then 
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SNR = \() 



log, 






|- 10 log,,, 



n~ 

T 



+ 301og,„(05/?)(t/S) 



( 2 . 15 ) 




(a) 




(b) 



Figure 2-8. Magnitude spectra for NTFs: (a) linear, (b) log scale 

For a second-order accumulator based modulator, the NTF is 
N{z) = {\-z 'J . The magnitude spectrum of this second-order NTF is shown 
in Figure 2.8 along with first and third-order ones for comparison. Notice 
that higher-order NTF provides more attenuation over low frequencies and 
more amplification over high frequencies. Note that the plots in Figure 2.8 
are for the ideal case. In practice, these NTFs may not be achievable 
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especially when a single-bit quantizer is used 123]. The SNR performance of 
the second-order system can be found in the same way as was performed for 
the first-order case: 



SNR = lOlogio 



|“10Iog,o 






+ 501og,o(05’i?)(c/5) 



(2.1) 



For every doubling of the OSR, the SNR increases by 9 dB and 15 dB (or 
equivalently 1.5 bits and 2.5 bits of resolution enhancement is obtained) in 
the first-order and second-order AZ modulation, respectively. In general, for 
an U -order modulator the SNR improvement is (6L-I-3) dB (or L -1-0.5 bits) 
per doubling the OSR. 

The system architecture for a typical AZ based A/D converter is shown in 
Figure 2.9(a). The first element in the system is the continuous-time analog 
anti-aliasing filter that band-limits the input signal to frequencies less than 
one half the oversampling frequency. For large oversampling ratios, the anti- 
aliasing filter can be as simple as a first-order RC low-pass filter. The 
filtered signal is then processed by a mixed-signal AZ modulator, which 
converts its analog input signal into a low-resolution digital signal. The AZ 
modulator is usually implemented using discrete-time switched capacitor 
(SC) circuitry^. An internal D/A is required to convert the encoded digital 
output back to the analog domain. The analog representation of the output is 
then subtracted from the input signal before being fed into a sampled-data 
filter. Therefore, any non-linearity in the internal D/A would essentially 
appear in the digital output of the modulator. As a result, the overall linearity 
of AZ converter is strongly dependent on the linearity of its internal D/A 
converter. This is the reason why most of the applications have exploited a 
single-bit quantizer and its associated single-bit D/A converter. A 1-bit D/A 
converter is by definition linear [24]. 

This linearity comes from the fact that a 1-bit D/A has only two output 
levels and two points always define a straight line in the input-output 
transfer characteristic. This inherent linearity of 1-bit converters is one of the 
major factors for making use of oversampling AZ techniques with 1-bit 
quantization. The low-resolution digital signal at the output of the modulator 
consists of the input signal plus a large amount of high frequency 
quantization noise. The out-of-band quantization noise is then removed by a 
digital decimation filter, which produces a multi-bit digital representation of 
the input at a sampling rate close to twice the signal bandwidth. 



’ In the switched capacitor realization of AL modulators, a separate sample-and-hold is 
normally not required as the continuous-time signal is sampled by the switches and 
capacitors. 
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mixed signal A£ modulator digital decimation filter 




discrete-time analog processing ^ ■ ► digital processing 



(a) 



digital interpolation filter digital modulaior 




digital processing 4 ► analog processing 



(b) 



Figure 2-9. System architecture for AI based converters: (a) A/D, (b) D/A 

A AS based D/A converter can be realized as shown in Figure 2.9(b). The 
input signal to the system is a multi-bit digital signal, which is sampled at or 
slightly higher than the Nyquist rate. This signal may then be up-sampled to 
an oversampling frequency by padding (OSR-i) zero valued samples 
between adjacent samples, or employing (05/?- 1) zero-order hold. This up- 
sampling process generates images of the original signal spectrum in the 
frequency domain. A digital low-pass filter is then deployed to reject all 
other images, except the original one. The filtered oversampled signal goes 
into an all-digital AZ modulator that produces a 1-bit encoded digital signal. 
As discussed earlier, the reason for choosing single-bit truncation (or 
quantization) is so that a 1-bit D/A converter with excellent linearity can be 
used to convert the digital modulator output into the analog domain. The 
high frequency quantization noise is then attenuated by using a simple 
analog low-pass filter. 
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2.5 Limit Cycle Oscillations and Tones 

Like digital infinite impulse response (HR) filters that operate with finite 
precision arithmetic, A2 modulators (both digital and analog) are also prone 
to display limit cycle oscillations and tones that result in the presence of 
periodic components in the output spectrum. This is not surprising because 
AS modulators are non-linear systems that employ negative feedback. 
Although the linear approximation models the quantizer error as a signal- 
independent white noise process, the quantization error from a AS modulator 
is usually colored and non-white [25], [26]. For dc (constant) inputs the 
quantization error is periodic generating what is commonly called idle 
channel tones. For sinusoidal inputs the quantization error is also periodic 
causing the output spectrum to contain harmonic tones, which are dependent 
on the input amplitude and frequency in a very complicated way [27]. Such 
periodic structures in the quantization error cause severe problems in many 
data conversion and other signal processing applications, and they are 
strictly objectionable as some tone components may well lie within or beat 
down into the baseband. 

At present there are a number of techniques in the state-of-the-art to 
overcome the tone problems. The most commonly used two techniques are 
dithering and chaos [7], [28]. Each of these methods is now explained in 
some detail in the following subsections together with their advantages and 
disadvantages. We also note that in some A/D and D/A conversion 
applications, the tones are shifted or forced to be placed where they are least 
harmful, by superimposing a very small dc offset on the input signal. The 
details of this technique have not been published yet, but the author is aware 
of this technique through private communication with his supervisor [29]. 
Also note that this technique may not be effective for dc inputs. 

2.5.1 Dithering 

Dithering in AX modulation is the act of adding a random (or pseudo- 
random) signal to the input of the quantizer (not to the modulator input) as 
shown in Figure 2.10. Assuming that the added random signal is white, it 
also becomes noise shaped like the quantization error so that the additional 
noise power in baseband is minimized [30]. The purpose of dithering is to 
effectively decorrelate and whiten the quantization error. Typically, the 
dither signal is generated using some sort of pseudo-random generator with 
only a few bits of resolution. In D/A applications, dithering can be 
accomplished by means of an additional digital adder along with Linear 
Feedback Shift Registers (LFSR) that can be used to generate pseudo- 
random signals. For A/D applications, however, a low-resolution D/A sub- 
converter is needed to convert the digital random signal into the analog 
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domain. Therefore, in A/D conversion dithering significantly increases the 
complexity of the implementation. 



dither signal 
d[n\ 




iwn 



quantizer 



H(z) — 






Figure 2-10. Noise shaped dithering of AZ modulators 

The optimum amplitude and PDF for the dither signal in order to 
eliminate tones and oscillations were investigated in [28] and [31]. In [28] it 
was shown that for a second-order single-bit modulator a dither signal with 
± 0.5A (where A is again the bin width of the quantizer) amplitude and 
rectangular PDF (uniform distribution) eliminated the idle tones at the 
expense of the decreasing the maximum SNR by approximately 6 dB. For 
higher-order modulators, however, the concession in peak SNR might be 
much greater, for instance as high as 26 dB in a sixth-order system [31]. 
Such huge reductions in higher-order modulators are due to the reduced 
system stability. In general, the higher the dither amplitude the more 
successfully tones are removed, but at the expense of increased baseband 
noise. In higher-order modulators large amplitude dither signals limit the 
useable input amplitude range. While correctly dithered modulators do not 
suffer from tones and oscillations problem at least to some certain extent, the 
in-band noise level as mentioned earlier is increased by the injection of the 
additional random signal. In addition the maximum input amplitude for 
which the modulator is stable is reduced. In summary, there are three main 
shortcomings in the dithering technique. Firstly, the hardware complexity is 
increased by the additional dither logic. Secondly the modulator has a lower 
peak SNR and reduced dynamic range. Finally, complete elimination of 
tones requires large amplitude dither signal reducing the stability of the 
modulator in a way, which is often difficult to predict, and requires re- 
checking the stability by means of extensive simulations, and in some cases 
by rapid Field-Programmable-Gate- Array (FPGA) prototyping. 
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2.5.2 Chaotic AS Modulators 

Chaotic operation of AS modulators is suggested as an alternative 
method to dithering, and can be accomplished by moving one or more of the 
NTF zeros outside the unit circle [32]. Such an approach is simply 
implemented by modifying the transfer function of the discrete-time 
integrators within a AS modulator as z“7(>-ce“') such that a>l instead of 
unity The advantage of chaos over dithering is that neither a random 
number generator nor a D/A sub-converter is required, making this 
technique especially attractive in AS based A/D applications [33]. 

Introducing chaotic operation breaks up any periodicity in the output, 
which is very desirable in order to reduce tonal behavior. However, the 
degree at which the tones are suppressed is dependent on how far the zeros 
are located outside the unit circle. Although the output spectrum of a chaotic 
modulator cannot be perfectly periodic, in order to completely eliminate the 
tones the zeros should be located far outside the unit circle jeopardizing the 
stability. Unfortunately, higher-order modulators cannot be made totally 
tone-free without incurring severe SNR penalties. Several researchers [281 
and [32] have performed extensive simulations and reported that reduction in 
tonality is obtained at the expense of increased in-band noise power, rather 
than decreasing the average tone power. For example, in a fourth-order 
modulator, the SNR decreases as much as 66 dB by introducing chaos [28]. 
In general, chaotic operation of modulators is less efficient in whitening the 
quantizer error when compared to dithering. 

2,6 Higher-order AZ modulators 

As with any other non-linear systems employing feedback, AS 
modulators have the potential to become unstable. In this section, we will 
shortly explain two different architectures for realizing higher-order 
modulators. The first one is the single-stage architecture in which there is 
only one quantizer in the feedback loop [34]. The second approach is the 
multi-stage architecture where several independent lower-order modulators 
are cascaded to construct a higher-order modulator. The final output is 
obtained through a noise cancellation circuit that operates on the output of 
each stage to cancel out all quantization error terms with the exception of the 
last stage [35]. 



Note that when the feedback coefficient of an integrator a is smaller than 1, it becomes a 
leaky integrator. 
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2.6.1 Single-stage Architectures 

It was long believed that higher-order modulators were impossible to 
stabilize when single-bit quantization was used. Fortunately, Lee [36] has 
proposed one of the first approaches for realizing higher-order single-stage 
modulators in 1990. A linear model for a third-order modulator based on his 
topology is given in Figure 2.11. By performing linear analysis, one can 
show that 



Y(z) = X(z) ^ + E(z)~ 



(i-z-7 



(2.17) 



Note that in this case, both the STF and NTF are HR transfer functions 
rather than finite impulse response (FIR) transfer functions. The purpose of 
HR type filtering is to reduce and flatten the out-of-band gain of the NTF. 
Unfortunately, the linear analysis gives an inaccurate estimation on the 
stability. It should be noted that a rigorous stability theory for single-stage 
higher-order single-bit modulators does not exist in the literature, leaving no 
choice to the designers, but to resort to computer simulations. Extensive 
simulations showed that in higher-order single-stage modulators the stable 
operation can only be achieved by a clever choice of modulator coefficients 
as well as restricting the maximum input amplitude [34]. In order to order 
obtain a stable operation, as a general rule of thumb, it was suggested to 
keep the peak gain of NTF to less than 2, i.e., max{A^(/)|}< 2 for all 

f e[- f j2,fJ2] [36]. The stability is also dependent on the choice of the 
NTF filter type. Possible choices are Butterworth or inverse Chebyshev filter 
types [37]. For a system that uses Butterworth pole alignments, the input 
range over which the modulator is stable is typically between 60-80 % of the 
reference feedback voltage. Throughout this research monograph, the 
maximum input signal amplitude with respect to the feedback voltage 
beyond which the modulator is no longer stable is referred to as the 
“maximum depth of modulation”. Other examples of stable higher-order 
single-stage modulators with a single-bit quantizer are described in [37]-[39] 
by Adams together with design methodologies for choosing the modulator 
coefficients. 
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Figure 2-11. Linear model for a third-order single-stage modulator 

2.6.2 Multi-stage Architectures 

In multi-stage approach, low-order AS modulators are connected in 
cascade each with its own quantizer. The advantage of this method is that as 
lower-order modulators (first or second-order) are unconditionally stable, the 
stability of the overall system is guaranteed. Such arrangement has been 
called Multi-stAge-noise-SHaping (MASH) modulator after it was first 
proposed in [40]. As an example, a third-order modulator constructed by 
cascading three first-order modulators is illustrated in Figure 2.12. The input 
is fed to the first modulator and the negative of the quantization error from 
the first and second stage is used as the input to the second stage and third 
stage, respectively. Finally, the outputs from each stage are combined 
together to form a final output after passing through the appropriate filtering. 
The output from each stage can be written in the z-domain as 



Y,(z) = X(z) + E,(z){\-2-'] 


(2.18) 


Y,(z) = -E,(z) + E,{z)(\-z-') 


(2.19) 


Y,(z) = -E,(z) + E,(z){l-z-') 


(2.20) 



where X{z), Ey{z) , EyO) z-transforms of the input and 

quantizer error sequences from the first, second and third stage, respectively. 
The final output is computed by 



y( 2 ) = r,(z)z-' + (y,i 2 )z-' + r,(z) - r,(z)z-' )(i - z-‘ ) 



( 2 . 21 ) 
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This sum results in the cancellation of the quantizer error terms from the 
first two stages. In order to see this, we substitute (2.18)-(2.20) into (2.21) to 
obtain 

Y{z) = x{z)z-^ + £ 3 (z)(i - z "' ) (2.22) 




Figure 2.12. Linear model for a third-order MASH modulator 

Although MASH type modulators are appealing in their improved 
stability and high dynamic range performance, there are some disadvantages 
in using these architectures. A brief comparison between MASH and single- 
stage architectures is listed in Table 2.2. In A/D applications, the noise 
cancellation network is implemented in the digital domain. Unfortunately, 
any mismatches between analog integrator and digital differentiator result in 
the incomplete cancellation of the quantizer error terms from earlier 
sections. This causes some poorly shaped quantization error to leak to the 
final output increasing the in-band noise power. Additionally because of 
imperfect cancellation the tonal performance of a higher-order MASH 
modulator is determined by those of the lower-order modulators in the 
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cascade structure [7]. In D/A applications, the cancellation network can be 
perfectly implemented in the digital domain, but the multi-bit nature of the 
final output calls for a multi-bit D/A converter, which has its own linearity 
disadvantages in VLSI implementations. Consequently, the performance of 
MASH modulators is more sensitive to the analog imperfections compared 
to single-stage architectures. Nevertheless, these sensitivity requirements are 
far below than those required in a Nyquist rate converter, and MASH 
architectures have been very popular in realizing high-resolution A/D and 
D/A converters for a number of aggressively specified applications [5], [6], 
and [8]. 

Table 2-2. Comparison of higher-order single-stage and MASH AI architectures 





advantage 


disadvantage 


Single-stage 


• single-bit output 


• only conditionally 




• relatively simple circuit 


stable 




design 


• difficult loop filter 




• less sensitive to 


design 




matching 


• no rigorous theory for 
stability 


MASH 


• guaranteed stability 


• requires matching 




• maximum input 


between analog and 




amplitude can be made 


digital transfer 




almost equal to 


functions 




reference feedback 


• multi-bit output 



2.7 Multi-bit AZ Converters 

As noted before, the reason for preferring single-bit quantization in AS 
modulators is that its associated single-bit internal D/A converter does not 
require precision component matching, making the implementation very 
robust. However, the resolution at the single-bit output is limited at a given 
oversampling ratio. The resolution could be improved by utilizing a higher- 
order loop topology, but these improvements diminish very rapidly because 
of the instability problems. Furthermore, these improvements in a higher- 
order modulator are most beneficial when the OSR is very large. At low 
OSR, there is a little increase in SNR performance when the loop order is 
increased [41]. In addition, because there is large amount of out-of-band 
quantization noise present in the output spectrum when a single-bit quantizer 
is used, the design of the decimation filter in A/D applications or the analog 
low-pass filter in D/A applications is further complicated and made more 
difficult. 

One of the solutions to these problems is to use a multi-bit quantizer. The 
primary advantage of multi-bit quantization is that the total noise power at 
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the output is substantially reduced compared to that of single-bit quantizer, 
as a rule of thumb typically 6 dB per additional increase in the bit count of 
the quantizer. Therefore, the resolution of AZ based converters can be 
increased without increasing the OSR or utilizing a higher-order loop. 
Besides this, modulators using multi-bit quantization display less tonal 
behavior. Naturally, the behavior of the quantizer error sequence in a multi- 
bit modulator is more accurately predicted by the white noise model [42]. 
Stability analysis based on the linear model also becomes more accurate. 
The difficulty associated with a multi-bit quantizer, however, is that it 
cannot be easily fabricated with sufficient linearity essential for high- 
resolution conversion. 

Fortunately, there are some methods present that alleviate the linearity 
problem encountered when a multi-bit quantizer is used. These include 
dynamic element matching (DEM) [42]-[44], mismatch shaping [45] and 
digital correction techniques [46]. 

2.8 State-of-the-art Analysis 

We have seen that the white noise model for the quantizer error sequence 
greatly simplifies the analysis of AZ modulators, thereby provides very 
important insights on how the system works. Because the model replaces the 
deterministic quantizer error sequence by an input-independent stochastic 
random process, conventional linear system theory can be used in the 
analysis. While it delivers astonishingly accurate results for the SNR 
performance, this linear model, however, completely fails to explain the 
tones and oscillations present in the output spectrum. Furthermore, higher- 
order modulators exhibit stability problems, which are not predicted by the 
linear model [47]. 

The white noise model has become very popular for the analysis of 
quantization error in the pioneering work of Bennett [17], who first 
developed conditions under which the white noise approximation can 
reasonably represent the true behavior of the quantizer error. Specifically, 
Bennett proved that if the following conditions hold then the input- 
independent white noise model could be safely applied: 

- the quantizer input is in the no-overload region 

- the quantizer bin width A is asymptotically small (i.e., the resolution of 
the quantizer is very large, for example more than 8-bits) 

- the joint PDF of the quantizer input signal at different sample times is 
smooth 

While these conditions can be easily secured in Nyquist rate data 
converters [48], where the bit resolution of the quantizer is very high, they 
are strictly violated in oversampled AZ modulators, because [15] 
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- it is usually not known whether the maximum quantizer input amplitude 
remains in the no-overload region 

- the quantizer usually has a few bits of resolution (one-bit most desirable), 
and therefore the quantizer bin width A cannot be regarded very small 

- the quantizer input at a specific time is also dependent on its previous 
value due to the feedback action, therefore the joint PDF will not be 
smooth 

The invalidity of the signal-independent white noise approximation of 
the quantizer error sequence was the major reason that forced researchers to 
develop more accurate models for AX modulators. In fact. Gray was the first 
scientist who developed an exact analysis method for the first-order AX 
modulator with dc inputs in [49]. In this approach, a non-linear difference 
equation describing the first-order modulator was derived and then using this 
equation long-term statistics of the quantizer error sequence were obtained. 
Gray’s work lead to an important demonstration that in the first-order 
modulator driven by a constant (dc) input, the quantizer error sequence is 
nonwhite and the white noise approximation is completely invalid. Gray et. 
al. [50], [51] then performed exact analysis for higher-order MASH AX 
modulators by combining the non-linear difference equation approach with 
the characteristic function method, where they assumed that the input to the 
modulator is an irrational fraction of the reference feedback level. Using the 
same approach, higher-order single-stage multi-bit modulators were also 
analyzed by He et. al. in [52], [53] with irrational dc inputs. The reason for 
the irrational dc input assumption is that if the amplitude of the signal is 
chosen according to some continuous PDF, the resulting amplitude will be 
an irrational number with probability one. Of course, this can only happen in 
an analog implementation, where there is always an uncertainty in voltage 
levels, and hence rational numbers cannot occur. 

Another method is to use a Fourier series expansion for the error 
sequence. Using this approach, Rangan [54] has developed exact formulas 
for the quantization noise in a second-order two-bit modulator. Unlike the 
previous exact analysis, this method directly determines the quantization 
noise spectrum by using an open loop model and using a Fourier series 
representation of the quantizer error function. One improvement achieved in 
this method is that it is valid for both rational and irrational inputs. However, 
traditional Fourier series expansions are not well suited to the discontinuities 
of the quantizer transfer function, and may give inaccurate results. Most 
importantly, this method is restricted to only the second-order case, because 
the algebra involved in higher-order integration and series expansion seems 
very difficult to solve if not impossible. Hence, to the author’s best 
knowledge no one has been able to extend this method to the higher-order 
cases. Galton [55], [56] has developed rigorous mathematical formulas for 
the autocorrelation of the non-overload quantizer error sequence in a class of 




30 



BASIC PRINCIPLES OF DELTA-SIGMA MODULATION 



modulators with arbitrary inputs. His formula covers more general inputs, 
but the shortcoming of this approach is the assumption of a small Gaussian 
noise superimposed on the input, which is only reasonable in analog 
implementations. While in the analog implementations the inevitable 
thermal noise or some other circuit noises at the analog front-end of the 
modulator could be well approximated by a very small amplitude Gaussian 
random process, in the digital implementation on the other hand such 
assumptions are improper. Poorfard et. al. [57] also developed a theoretical 
frame-work to analyze AZ modulators with zero mean stochastic inputs. This 
approach was based on the assumption that the input to the single-bit 
quantizer was a Gaussian random process with zero-mean. This is obviously 
only valid especially when the input to the modulator is a stochastic signal 
and the modulator order is very high. 

Most of these theoretical developments concentrated on analog 
implementations, where either there exists some circuit noise or the 
amplitude of the input signal is irrational. On the other hand, digital 
implementations of modulators severely suffer from poor randomization of 
the quantizer error sequence. This comes about from the absence of circuit 
imperfections (such as thermal noise or input referred circuit noise) that 
helps to partially decorrelate the successive quantizer error samples. 
Additionally, the amplitude of the input signal as well as all the state 
variables can only be represented by rational numbers relative to the 
reference feedback level, which result in strong idle tones and limit cycle 
oscillations in the output spectrum. To date, the exact analysis of modulators 
in their digital implementations has been unexplored. The next two chapter 
of this research monograph focus on this open issue and provide the first 
rigorous analysis (according to the best knowledge of the author) of AZ 
modulators (both single-stage and MASH) with rational dc inputs. 




Chapter 3 

ANALYSIS OF MASH DELTA-SIGMA 
MODULATORS WITH DC INPUTS 



3. ANALYSIS OF MASH DELTA-SIGMA 
MODULATORS WITH DC INPUTS 

In this chapter, we have performed rigorous analysis for higher-order 
MASH AS modulators, which are built out of cascaded first-order stages, 
with rational dc inputs and non-zero initial conditions. Asymptotic statistics 
such as the mean, average power and autocorrelation of the binary quantizer 
error are formulated using a non-linear difference equation approach [58]. 
An important topic of interest considered here is the Fractional-N PLL 
frequency synthesis application, where the input to the modulator has to be a 
rational constant. It has been mathematically shown that, regardless of the 
initial conditions, first-order and second-order 1-1 MASH modulators with 
rational dc inputs cannot sufficiently randomize the quantizer error samples, 
and therefore are not suitable for Fractional-N synthesis applications. An 
irrational initial condition imposed on the first accumulator of a third or 
higher-order MASH modulator, on the other hand, annihilates the tones 
throughout the whole output spectrum and provides a very smooth noise 
shaping. Simulation results are provided to support our theoretically derived 
results. Implementation issues of the irrational initial condition in the digital 
domain are also discussed and investigated together with the effect of finite 
accumulator size on the Noise Shaping Quality factor. 
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3.1 Introduction 

The exact analysis of AX modulators is not a straightforward task, 
because of the severe non-linearity introduced by the single-bit quantizer 
[49]-[58]. As explained in the previous chapter, the commonly employed 
approach that attempts to overcome this difficulty is the linear 
approximation method, which involves modeling the quantizer error as an 
additive white noise source. Despite the fact that the linearized model can 
provide reasonable results with respect to SNR performance, it completely 
fails to determine the exact spectral shaping of the modulator, which is 
highly likely discrete and colored. 

A more advanced linearized model for analyzing AX modulators is the 
so-called describing function method, in which the quantizer is modeled as a 
linear gain (chosen in a minimum mean-square-error sense) followed by an 
input dependent additive white noise source [59]. This improved model 
gives crude explanations on the input dependent stability and large 
amplitude limit cycle behavior of the modulators. However, it is still an 
approximation in the sense that, except the fundamental component 
(assuming a sinusoidal input), all other harmonics generated by the quantizer 
are perfectly removed by linear filters within the modulator [15], [57]. The 
exact quantization noise shaping performance of AX modulators is an 
important characteristic, as in certain applications the spectral spikes 
contained in the output spectrum may be strictly objectionable, for instance 
in the Fractional-N PLL frequency synthesis applications [60], [61] or in the 
audio applications [7], [39]. 

Fractional-N PLL frequency synthesizers utilize an all-digital AX 
modulator to generate accurately defined frequencies and desired channel 
selection in wireless transceivers [62]-[64]. Figure 3.1 illustrates such a PLL 
synthesizer where the output of the AX modulator is used to modulate the 
instantaneous division ratio of a multi-modulus frequency divider. The 
division ratio alterations take place very rapidly and in a random/pseudo- 
random fashion so that the phase noise and spurious contents are pushed to 
higher offset frequencies, where the existing PLL loop filter can easily cut 
them off [64]. Consequently, the phase error (analogous to the quantizer 
error in data conversion applications) spectrum is shaped by the high-pass 
noise transfer function of the modulator, so that the phase noise in the 
vicinity of the carrier (analogous to the baseband frequency in data 
conversion applications) is small. This technique permits narrow step sizes 
compared to the reference frequency and fast settling times while improving 
the phase noise performance of the PLL frequency synthesizers [64], [65]. 
The Fractional-N PLL concept will be explained in greater detail in Chapter 
5. 
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Figure 3-1. Fractional-N PLL synthesizer using AX modulation 

The difficulty involved in the realization of Fractional-N frequency 
synthesizers is that the A£ modulator is driven by a worst case input, i.e. a 
constant input that is the desired fractionality. The synthesized frequency at 
the PLL output is thus given by [60]: 

/»=/^(A' + -y) (3.1) 

where N is the integer part of the division ratio of the multi-modulus 
divider and X € [0,1) , the input to the modulator, is a constant number 
indicating the desired fractionality. It is a well known fact that constant 
inputs are the worst case inputs to the oversampling A£ modulators as the 
successful operation of these modulators relies on the assumption that the 
successive input samples are uncorrelated (i.e., busy inputs), and under static 
input condition this is definitely violated [49]. These complications are 
worsened in the context of Fractional-N synthesis, because the modulator is 
implemented in the digital domain and due to the finite-word-length of the 
arithmetic, the input as well as all the internal states are represented by 
rational numbers. It is well known in the art that simple rational dc inputs to 
the modulator produce purely limit cycle oscillations and tones, and hence 
no noise shaping is present at all [25], [26]. 

Another crucial point here, is that, unlike data conversion applications, in 
Fractional-N synthesis it is very dangerous to generate out-of-band^ limit 
cycles and tones as these components are very likely to be folded back to the 

^ Note that in Fractional-N synthesis the baseband frequency is defined as the pass-band of 
the low-pass loop filter in the PLL. Out-of-band frequencies correspond to the range of 
frequencies starting from the edge of the pass-band cut-off frequency to the Nyquist 
frequency of the modulator. Similarly, an equivalent oversampling ratio is also defined as 
the sampling rate divided by two times the baseband frequency. 
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baseband frequencies, because of the non-linearity in the implementation of 
the phase detector [61], [66]. Consequently, suitable AZ modulators for 
Fractional-N synthesis should produce tone free output spectra throughout 
the whole of the spectrum, not only in the baseband [67]. This fact draws 
attention to the importance of obtaining a very smooth noise shaping (i.e., a 
noise shaping without discrete spectral spikes throughout the whole of the 
spectrum) from a AZ modulator when deployed in a Fractional-N synthesis 
application. This is the main motivation for the exact analysis of AZ 
modulators with rational dc inputs performed in this research monograph. 

As cited before, exact analysis of the first-order AZ modulator with dc 
inputs was performed by R. M. Gray in [49], where it was shown that the 
output spectrum consists of discrete spectral spikes whose locations and 
amplitudes are input amplitude dependent. The technique was then extended 
for the higher-order MASH [35] AZ modulators in [50] and [51], using 
concepts from the ergodic theory [68]. All of these exact analysis studies 
were based on obtaining a non-linear difference equation from which the 
asymptotic autocorrelation function of the binary quantizer error was 
determined. The results in [50] and [51] have mathematically proved that in 
a higher-order MASH modulator the white noise assumption for the binary 
quantizer error from the last stage is indeed correct, when irrational dc 
inputs are employed. The same approach performed for the exact analysis of 
MASH modulators are also used in [52] and [53] in order to analyze higher- 
order single-stage modulators with irrational dc inputs. 

All of the above-cited studies on the exact analysis of modulators assume 
that the initial conditions in the accumulators (for a digital modulator) are 
zero. In fact, for the case of irrational dc inputs the asymptotic behavior of 
the quantizer error is not affected by the initial conditions [15], [49]-[53j. 
However, when the input is a rational constant (as it is the case in Fractional- 
N synthesis), such initial conditions may be very effective in randomizing 
the binary quantizer error. It is, therefore, the aim of this chapter to provide 
exact analysis results for the higher-order MASH modulators with rational 
dc inputs and non-zero initial conditions. In our analysis, we follow similar 
steps to those performed in [51] and [53] and make use of the theorems from 
ergodic theory as well as the characteristic function method combined with 
the solution of the non-linear recursive difference equation. Our results show 
that an irrational initial condition imposed in the first accumulator 
guarantees a tone free output spectrum for a third or higher-order MASH AZ 
modulator when driven by rational dc inputs. 

The rest of this chapter is organized as follows. Section 3.2 presents the 
recursive non-linear difference equations for the first-order and higher-order 
MASH modulators. Section 3.3 involves the derivation of the 
autocorrelation function of the binary quantizer error using the non-linear 
difference equation obtained in the previous section. Simulation results, 
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which justify our theoretical developments, are presented in Section 3.4. 
Implementation issues for the realization of an irrational initial condition in 
the digital domain are considered in Section 3.5. Finally, conclusions are 
given in Section 3.6. 

3.2 Non-linear Difference Equations 

3.2.1 First-order AZ Modulator 



Consider the first order modulator as shown in Figure 3.2, which is 
described by the following equations: 



v[n] = 



|vo; « = o 



n 






(3.2) 



w[«] - A‘[«l- y[w]; n ~ 1,2, ••• 



(3.3) 



y[w] = C(v[«l) = 




ifv[n]>0 
if v[«]< 0 



(3.4) 



where Vq is the initial condition in the accumulator, ±b is the single-bit 
(or binary) quantizer output levels, and Q(-) is the quantization function. Let 
us define the binary quantization error sequence such that e[n] = y[n]-v[n]. 
We will now derive a non-linear difference equation, which relates the 
binary quantizer error samples to the input of the modulator. 



■*|nl 







Figure 3-2. Discrete-lime model for the first-order AI modulator 



Using (3.2) and (3.3) one can write 
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= = (3.5) 

Using the definition for the quantizer error, equation (3.5) can be written 
as v[n] = x[n~Y\-e[n-\\ and hence an alternative representation for the 
binary quantization error sequence is: 



e[n] = 



jc[w]- v[« + l]; rj = \,2,- 



(3.6) 



We define here a scaled and shifted version of the binary quantizer error 
such that: 



w[«]s n = \,2,-- 

2 2b 



(3.7) 



which can expressed as follows by using (3.6) 



, , 1 v[n] 

w[„] = - + -LJ. — i i- w = i,2,- 

2 2b 2b 



(3.8) 



Note that one can easily recover each sequence of our interest from w[n] 
as follows: 



2b 2 2b 2b 



(3.9) 



v[«] . T 1 x[n~\] 

2b 2 2b 



(3.10) 



y[n] _ v[n] ^ e[n] v[rt] ^ x[n] v[« + l] 

2h 2b 2b 2h 2b 2b 



(3.11) 



Using (3.10) we can write = h{« + 1 ]-- + ^^, substituting this 

2b 2 2b 

value and (3.10) into (3.11) gives: 



= + 

2b 2b 



(3.12) 



Property 3.1: If 4«]€ [-^,+6]; V« and € [-26,+2^) , then 
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v[«]e [-26,26); V« (3.13) 

Proof: The proof is by induction as follows: We know that by definition 
-26^V() <26 and assume that (3.13) holds for n-\, i.e. v[w-l]€ [-26,26) . As 
stated in (3.5) we also know that v[n] =x[n-l] -y[n -1] + v[n-l]. Now assume 
that v[w-l] >0=^ y[«-l] = +6 , and hence v[n] =x[n~l] - b + v[n -1] , which can 
be written as 

x[n - 1] - 6 < v[w] < x[fi - 1] - 6 + 26 (3.14) 

Since |x[«]| S 6; Vn we deduce that 

-26 = -6-6 <v[«]< 6-6 + 26 = 26 (3.15) 

For v[«-l]<0=>y[w- 1] = -6 and v[n] =x[n-\] + 6 + v[« -1], which follows 

x[M-l] + 6-26<v[«]ct[«-l] + 6 (3.16) 

-26 = -6 + 6 - 26 < v[h] < 6 + 6 = 26 (3.17) 

This property indicates that given both the input of the modulator and the 
initial condition imposed on the accumulator is bounded within [-6, +6] and 
[-26, +26), respectively, then the accumulator state v[n] will be bounded and 
be within two times of the quantization level. This ensures that under these 
conditions the quantizer will never be overloaded. 

Property 3.1' The shifted and scaled version of the error sequence w[n] 
satisfies 

0<w[/7]<l, (3.18) 

Proof: In order to prove this property, we first show that 

Vn and (3.18) directly follows by application of this result to 
(3.7). By using the fact that -26 < v[n] < 26; Vn one can write 

for v[«] > 0 =^> y[rt] = 6 , -6 = 6 - 26 < e[n] ^ 6 (3.1 9) 



for v[n] < 0 => y[n] = -6 , -6 < e[n] < -6 + 26 = 6 



(3,20) 
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thus we conclude that - — . 

2 2h 2 

After having established these properties, we rewrite (3.12) as 
-= w[n - 1] + - ~^2b ~ substituting (3.3) and (3.10) into this 

equation results in the following relation: 



w[n\ = w[n - IJ + 



x[w - 2] Q{2bw[n -]]- h + x\n -2]) 



(3.21) 



Putting the initial condition on w[n] 



ir 1+^ ^ = 1 

wf/}1 = < 

' ^ „=2,3, 

2h 2b 



(3.22) 



Property 3.3: w[n] can be expressed as: 



VV[/7] = / wIa? - 1]+-^+ I' ^ ~ 



(3.23) 



where the notation {w{n\) stands for the fractional part of w[n] (i.e., 

(^["l) “ M^lmodl ). 

Proof: From (3.22) and the definition for the single-bit quantizer output 
stated in (3.3), we can write 

. x[n-2] 1 1 x[n~2] 

w’[w-l] + — if w[n-\]>- — 

M«] = i « = 2,3,-(3.24) 

r ,1 x[n-2\ 1 , . 1 x[n-2] ^ ’ 

w[^-\]+ +- (/ w[n-\]<- — 

2b 2 2 2b 

First it is necessary to find the upper limit for the value 

+ + Using the well-known Cauchy-Schwarz inequality [22] 

2 2h 

and by making use of (3.18) we can write 



r 1 x[n-2\ , ,,| 1 x[n-2\ ,11., 

+ - + W[/7~1] + -+ < I+- + -<2 

2 2b I I 'I 2 2h 2 2 



(3.25) 
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2 rf/f 51 

First we assume that iv[w-l] + — + — ^<1. Since the left-hand side of 

this equation is non-negative (see Property 3.2) and its value is confined to 
be within one, its fractional part is directly itself. Therefore, 

w[»] = l^n - 1] + Y + ^ ^ when n = 2,3,--- 

Now assume that 1< w[w-l] + — + — ~<2. In this condition, w[n] is 

given by the first part of equation (3.24), that is: w[n] = vv[n-l] + ^t^ — 

2b 2 



Also observe 



that 



from (3.18), under this assumption 
0 < w[« - 1] + ^ < 1 . Adding “1” and taking the fractional part of the 

final value does not change the final result: 

w[w] - (^n - 1] + j ~ ^ ^ ^ ; w = 2,3,--- 

The proof is now completed. 

Theorem 3.I.- For a first-order AS modulator shown in Figure 3.2 with 
arbitrary input x[w]g [~b,+b]; V« and arbitrary initial condition in the 
accumulator v^, e [-2b,+2b ) , the normalized binary quantization error 
sequences is given by 



2* 2 \2b ^[2 2b} 



(3.26) 



Proof: In order to prove this theorem, first we must show that 

wf«l = (— + Vf— n for n = l,2,---. The summation is equal to zero 
\2b ^[2 2ft j/ ’ 

when the upper bound is negative (i.e., n = 1). The proof can be 



accomplished by the method of induction. Note that w[\] = . In modulo- 

one arithmetic we have (~) = l+~ if — <0. Using Property 3.3 stated in 

\2ft/ 2ft 2ft 

(3.23) one gets ‘ ^hat for all real 



numbers 



b, {(^a) + (ft)^ = {a + b) . Now assume that 
{n-l] = /-^+^|^-j+-^^j\ is true. According to (3.23) and using the above 



a and 

«-3 



equality we can write 
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vw[«] = 





2 2b 



+ — + 
2 



+ y f i + iLti l\ 
2^ / \2^ 2/. J/ 



(3.27) 



Equation (3.26) directly follows by inserting (3.27) into (3.9), which 
completes the proof. 

For a dc input x[n] = X for all n (i.e., the input is held at a constant for 
all time), hence the normalized binary quantizer error sequence can be 
written as: 



e[n] = 



1 

2 





X 

H 

2h 



n = 0,\,-- 



(3.28) 



As a matter of convenience, we follow the same notation as used in [49] 
and define 

j}s{]/2 + X/{2b)) (3.29) 



?[«]s(v„/(2A) + «/3) (3.30) 

and hence £[«] = l/2-f[«] . 



3.2.2 Higher-order MASH AZ Modulator 



A discrete-time model for the I"' -order MASH AZ modulator, which is 
formed out of first-order stages, is shown in Figure. 3.3. 

The input to the modulator goes into the first stage, and the negative of 
the binary quantizer error sequence from the first stage is fed-forward into 
the second stage as its input, and visa versa. The outputs from all stages are 
processed by an error cancellation network, which removes the binary 
quantizer error components from all the stages except the last one [7], [35], 
and [40]. The difference equations describing this Z,'* -order MASH 
modulator are given as follows: 



v,[n] = 



fv,,; n = 0 
|M,[rt-l] + v,[rt-l]; 



n = \X- 



;/ = l, 2, 



(3.31) 



u^[n] = n = 1,2,--; / = 1,2, •••,/. 



(3.32) 
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x[w] = 0(v,[w]) = 




v ,[«]>0 
< 0 






(3.33) 



w,.[n] = 



' = 1 ,y„ 



(3.34) 



^il^J ~ ;/ = 1,2, (3.35) 

where and y. {i-\,2,---,L) are the initial condition on the 
accumulator and the output from the i"' stage, respectively. 




Figure 3-3. Discrete-time model for an Z."' -order MASH AI modulator 

As it is well known in the art, for an ideal U' -order MASH modulator the 
only noise term that appears in the output is the binary quantizer error from 
the last stage [35], [40]. The reason for this is that, as mentioned before, the 
quantizer error from all stages except the last one is cancelled in the error 
cancellation network. Hence, all that is needed for an exact analysis of a 
higher-order MASH modulator is the behavior of the binary quantizer error 
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from the last stage. Since every stage of the MASH modulator is built out of 
first-order ones, we can use Theorem 3.1 as long as the initial condition on 
each of the accumulators is bounded within the range [-2b,+2h). Namely, if 
e [~2^),+26);/ = 1,2, , then the normalized binary quantizer error from 
each stage is given by 






lb 



2 

2 




lb 




lb 



/ = ] 

;Vff 

/■ = 2,3,-,I 



(3.36) 



Since x[n]e[-b,-\-h] for all n and v, ^ e [-26,+26) , (3.36) is valid for / = 1 
due to Theorem 3.1 and furthermore it implies that for all n . 

With this conclusion, if v^^e[-2b,-\-lb) then (3.36) is also valid for / = 2 . 

Similar reasoning can be applied by induction to show that it is valid for all 
stages [51]. The following theorem gives a general description for the binary 
quantizer error from the last stage. 

Theorem 3.2: For an Z,'* -order MASH modulator formed out of first- 
order stages with any arbitrary input [-6,+6] for all n and arbitrary 

initials conditions on the accumulators such that (-2Z>,+26); / = 1,2, , 
the normalized binary quantizer error from the last stage is given by 




Proof: The proof is by induction on L . Assume that L = 1 then (3.37) 
coincides with (3.36). Assume that (3.37) holds for L, and write £;,.,[«] 
using (3.36) 






(3.38) 



Inserting (3.37) into (3.38) we obtain 




Simple algebraic manipulation yields 
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which complies with (3.37). 



Note that (3.37) applies to any arbitrary input and initial conditions. In 
this research monograph, we consider only rational dc inputs in the analysis, 
as it is the case in Fractional-N PLL synthesis applications. Once again let 
x[n] = Xe [~b,+b] for all n , and using the definition for p stated in (3.29) we 
write 



(3.39) 

that eJn] = l/2-Ci[«], 



(3.40) 

and (3.39) relate the 
normalized error sequence to the input of the modulator. These equations 
will be used in the following section to determine the statistics of the binary 
quantizer error. 

3.3 Statistics of the Quantizer Error Sequence 

3.3.1 Preliminaries 

Our primary goal is to evaluate the sample mean, average power and 
auto-correlation of the binary quantizer error sequence. For a stationary 
ergodic random process i[«] (« = 0,1,--*), the sample mean, average power 
and auto-correlation are defined by 




For simplicity purposes we define such 
where 




The non-linear difference equations in (3.28) 






(3.41) 




(3.42) 
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1 . 

/•;. ( 0 = Jim — ^ x[n ]jr[« + /] 



(3.43) 



if the limits exist. Let us define another ergodic random process 
y[n\ = l/2-x[n] and observe that the statistics of y[n] can be recovered from 



those of i[«] as follows: 




M{y)=\H~ M{x) 


(3.44) 




(3.45) 


r-{l)^\lA-h4{x) + r-{l) 


(3.46) 



The following theorem, which states important facts on the distribution 
of the fractional numbers, is well-known in the ergodic theory and will play 
a fundamental role throughout our higher-order MASH analysis, as it did in 
[ 50 ]-[ 53 ], 

Theorem 3.3 (WeyLs Theorem [15]): Let c[n] = a„+a^n + a:n' +'■■ + a^n^ 
be a polynomial with real coefficients. If among at least one of 

them, except the constant term , is an irrational number, then for any 
Riemann integrable function /(•) the following equality holds. 

Jim— (3.47) 

r»=(l ,) 



For the purpose of later use, several useful properties of modulo- 1 
arithmetic are also given here. Let m be an integer and a and b real 
numbers. It follows from the definition of the fractional operator that 



{a±m) = {a) and i^{a) ± = {a± = {a±b') 


(3.48) 


(fl + 6) = (a) + (A}-j;,, ,,„((/.)) 


(3.49) 



where the indicator function is defined as 
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1 ify^[a,c) 
0 otherwise 



(3.50) 



with 0 < fl < c where c < 1 . 

3.3.2 First-order AZ Modulator 

Now assume that the input X to the modulator is a rational number, this 
implies that p is also rational. We make use of the following lemma 
adopted from [49]. 

Lemma 3.1 [49]: If /? = K/N where K and N are relatively prime 
integer numbers, then the collection of numbers = 

corresponds to the same collection of numbers = 0,l,"-W-l}. In other 

words, there is a function /mapping O = {n = 0,l,"-,A^-l} onto itself with the 
property that 



nK\ /(«) 



N / N 






(3.51) 



Observe that in this case ^[n] is periodic with in n (see 3.30), and the 
infinite summation in (3.41) reduces to a finite summation. Using this fact 
and the property in (3.49) the sample mean of ^[«] can be written as: 



NAf\2b N N^\lb \n 



(3.52) 



Assuming > 0 (recall that e [-2b, lb ) , and note when > 0 , 

(vo/(26)) = Vo/(26)) and using Lemma 3.1 we get 



A^(^)=2:L+-Ly2i__Ly;^ („/a^) 



(3.53) 



The above equation can be easily evaluated when the appropriate limits 
are used for the conditional summation (rightmost term) such that 



M(0=A+lviL_l y, 

2b N^N N , 



(3.54) 
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where |”a:'| is the smallest integer greater than or equal to . As a matter 



of convenience let us call r|o - 



2b 



. With this definition 



Substituting (3.55) into (3.44) we obtain an expression for the mean 
value of the normalized binary quantizer error sequence. 

= + (3.56) 

2N 2b N 



The second moment of ^[n] is given by 




Combining (3.57) with (3.49) and using Lemma 3.1 

1 [ 

M(0 = -I 



2h N '' 



(3.57) 



(3.58) 



By expanding the square term in (3.58) we get 



2h N V ' 2h^v-f^.u ' [lb] ’ 

(3.59) 



Each term in (3.59) can be evaluated independently as follows: 



_L vf-lt 

J ~ 



(3.60) 



l ^ 2 Vq N -\ 

~N N 

»»o 



(3.61) 
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(3.62) 



-Lf^y („/w)=^-Lyi=^fi-Ik 

AtZ- 26 ^ IbN^ 2b N 



2v. 1 2vn 



n=DO 



(3.63) 



i / \2 /• 

J_Y(^1 - ^2. 

Af 26 ) “ 2b 



(3.64) 



n =0 26 

By substituting (3.60)-(3.65) into (3.59), one can achieve 

= + + (3.66) 

6N^ 2b[N 2b N ) n[ N ) ^ ^ 

Hence the average power of the binary quantizer error is found using 
(3.45), (3.55) and (3.66) followed by some algebraic manipulations. 



12 26 26 N j N N TJo 



(3.67) 



In order to find the autocorrelation function of the binary quantizer error 
we write 



y(— +— 

'' N N 26 



InK ^ vA/ {n + l)K ^ Vq 
Nj^\N 2bl\ N 26 



(3.68) 



iv4^\iV 2b 
«»0 ' 



nK . v« 



IK 



N 2bl'^\N ' 



InK 




\N 


2b, 



(3.69) 



The last equation can be further simplified to 
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r, (/) = ^[«D + (/^ >M(f 

«=0 ^ * \ 



I nK Vq 

).u\ 



(3.70) 



In order to obtain a compact formula for the rightmost term in (3.70) we 
need to determine the set of specific integer numbers in = = 

satisfying the following inequality 



IK\ nK Vq, , . 

— ) + ( — + — )>1 ;weQ 
N \n 2h 



(3.71) 



which can also be written using Lemma 3.1 as follows: 

;„,n 

N \N 2b/ 



(3.72) 



In determining the specific integer numbers in LI which satisfies (3.72), 
a very useful lemma is introduced here. 

Lemma 3.2: If K and N are relatively prime positive integer numbers 
such that^** A'/7ye[0,l] and v„/( 2 ^))e [- 1 , 1 ) , then there are only and only /(/) 
different integer numbers in the set of Q = {0,l,---.A^-l} satisfying the 
inequality in (3.72). 

Proof: With analogy to Lemma 3.1, we note that 



N 2b 



N 



n-0X--,N-\ 



where T(n,vO is a real number in the set of [0,A^ ). In fact, when v^/2b is a 
rational number then the value / («,Vo) is an integer belonging to the set of 
U . However, the above argument should be included in order to cover the 
cases of irrational and rational initial condition values. The next step is to 
prove that for a specific initial condition each n number results in different 
^ values. Mathematically this is written, 

/'(«,Vo) I{k,vA , , 

and n*k 

N N 

In order to prove this statement, we assume that the contrary is true, i.e. 
assume that there are two different integers n and k with n<k<N(k = n+ j 
for some non-zero integer 7 < ) forwhich I {nyVO = T{k,vO . This implies 



^ ^ - b^K/N + v,/(26)J = ^ + ^ - IkK/N + v,/(2A)J 



Note that by definition P=K/Ne[0,\] since X €[-b,b]. 
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where [crj represents the greatest integer which is less than or equal to a 
so that {a) -a- [o J . Thus, 

\_kKjN + vj(2b)\- \nKjN + vJ{26)J= ^ 

The left-hand- side of the last equation is definitely an integer, say L, 
L K 

then — = — with j < N , violating the assumption that K and N are 

j R 

relatively prime. 

For simplicity purposes, let us define = observe 

that ; VneQ. By using the fact that each n number 

corresponds to a different I (n,vO value and hence different F{n), one can 
establish an ordered sequence such that 

-^<F(«o)<F(«,)<--<F(«^_,)<^ + l where e a ;;? = -1 . In 

order to proceed further we determine a lower bound for the difference of 
and to this effect we write 





^ Vo 


, _L 






N 2b 


j T 


L N 2b 


Calling the integer -\n^^K/N + v^/(2b)_ 




In^K/N + v 



. The numerator in the right-hand-side of this 



equation should be a positive integer (recall that F(Wp^,)>F(»^) ) and thus a 
lower bound for the difference of F(Wp+,)-F(/7p) can be given as 

After having established these properties, the proof of Lemma 3.2 
follows by the method of contradiction. First we assume that there are 
/(/) + M numbers in the set of Q , where M is a non-zero positive integer, 
satisfying the inequality in (3.72). Hence, one can write 

1 ^ F(rtff_y) ^ 



Writing the differences 

)i\/N 

F(n,_,)-F(n^_,)i\/N 
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Summing up these terms yields 






(/(/) + AY-l) 
N 



and 



finally we get 






N 



+ 1 



/(/) 



which does not comply with the fact that + 1 . 

N 

Assume this time that there are /(/) -M numbers in the set of Q , where 



M is again a non-zero positive integer, satisfying the inequality of equation 
(3.72). Hence, 

^ < F(«„ ) < f («, )<•■•< F(n , ) < 1 

Writing the differences in a similar manner and summing up, results in 
^(^w-nt)^M-\)-F(no)> ^ fj-oni which it is deduced that 



which leads to a contradiction due to the fact that \ y and 

hence the proof is completed. 

We expand (3.72) using (3.49) 



M + 2L 

N N 



+ -^-X 

2b i' 



{nlN)>\\ 






(3.73) 



This inequality is investigated in two different cases: 

Case 1; The first case involves the condition that /]□>/(/). Note that as 

n is traced from 0 to N- 1, the value of x vo is altered depending 

^ 2b 

on whether n is greater or smaller than r)^ . Now suppose that w < r/p , and 
note that in this case x ® » hence 

N N 2b 

This yields ws + l.--,7?o “!}• With these lower and upper 

limits on n there are /(/) different integer numbers in satisfying the 
inequality of (3.72). This means that when n>T]f^ there shouldn’t be any 
further numbers satisfying this inequality (see Lemma 3.2). In order to check 
this, now suppose that n>r]^ and hence x {n/N) = \ . 



N N 2b 10 y f 



This suggests an impossible situation, as n can not be greater than A^-1. 
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Case 2: The second case is the condition when </(/) . As in Case 1, 
we follow a similar procedure and assume w < r/o , thus 

+ — + — >l=*n>7]o- /(/) 

N N 2b »o V / 

As is a negative number, we conclude that we With 

these limits, there are only t]o different integer numbers in satisfying the 
inequality of (3.72). According to Lemma 3.2, when «>rjQ there should be 
/(/)-Tjo further different numbers. Setting and x vo («/A^) = 1 , we get 

N N 2b 

The last equation implies that we {a^ + -/(/), + + -l}. 

As a result the conditional summation in (3.70) can be written as follows: 



2bl\<-{!^y{\ N 26/ J 

n=%-I(t) ' I 



2bnU^)M\N 2b 



I . Qqz! I .jr . \ 



(3.74) 



N^\N 2b N 2b 

/t=0 ' ' n=N*f}(j-UI) ' 



After some algebra, the summations in the right-hand-side of (3.74) can 
be evaluated as: 



1 ^ 

Zd \ N 1h 



M \ N 2b 



/(/)(2t;,- 7(0-1) ^ /(/) Vq 



(3.75) 



N Ad\ 2b N ^ \ N 2h 





1 

Vo 


u 


2b, 



/(/)(2t?„- 7(0-1) ^ 7(0 v„ 
2N^ N 2b 



(3.76) 



The above results show that the conditional summation in the left-hand- 
side of (3.74) is in fact a one-valued function. Noting that 
[l^) = {lKjN) = I(l)lN and inserting (3.75) and (3.76) into (3.74) we get 



NZd\u 2h [U‘4 :)m\n 2h ' ' 



nK Vn 



u. m I 
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Combining (3.77), (3.70), (3.66) and (3.55), we obtain 

M-il 

2 2 

(3.78) 

which, in conjunction with (3.46) and (3.55), gives an expression for the 
auto-correlation function of the binary quantizer error sequence that is given 
by 



^.(0 = 



-4^^+3^ + l ^ 

6N^ A/ 



■ + — -1 


I--1 


-!k' 


'-^0 


2h 


Nj 


N 

V 





+ 1 +l + (/;3) 





^ ^^0 1 2 L 


Lik 




12 6N^ 2b\ 


, N 2b N 

K / 


1 N 


N 



+ 2 - 



(A^ + 1) 
Ho 



.M 

2 



(3.79) 



For zero initial condition (i.e., =0 and = n ), the formulas given in 

(3.56), (3.67) and (3.79) reduce to M(e) = — , A/(e^) = — + and 

2A^ 12 

= which are identical to the equations derived in 

[49, eqs. 3.11a, 3.11b and 3.14], respectively. Note that in the analysis 
performed here, it was assumed that v<,>0. Similar procedures, however, 
were undertaken to extend these results for the negative initial condition by 
using the fact that (v^/(26)) = l-|vo/2/)| when Voe[-2^>,0). The details are 

omitted here for the sake of brevity. Nevertheless, the results are the same 
for negative initial conditions, and consequently the formulas in (3.56), 
(3.67) and (3.79) are still valid for this case. This fact will later be justified 
by computer simulations. It should also be emphasized that the above 
development does not put any restriction on the nature of the initial 
condition, hence the results are valid for both rational and irrational initial 
conditions. 



3.3.3 Higher-order MASH AX Modulator 

In order to evaluate the moments of the normalized binary quantizer error 
sequence in a MASH modulator, one needs to determine a closed form for 
the result of nested sums in (3.40). In achieving this, we follow similar 
notations as used in [51] and [53]. Let P^(«) represents “genetically” a 
polynomial of degree kin n . The qualifier “genetically” indicates that P^(n) 
represents any polynomial with degree k [51]. Effectively, one can write for 
any positive integer numbers k and m with m<k 
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= + (3.80) 

where a and b are two real numbers with a ^ 0 . This statement implies 
that the superposition of two polynomials of degree k and m with m^k y 
respectively, results in another polynomial of degree k . The reason for this 
notation that should be incorporated in the analysis of higher-order 
modulators is to pave the way for a streamlined solution to the problem 
without necessitating polynomial expansions. One only needs to know the 
form of the polynomials to derive the desired results. Similarly, we define 
and 3*(«) to represent generically polynomials of degree k with all 
rational and irrational coefficients, respectively. Note that (3.80) is also valid 
for rational and irrational polynomials. We also note that for any real 
numbers a and b, where , and positive integer numbers k and m with 

m^k y 

3((n) = a3*(n) + 6P„(n) (3.81) 

This fact implies that the superposition of an irrational polynomial of 
degree k and a generic polynomial (rational or irrational) of degree m 
where m<k , gives an irrational polynomial of degree k . This easily follows 
from the fact that addition of two numbers results in an irrational number if 
at least one of the numbers is irrational. 

Lemma 3.3: The Z/'* -order sum of the constant of “1” is a polynomial of 
degree L with all rational coefficients. Mathematically, this is given by 

(3.82) 

* j»0 * 1=0 



Proof: The proof is by induction on L. For L = 1 we have 

/ i-i 

= « = . Assume that (2.82) holds for L, i.e., 



*1=0 



n-| ijM h-\ 



ki=<i * 2*0 

where each coefficient c,. (/*= 0,1, 2, •••,!) is a rational number. Then, 



n-l */4.1-l ^-1 ^-1 n-l 



(3.83) 



* i + t =0 * 2=0 * 1=0 * 1 + 1=0 



We note a very useful identity found in [22] 
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^ /.+i 

=c,,-^ + 9t,(«) = 9?,,„(/r,„),Vi>0 (3.84) 



Substituting (3.84) into (3.83) proves the lemma. Some useful 
corollaries, which are a direct consequence of Lemma 3, follow as below. 

Corollary 3.1: The iT -order sum of a rational or irrational constant is a 
polynomial of degree L with all coefficients rational or irrational, 
respectively. This can be mathematically stated as follows: 



il =0 * 2=0 * 1=0 



(3.85) 



where is replaced with 51, («) and 3,(w) when a is a rational or 

irrational number, respectively. The proof follows as a direct consequence of 
Lemma 3 combined with simple algebraic manipulations. 

Corollary 3.2: A closed form for the L"' -order nested summation in 
(3.40) is given by: 







2h 




A 



(3.86) 



Proof: From the linearity of summation we can write 



2b ^ 



*/.=o 



2b 



t 



— =— +-^y i+i^y ^i+...+ 



*^I.O 

2b 



*>=*' *i.=0t7..|=C ti=l)i,=a 

An application of (3.82) to the above equation gives (3.86) proving the 
corollary. 

Using (3.86), the mean value of the sequence is given by 



lim— y + + 



lb 



(3.87) 



Throughout the developments in this subsection, we consider the case 
where the constant input is a rational fraction of the quantizer step size (i.e., 
)3 is rational) and the initial condition on the first accumulator is an 
irrational number ( v, „ is irrational). With this irrational initial condition 
assumption the infinite limit in (3.87) is evaluated with the aid of Wely’s 
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Theorem as stated in (3.47). Applying (3.80) and (3.81) to (3.87) and using 
(3.47) one can write 

A'-l * 

rt=0 0 

provided 1^2 (note that the Wely’s Theorem can only be utilized when 
L > 2 , i.e. when the order of the irrational polynomial is at least one), 
Similarly, the second moment of is written as; 

= i (3.89) 

n=0 0 

if L'tl. The sample mean and average power of the binary quantizer 
error from the last stage directly follows from (3.44) and (3.45), 
respectively: 



M(e,) = 0, = ^ 



(3.90) 



The following subsections will derive the autocorrelation of the binary 
quantizer error sequence for second-order and higher-order MASH 
modulators. Although the underlying ideas are the same, the development of 
the autocorrelation of the binary quantizer error in a second-order MASH 
modulator needs some special treatments, therefore the second-order case is 
separated here from the third or higher-order ( I ^ 3 ) cases. 



3.3.3.1 Second-order MASH: Using (3.40) and (3.43) the autocorrelation 
of^ 2 [«] is given by 



r^M)= -■^n)l^+^{n + l) + ^in + lY + 

' N-.-N^\lh 2h 2 2 \2b 2b^ ’ 2^ 2^ ’ 



(3.91) 



In order to be able to evaluate the above limit, we make use of the 
characteristic function method as in [15]. Let us define g( 7 )s(y). Note that 

this function is a well-behaved periodic function and it can be expanded into 
a Fourier series as 



g(r) = 



(3.92) 
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where the Fourier series coefficients are given by [12] 



i(0 = 



/ = 0 

2 



(3.93) 



Ini 



i*Q 



Substituting this Fourier series representation into (3.91) gives us an 
exponential expansion for the autocorrelation of l^^[n] as: 






' 2h 2h 1 



nsO 






(3.94) 



Changing the order of taking limits, then 



^ \ I- ^ ^ \ 2h 2h 2 2 I \ 2t) 2h 2 2 

'■c2(l^ = LLs(0g(k)hm~2,‘^ " 



(3.95) 



We continue by noting that for any integer i and any real number a we 
have the following equalities: 



gJ2ltia _ g>2s<(fl) _ 



(3.96) 



These equalities easily follow by the periodicity property of the 
exponential function with respect to 2;r . Now let us define 4>(/) as the 
sample average characteristic function of ^ 2 [«] [15], which is given by 



. N-[ 

0(i)= lim-i-Y 



;2w'' 



‘l.O . '’1,0 . /J 2 

— — - n*— ff • 

Ih U 1 



(3.97) 



where i is an integer. Similarly for two integer numbers i and k , we 
define thejoint characteristic function of as [51] 



1 7^®' 



>e ' 



(3.98) 
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Observe that ^(0 = <!>,(/, 0) for all integer numbers /. By using the 
periodicity property of the exponential function stated in (3.96), one can 
write 



Reordering (3.98) and using (3.96), we obtain 
0,(a)=lim— >e . 



( 3 . 99 ) 



( 3 . 100 ) 



Once again, an application ofWely’s Theorem to (3.100) results in 









0 , 






j + /t = 0 
otherwise 



( 3 . 101 ) 



The infinite summation in (3.101) is given by [22] 

lim ± y = I*’ = '"'‘'2^'' 

1 0 otherwise 

n=0 L 

Hence (3.101) can be re-written as follows: 



e '' i + k = 0 and pkl - integer (3.102) 

0, otherwise 

The limit interchange in (3.95) is justified in (3.102), which clearly 
implies that both 0(/) and <!>,(/, A) exist and lim<t,(i»0) approaches zero as i 

goes to infinity [15], [50] and [51]. The autocorrelation of fjL"] is found by 
substituting (3.102) into (3.99) 









^ >2/*/ 



* 1 . 0 *^ 



(Ni) 



*=- 

0 , 



= integer 
otherwise 



( 3 . 103 ) 



Using the Fourier series coefficients given in (3.93), we can write 
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^ 



> I 2o 2 



/ = 0 



/ ^ 0 anc/ f3kl = integer 



(3.104) 



Note that when / = 0 , then is already an integer number. Recall that 
'« 1 









-=1/12 [69]. The other summation with the exponential function can 



be simplified as follows: 

■ I - 



k~- 

(** 0 ) 



An k 



rl cod 



2;cfe/| 



2b 2 



(3.105) 



Hence, 






i _L 

4 '^ 12 ' 



— + — -t^osI iTtkt 
4 2n^4ie 

k=\ 



2b 2 



/ = 0 

(3.106) 

I ^ 0 and = integer 



In order to evaluate the above equation, one needs to determine the set of 
integer numbers k , which makes f5kl an integer. Alternatively, the question 
under investigation can be stated as; for a fixed rational number and an 
integer / ^ 0 what is the set of integer k numbers, which satisfies {^pkl) - 0 . 
Assume again that /? is given by a ratio of two integers numbers such that 
= K/N where K and N are relatively prime numbers with K <N . The set 
of integer k numbers which makes = 0 (where 1^0) is then expressed 
by = w;U(iV,/);w:l,2,-'-} in which ^{N,l) is a positive integer number which 
is a function of N and / , and is given by the following relationship: 



when = 0 and |/| ^ N 
NjGCDiN, |/|) , when (a^/|/|) ^ 0 and |/j < N 
1, when = 0 and |/[ > N 

NjGCDi\l\, N), when (j/j/A^) 0 and |/j > N 






(3.107) 
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where GCD(a,b) represents the Greatest Common Divisor of positive 
integer numbers a and b . The proof of (3.107) can simply be accomplished 
by using results from the elementary number theory [22], and therefore the 
details are omitted. From now on we will use // instead of n{NJ) for the 
sake of brevity. Using (3.107), we can write the following equality 



[ 126 2 'Jj [ ^1^26 2 ' ^jj 

for /^O and {j5kl) = 0. We note here a useful identity found in [69, eq. 



A360], 




(m9)=—{9-7cy — —7C^ 
^ 4^ ^12 



(3.108) 



where 0<$<27t . Let us call y s 2k{ 



we write 



|)>then using (3.108) 



-ir-y —7-YCos(my)= — wO 

2?r 27rVl,4 12 J 

Finally, the autocorrelation of ^j[«] can be found as 



(3.109) 






1 1 
— I — , 
4 12 



/ = 0 



J 1 / V 1 2) ; rx 

4 2?rV(4^ 12 / 



(3.110) 



From the above equation the autocorrelation of the binary quantizer error 
from the second stage is found in conjunction with (3.46) and (3.88) as 






12’ 



/ = 0 



2n 



^^1^4^^^ ^ 12 ; 



(3.111) 



3.3.3.2 Third or higher-order (Z, > 3 ) MASH: The autocorrelation of 
fiN is given by 



1 

= jlim-^ ,(«)+9t,(«))(3._,(« + 0+9?,('^ + 0) 



(3.112) 
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We again utilize the characteristic function method [15] as discussed in 
the previous subsection. Using the Fourier series representation for the 
modulo 1 operator, the autocorrelation of [n] is given by 






(3.113) 



where thejoint characteristic function is as follows: 









(3.114) 



The above equation can be evaluated using the following lemma. 
Lemma 3.4: Consider an irrational polynomial of degree(L- 1) 

in n . For all integer numbers i , k and I , the following equality holds; 






3/.,(«), if i + k^0,V( 

ifi + k = 0,kl^0 
0, if i + k =0,k/ = 0 



(3.115) 



Proof: Let Z,_i{n) = a,__^n‘-'' +ai_.n'-'‘ +-- + a^ where all the 

coefficients (/ = 0,1,2, ) are irrational numbers. Following this, we 

may write ^,_^{n + l) = a/_Xn + iy'~' + +a,_i{n + lY~^ + --- + a^. Each term 

of this polynomial can be expanded using the Binomial Theorem [22], which 
is given as follows. 



{a + hf = 




(3.116) 



where a and b are two real numbers and m is a positive integer. The 
binomial coefficients are given by: 
w'|_ ml 
k ~ {m-k)\k\ 

Using (3.116), one can write 



(/ + k)'^^a^n'" + ^ 4. kla/_2 ^ 



L~2 
m + \ 






(3.117) 
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Equation (3.117) clearly indicates that when i + k*0 the result is an 
irrational polynomial of degree (L-1) in n . On the other hand, / + ^ = 0 and 
kl^O results in an irrational polynomial of degree (L-2) in n . Finally, if 
/ + ^ = 0 and kl = 0 the result is zero. Hence, the proof is completed. 

Once again by making use of Wely’s Theorem and Lemma 3.4 we write 
(3.114) as 






1, i + k-Q and kl = 0 
0, otherwise 



(3.118) 



Note that in the above development, in order to be able to utilize Wely’s 
Theorem the order of the MASH modulator must be at least three (i.e., 
1^3). When L = 2 (3.115) is no longer valid, thus this case needs special 
treatment as was performed in the previous subsection. Also note that 
(3.118) implies that both 0(;) and 0/(/,A:) exist and limO,(/,0) approaches 



zero as i tends to infinity. So the limit interchange is again justified [15], 
[50] and [51]. The autocorrelation of ^^[«] is found by substituting (3.118) 



and (3.93) into (3.113) as 






4 2itk 2jz{-~k) 

{k*0) 

1 



/ = 0 



1^0 



By making use 






= 1/12 , we finally obtain 



'ft 



(0 = 



; = o 

4 12 

/;£0 

4 



(3.119) 



In order to find the autocorrelation of the binary quantizer error from the 
last stage, we combine equation (3.119) with (3.46) and (3.88), and the result 
is: 







/ = 0 
/;fi0 



(3.120) 
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3.4 Simulation Results 

In order to verify our theoretical results, we have undertaken extensive 
simulations on AL modulators employing single-bit quantization with output 
levels ± 1. A time-domain difference-equation model was run for 2'* number 
of output samples from which the sample mean, average power and 
autocorrelation of the quantizer error were obtained. Figures 3.4-3.T show 
the simulated autocorrelation function of the binary quantizer error along 
with its analytically derived/obtained result, for various modulators, inputs, 
and initial conditions. The simulated mean value and the average power are 
also listed in these plots. 

As seen from these plots, the theory developed in this research 
monograph and the corresponding simulations are in exact agreement 
justifying the derived formulas. Figure 3.4(a) and (b) illustrate the 
autocorrelation function of the binary quantizer error in the first-order AS 
modulator with a positive and negative initial condition, respectively, where 
X = 0.1. This result supports the fact that the derived formulas in (3.56), 
(3.67) and (3.79) are valid for both of the polarities of the initial condition. 
Figure 3.5 and 3.6 show the autocorrelation function in the second-order 
MASH modulator for various input amplitude and initial condition 
combinations. These figures justify the validity of the derived formulas in 
(3.88) and (3.111). As clearly shown in these figures, the autocorrelation 
function does now vanish for non-zero lag values, which is an indication of 
strong tone components in the output spectrum. On the other hand. Figure 
3.7(a) and (b) depict that the binary quantizer error sequence from the last 
stage of the third and fourth-order MASH modulators, respectively, has 
identical first and second-order statistics corresponding to those of uniformly 
distributed white noise, as predicted by (3.88) and (3.120). 
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mean value=-0.021 57870, average powersO.08359064 




^06 - 

-100 -50 0 50 100 

lag index (I) 
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Figure 3-4. Autocorrelation function for the binary quantizer error in the first-order modulator 
with 0.1 dc input: (a) v„ = 0.60684258354180 (b) = -0.60684258354180 
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Figure 3-5. Autocorrelation function for the binary quantizer error in the second-order MASH 
modulator: (a) A' = 0.I, v^^^=0.60684258354180 and V 2 ,,= 0 ,(b) .^ = 0.1, 

V, ^=0.60684258354180 and v^^=0.2 
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Figure 3-6. Autocorrelation function for the binary quantizer error in the second-order MASH 
modulator: (a) X = -0.i, v,_„ = 0.95012928514718 and „ = 0 , (b) X = -0.5 , 

V, ,^^-0.23113851357429 and =0.2 
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Figure 3-7. Autocorrelation function for the binary quantizer error in the third and fourth- 
order MASH modulators: (a) A' = 0.1, v^ = 0.60684258354 180 , Vi„-0 and (b) 

= -0.8 , V, 0 - -0.48598246870930 , v,,o = 0 , = 0 and 1 - 4,0 = 0 
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Figure i.8. The difference between the simulated and theoretically calculated autocorrelation 
function; (a) first-order modulator of Figure 3.4(a), (b) second-order MASH modulator of 

Figure 3.5(a) 
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Figure 3-9. Output noise spectrum for the second-order MASH modulator where X = 0. \ , 
V; 0 = 0 : (a) V, 0 = , (b) v, „ = 0.60684258334 J 80 
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(b) 

Figure 3-W. Output noise spectrum for the third-order MASH modulator where A' = 0. 1 , 
Vj,o = 0 and v, „ = 0 : (a) v, „ = 0 , (b) v, ^ = 0.60684258354180 

It should be also noted that only -100 and +100 lag indexes are shown in 
Figures 3.4-3.7 for clarity. For the sake of better accuracy, Figure 3.8(a) 
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demonstrates the difference between the theoretically derived and simulated 
autocorrelation function of the binary quantizer error in the first-order 
modulator for the case shown in Figure 3.4(a). This very small (on the order 
of le-7) discrepancy between the simulated and the analytically calculated 
autocorrelation may be attributable to possible numerical inaccuracies and/or 
drifts in the simulation software environment [67]. The difference between 
the analytically calculated and simulated autocorrelation function in the 
second-order MASH modulator for the case of Figure 3.5(a) is shown in 
Figure 3.8(b). This difference (on the order of le-5) is mainly a result of the 
finite number of samples that can be obtained from a computer simulation. 
Referring to (3.52), (3.57) and (3.68), the statistics of the quantizer error 
sequence in a first-order modulator are calculated over a finite number of 
samples. Therefore, unlike higher-order cases, the effect of the finite number 
of samples does not apply in this case [67]. 

3.5 Output Spectrum 

The output noise spectra (dc input removed) for the second and third- 
order MASH modulator with zero and non-zero irrational initial condition on 
the first accumulator are shown in Figure 3.9 and Figure 3.10, respectively. 
As shown in Figure 3.9(a), the second-order modulator exhibits strong tone 
components in the output noise spectrum when zero initial condition is used. 
Although these tones are substantially reduced when the initial condition on 
the first accumulator is set to an irrational number (see Figure 3.9(b)), there 
are still some periodic components in the second-order modulator as 
suggested by our theoretical findings. The effect of the irrational initial 
condition for the third-order modulator is made abundantly clear in Figures 
3.10(a) and (b), where big differences in tonal components exist between 
these two spectra. 

Throughout this research monograph we have used the Welch’s method 
[70] of averaging modified periodograms in power spectral density 
estimation. Specifically, 16 50% overlapping periodograms, each with 32K 
of length, were averaged before the Fourier transform was taken to produce 
an estimate of the power spectrum. Also note that each periodogram was 
processed by a Hanning window to reduce the spectral masking caused by 
the finite number of samples obtained from a digital simulation. The 
reduction in absolute power that is a result of the windowing was 
compensated by scaling the signal with the inverse of the norm of the 
window, which was 8/3 for the Hanning window [19]. The overlapping 
averaging involved in the Welch’s method is particularly advantageous for 
further reducing the spectral leakage while maintaining a good spectral 
resolution as well as providing a consistent estimate of the power spectrum 
[19]. 
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3.6 Digital Realization of the Irrational Initial Condition 

In the previous sections we have mathematically shown that the third or 
higher-order MASH AZ modulators provide tone free output spectra, when 
an irrational initial condition is used on the first accumulator. However, as 
mentioned before in Fractional-N frequency synthesis applications the 
modulator is implemented in the digital domain, and one might rightly 
consider how the irrational initial condition could be realized in the digital 
domain. The answer for this question relies on the fact that a rational 
number, which is represented by a ratio of two relatively prime integer 
numbers, will tend to an irrational number when the number in the 
denominator is made very large [49]. This suggests that provided the 
dynamic range of the accumulators is sufficiently large an odd^^ number 
initial condition can still act as an irrational number. 

Another point, which should also be considered here, is that the 
theoretical developments performed in the previous section for the mean 
value of the binary quantizer error sequence is for infinite time. This means 
that if the data was infinitely long the mean value of the binary quantizer 
error would approach zero (i.e., the mean value of the output would 
approach the value of the constant input). However, the Fractional-N PLL 
frequency synthesizer should settle to the desired frequency within a 
specified finite time (which is the settling time of the synthesizer). This 
implies that the number of samples that can be taken from the output of the 
MASH modulator is large, but limited. Arbitrarily choosing a large odd 
number for the initial condition on the first accumulator, therefore, may 
result in relatively large error between the mean value of the output and the 
constant input. This yields, of course, a relatively large frequency error in 
the synthesized PLL output. In order to minimize these frequency deviations 
it is desirable to utilize the minimum odd number, which is “1”, in the first 
accumulator as the initial condition. This is accomplished by setting “1” 
LSB in the first accumulator each time the circuit is reset [71]. A digital 
implementation of the 1-1-1 MASH modulator with the “1” LSB initial 
condition in the first accumulator is depicted in Figure 3.11. 

Figure 3.12 presents the output spectra of the three, four and five-stage 
MASH AX modulators, where the dynamic range of accumulators was set to 
20-bits, and “1” LSB irrational initial condition is used on the first 
accumulator. These spectra were obtained from the fixed-point simulations, 
where 16 50% overlapping modified periodograms (Welch’s method) were 
used to estimate the power spectral density. The results clearly indicate the 



11 



Note that here the reason for choosing an odd number is to represent the initial condition by 
a ratio of two relatively prime integer numbers, with the largest possible denominator. 




72 



ANALYSIS OF MASH DELTA-SIGMA MODULATORS 



effect of “1” LSB initial condition, as no tones are observed throughout the 
whole spectrum. 




Figure 3~I1. Digital implementation of the third-order MASH AI modulator 

In order to investigate the effect of “1” LSB initial condition for various 
accumulation size, here we define a new measure dubbed as Noise Shaping 
Quality {NSQ) factor. The NSQ factor of an Z,'* -order modulator is defined 
as follows: 

NSQ = \~~ ( 3 . 121 ) 

where SAD is defined as “the Sum of the Absolute Differences in power 
between the measured output quantization noise and the theoretical noise 
shaping performance of the modulator”, and TP is “the Total Power in the 
noise transfer function”. Mathematically, this is written as: 

fs 

NSQ = 1 --2 

|i^.(/)r df 

0 



( 3 . 122 ) 
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where /, is the sampling frequency of the modulator, 
is the theoretical noise shaping transfer function, and 
?£.(/) is the power spectral density of the quantizer error sequence. The 
NSQ factor gives an indication of how closely distributed the output 
quantization noise spectrum is with respect to the theoretical noise shaping 
performance of the modulator. 

It is clear, however, from the definition stated in (3.122) that for a white 
noise quantizer error the value of NSQ does not equal to unity (in fact it is 
much smaller than unity). This means that the maximum achievable NSQ 
factor for a real modulator is limited to that of the truly white noise. 
Therefore, before the extensive fixed-point simulations are undertaken, it is 
required that the NSQ factor of a white noise is determined and all the raw 
NSQ measurements are normalized to this value. 

For an ideal white noise quantizer error the power spectral density is 
constant over all defined frequencies. That is PEif)-o] for 0</</^, 
Therefore, 

|i -/’*(/)! ||^^e(/)|V 

NSQ„, = 1 ^ 2 

0 

We assume that the maximum output level of the quantizer is normalized 
to 1, so that for a uniformly distributed quantizer error the variance is 
smaller than 1. Hence, NSQ^,„=al^ where NSQ^ is defined as the Noise 
Shaping Quality factor of a AS modulator when the quantizer error sequence 
is assumed to be a sample random process of an ideal white noise. 

Figure 3.13 shows the spectral density for an artificially generated white 
noise, which is shaped by the third-order noise shaping transfer function. 
The theoretical noise transfer function is also overlaid on the same plot. The 
white noise^^ was uniformly distributed between -1 and 0, and has a mean 
value of-1/2 and a variance of 1/12. The NSQ calculation indeed shows that 
= 1/12 for this case, as predicted. 

The effect of the finite accumulation size on the NSQ factor of the three, 
four and five-stage modulators is illustrated in Figure 3.14. Note that in each 
of these fixed-point simulations, “1” LSB of the first accumulator was set at 



Note that in the digital implementation of MASH modulators there is no physical quantizer 
(see Figure 3.11). The quantization process is embedded in the overflow mechanism ofthe 
accumulators, which produces 0 or 1. The quantizer error is therefore bounded in the range 
[- 1 , 0 ]. 
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the beginning of the modulator start-up. Panel (a) of this figure was obtained 
by averaging the NSQ factors over the worst ten inputs (i.e., X was traced 
from 0 to 1 with 0.1 increments), while the NSQ simulations in panel (b) 
were averaged over all the specified inputs ( X was traced from 0 to 1 with 
0.01 increments) in a Fractional-N synthesis application^^ . 

The relatively low value of the NSQ at low number of bits is a result of 
the huge truncation error in the internal states. Four and five-stage 
modulators, however, demonstrate almost perfect operation (i.e., they 
produce a quantizer error which is really close to white noise), when the 
number of bits in the implementation is higher than 15. As shown in Figure 
3.14(a), in the three-stage modulator, as the accumulator dynamic range is 
increased the NSQ factor begins to deteriorate. The reason for this is that 
when the number of bits is large the third accumulator does not experience 
sufficient activity near the LSB positions, resulting in the tonal/structured 
output spectra. This deterioration is not clearly seen in Figure 3.14(b), 
because when the LSB of input word is set it behaves as an irrational 
number, and therefore produces tone free noise shaping regardless of the 
initial conditions [67]. In other words, like our irrational initial condition 
principle, the LSB set in the dc input word results in the proper operation. 
When the input amplitude approximates to an irrational number, there is no 
need to use an initial condition. However, the use of the initial condition 
ensures proper operation for all input amplitudes and should therefore be 
included in designing AS modulators for Fractional-N synthesis applications. 
This issue will further be elaborated in Chapter 5. Nevertheless, we conclude 
that four and higher-order modulators, on the other hand, are more robust 
with respect to the digital implementation of the “1” LSB irrational initial 
condition. 

3.7 Conclusion 

In this chapter, an exact mathematical analysis has been presented to 
study the effect of the initial conditions in eliminating the tones and 
oscillations in MASH AX modulators driven with rational dc inputs. The 
analytical solutions have shown that for the case of the first-order modulator, 
the quantizer error sequence contains periodic structures regardless of the 
choice of the initial condition. In the second-order case, the quantizer error 
sequence is partially decorrelated, but it still fails to comply with the white 
noise model. On the other hand, higher-order modulators may benefit from 
more randomization of the quantizer error. The more the accumulation 
stages, the more the quantizer error sequence is randomized. This chapter 



We assumed that f^j = IS^MHz and the desired frequency step was = 200^//z . 
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has mathematically as well as experimentally proved this argument. We 
have shown through a rigorous mathematical analysis that an irrational 
initial condition on the first accumulator produces a signal-independent 
uniform distribution white noise quantizer error sequence in a third or 
higher-order modulator. Additionally, the quantizer error sequence is 
uncorrelated with the dc input. To see this, we write the cross-correlation 
between the input and the quantizer error sequence in an V‘' -order MASH 
modulator as (/) = M(^ft)=XiVf(e^). Since the mean value of the quantizer 

error sequence is zero (see 3.90), it is concluded that r,^^(/) = 0 for all lag 
indices /. To support our theory, several simulation results were also 
presented. Specifically, the first and second-order moments of the quantizer 
error sequence in first and higher-order AZ modulators were simulated and 
compared against their corresponding analytic solutions. For each case, the 
results of our theory and simulations were in exact conformance, which fully 
justify our theoretical framework. 

One difficulty in the irrational initial condition principle is that it cannot 
be perfectly represented in the digital domain. The (approximate) realization 
of the irrationality in the digital domain has also been thoroughly 
investigated in this chapter. We have discussed that a rational number with a 
reasonably large denominator can behave like an irrational number in fixed- 
point implementations, if the numerator is a relatively prime integer with 
respect to the denominator [49]. While this is true in general, in the 
Fractional-N synthesis applications, however, the numerator cannot be 
chosen arbitrarily. In order to minimize the frequency settling errors at the 
synthesizer output, which are due to the finite length of the modulator 
output, we have opted to set “1” LSB on the first accumulator as an 
irrational initial condition. A new measure, NSQ factor, has been developed 
to investigate the effect of this irrationality assumption involved in the “1” 
LSB initial condition for varying accumulator sizes. Several fixed-point 
simulations were run for third, fourth and fifth-order MASH modulators 
with “1” LSB initial condition. The measured factors have shown that 
in the third-order modulator increasing the accumulator size more than 23 
bits lessens the irrationality assumption involved in the “1” LSB initial 
condition. This is because as the denominator (i.e., the accumulator size) 
gets larger the initial condition approaches zero. Fourth or higher-order 
modulators are more robust in terms of the irrationality assumption, because 
of the presence of more accumulation stages. These findings and results will 
be explained and exploited in greater detail in Chapter 5, which presents the 
details of the AZ controlled Fractional-N PLL frequency synthesizers. 
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(b) 
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(c) 

Figure 3-12. Output quantization noise spectra for MASH modulators from 20-bit fixed-point 
simulation: (a) three-stage, (b) four-stage, (c) five-stage 




Figure 3-13. Spectrum of an artificially generated white noise shaped by the third-order noise 
transfer function. Note that RSQ^„ « 1/12 
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Figure 3-14. Plot of the NSQ factors of third, fourth, and fifth-order MASH AS modulators 
(with the “1” LSB on the first accumulator set) for different number of accumulation size: (a) 
averaged over worst ten dc inputs, (b) averaged over all specified dc inputs 
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4. ANALYSIS OF SINGLE-STAGE DELTA-SIGMA 
MODULATORS WITH DC INPUTS 

In this chapter, we provide rigorous analysis of single-stage higher-order 
AZ modulators with rational dc inputs. Higher-order modulators with pure 
differentiation noise transfer function can only be stabilized (in the sense we 
define later) utilizing a multi-bit quantizer. A non-linear difference equation 
is derived describing the generic Z,"" -order modulator, from which no- 
overload stability conditions are obtained. It is shown that the long-term 
statistics of the no-overload quantizer error sequence are the same as that of 
the identical distribution uniform white noise process, provided that the 
modulator is at least a third-order one and an irrational initial condition is 
imposed on the first accumulator. 

4.1 Motivation behind the Work 

Higher-order AZ modulators achieve higher resolution by pushing more 
quantization noise outside the band of interest. In the previous chapter, we 
provided an exact analysis and results for higher-order MASH type 
modulators. Another important class of modulators, called the single-stage 
or interpolative structures, employs only one quantizer within the system, 
and constitutes a vast majority of practical designs. In general, an Z,'* -order 
single-stage modulator is the logical extension of the first-order one, and 
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thus contains L cascaded accumulators which result in a noise transfer 
function of the form of a pure differentiation (i.e., (l-z'‘)^). 

Unfortunately, higher-order single-stage modulators are not stable when 
a single-bit quantizer is used, due to the high degree of non-linearity. In 
particular, if one uses a single-bit quantizer with output levels ±b , then the 
internal states within the modulator grow indefinitely and the quantizer 
becomes significantly overloaded, even if the input amplitude is kept in the 
range [-b,b]. The use of multi-bit quantization can overcome this stability 
problem in higher-order single-stage modulators. Stability is often described 
in the sense that the quantizer is not being overloaded^"^ . For a non-overload 
stable modulator, all the internal states (especially the quantizer input) is 
bounded such that the quantizer error amplitude does not become greater 
than half of the bin width. 

This chapter is devoted to the exact analysis of higher-order single-stage 
modulators, which employ multi-bit quantization. First, the non-linear 
difference equation describing a generic L"' -order modulator is derived, and 
using this relation the conditions that lead the quantizer to operate in the 
non-overload region are determined. Long-term statistics of the granular 
quantizer error is then theoretically obtained and thoroughly verified by 
means of computer simulations. 

4.2 Uniform Quantizer in the No-Overload Region 

A basic component common to all AL structures is the quantizer, be it 
either single-bit (binary) or multi-bit. The input-output transfer 
characteristics for two-bit and three-bit uniform quantizers are shown in 
Figure 4.1 and 4.2, respectively, where ±b represent the maximum output 
levels, V is the input to the quantizer, y = Q{v) is the output of the quantizer, 
and A is the bin width. Note that, A = 2b/3 and A = 2^/7 for a two-bit and 
three-bit uniform quantizer, respectively. The quantizer maps the amplitude 
of its input to a finite number of output values, which are represented by a 
digital code word. Clearly only input values, which are separated by at least 
A can be resolved to different and unique output values. 

For an M-bit quantizer, the input-output transfer characteristic can be 
generalized as: 



y = Qiv) = 



= -i)a/ 2; (2"-' -l)i<v 

/ , , (A:-1)a<v</:A; 

(*-1/2)A; i = (-2*'-'+2),-,(2*'''-l) 

-6 = (-2"' +i)a/2, (-2''“'+i)4>v 



( 4 . 1 ) 



The no-overload quantization error is also referred as granular quantization error. 




ANALYSIS OF SINGLE-STAGE DELTA-SIGMA MODULATORS 



81 



where A = 26/(2“^ -l) is the quantizer bin width. The no-overload region 
of the quantizer is defined as the range of input amplitudes, for which the 
absolute value of the quantizer error e = y-v does not exceed half of the bin 
width (i.e., |e| < A/2 ). Figure 4.3 and 4.4 depict the quantizer error for a two- 

bit and three-bit quantizer, respectively, as a function of the input amplitude. 
It is clear from these figures that for a two-bit and three-bit quantizer, the no- 
overload region is given by ve[-2A,2A] and v€ [-4A,4A], respectively. In 
general, if the input of an M-bit quantizer is bounded in the range 
ve [-2^“‘A,2 *^“'a], then the quantizer operates in the no-overload region and 
consequently |e| < A/2 . 



y=^v) 

4 





i 


k 


A/2 

-2A -A 1 










A 2, 

-N1 

/•= -3A/2 

► 


A ^ 




4 



no overload region 



Figure 4-1. Input>output transfer characteristics for a two-bit uniform quantizer 



y=Q(v) 




Figure 4-2. Input-output transfer characteristic for a three-bit uniform quantizer 
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e=y-v 




Figure 4-3. Quantization error of the two-bit quantizer 



e=y-v 




Figure 4-4. Quantization error of the three-bit quantizer 

Normalizing the quantizer output and input by the quantizer bin width A , 
and using (4.1), 



e V y 

-e = — = — — 
AAA 




(2"-' -l)<v/A 
{k-\)<v/A<k; 

(_2'^-'+i)>v/A 



0 



( 4 . 2 ) 



If the quantizer input satisfies ve [-2" ‘A, 2*” 'a] or v/Ag [-2^ ',2^ '], then 
the above formula simplifies to 



e 

A 



= |_(4_|)_1; (i-l)<v/A<*; 

k = {-2^'' 



-£ = 



( 4 . 3 ) 
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This can be further simplified using the fact that every real number r can 
be uniquely written in the form r = [rj+ (r ) , where [rj is the greatest integer 

less than or equal to r, and 0 ^ (r) < 1 is as before the fractional part of r (or 

rmodl).Then, 



-e = 



e 

A 




2 




2 



Finally, we obtain 



e 




(4.4) 



The above equation gives a mathematical relationship between the 
quantizer input and normalized error. Note that this relation is only valid in 
the no-overload region. In other words, if the amplitude of the quantizer 
input is confined in the range v€ [-2^“'A,2 ^"'a], then it is guaranteed that the 
absolute amplitude of the quantizer error is within half of the bin width 
|e| < A/2 , and (4.4) holds. 

4.3 Non-linear Difference Equations 

In this section, we derive non-linear difference equations governing the 
higher-order AS modulators. We start our analysis with the simplest form of 
the higher-order single-stage stage modulator, which is the second-order 
system, and inductively apply the results to obtain a generic formula. 

4.3.1 Second-order AS Modulator 

First, consider the second-order modulator shown in Figure 4.5 
employing an M-bit quantizer having maximum output levels ±b . 



M-h'w 




Figure 4-5. Second-order single-stage AI modulator 
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The modulator input, output and accumulator outputs are represented as 
time sequences x[ti\, y[n\ Vi[«] and V 2 [«], respectively. We assume that the 
system is initiated at time n = 0, or in other words, the input to the system is 
applied at time zero. The initial conditions on the first and second 
accumulators are defined as v,[-l]=Vi_, and Vj[0] = v 2 o, respectively. Note 

that since the first accumulator is non-delaying, we impose an initial 
condition at n = -1, in comparison to the second accumulator on which the 
initial condition is specified at n = 0. The equations describing the second- 
order modulator are as given below: 



«,[«] = 4«] - y[«]; « = o,i.2,--- 


(4.5) 


v,[«] = v,[« - l] + w,[/j]; n = 0,1,2,- •• 


(4.6) 


1'*] = ^1 [«] ~ ; « = 0,1,2, • • ■ 


(4.7) 


V'2[«] = V2[«-l] + W2[n-l]; « = 1,2,--- 


(4.8) 


y[n] = Q(vAn]}; « = 0,l,2,-- 


(4.9) 



where Q{-) is the quantization function of an M-bit quantizer, and is given 
by (4.1). 

We define the quantization error at any time n as follows: 
e[n] = y[fj]-V2[n]; « = 0,1,2,- ■■ (4.10) 

and with this definition, (4.5)-(4.9) could be summarized as 
V( [«] = V| [rt - 1] + jc[fl] - y[n] ; n = 0,1,2, ■ • ■ (4.11) 



Vjlw] - Vjlw- 1] + V|[w-2] + jc[ff-l]-2y[w- 1]; « - 1,2,- •• (4.12) 



>’t«] = Vj [/I - 1] + V, [n-2] + x[n - 1] - 2y[n - 1] -l- e[n ] ; n = 1,2, • • • 



(4.13) 
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Note that by definition stated in (4.10), the output atn = 0 is 
y[0] = Vj 0 + e[0] . Using the fact that -y[n - 1] = -Vj [« - 1] - e[n - 1] ; n- 1,2, - - • , the 



last equation can also be written as 

y[«] = [n ” 2] - Vj [n - 1] + x[n - 1] - 2e[n - 1] + e[n] ; n = 1,2, ■ (4.14) 

According to (4.14) aX n = I 

y[l] = - V 2.0 + 40 ] - 2e[0] + e[l] (4. 1 5) 

and for the rest of the time (4.14) is still valid 

y[n] = v,[ff-2]-v2[/i-l] + ;c[/i-l]-2e[n-l] + e[«]i w = 2,3,-" (4.16) 

Using (4.11) and (4.12), we write 

v,[n-2] = v,[«-3] + jc[«-2]-y[«-2]; h = 2,3,-- (4.17) 

V2[«-l] = V|[n-3] + V2[«-2] + x[«-2]-2y[n-2]; w = 2,3,--- (4.18) 

By subtracting these two equations we get 

v,[n-2]-V2[n-l] = y[n-2]-V2[n-2] = e[«-2]; n = 2,3,”* (4.19) 

Substituting (4.19) into (4.16) 

y[«] = x[«-l]-2e[n-l] + e[«] + e[n-2]; n = 2,3,-*- (4.20) 



Using (4.15) and (4.20), the output of the second-order modulator can be 
given in terms of the input, quantizer error and initial conditions: 

V2.o+<0]; « = o 

>'["] = ■ - V2_o + 40] + e[l] - 2e[0] ; n-l 

x[m - 1] + e[«] - 2e[n - 1] + e[rt - 2] ; n = 2,3, * • ■ 



(4.21) 
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4.3.2 Higher-order AS Modulator 

By proceeding in the same manner, for a third-order modulator the 
relationship between the output, input, quantizer error and initial conditions 
can be found as follows: 



>-[w] = 



V 30 +e[0\; 

‘'i.-i - 2 V 30 +401 + e[l]-3e[0]; 

- Vj,! +^30 +Jt[lJ + £^[2]-3e[l] + 3e[0]; 
jr [/7 - ]] + e[n]-3e[n - 1] + 2e[n~2]-e[n-3]', 



n = 0 



n = \ 
n = 2 

« = 3,4,-- 



(4.22) 



where v, , Vj.-i are the initial conditions on the first, second and 

third accumulators, respectively. The following theorem will now state a 
general description for an L'‘ -order AS modulator. 



«,lnl «jl«| Vj|n) 







► >i«| 



Figure 4-6. 11^ -order single-stage AI modulator 

Theorem 4.1: Consider the U' -order AS modulator shown in Figure 4.6 
employing M-bit quantization having maximum output levels of ±b . 
Assume that the input to the system x[n] is applied at n = 0. The initial 
conditions on the non-delaying accumulators at n = -i are defined as 
(where i =l,2,...,L-l),and the initial condition on the delaying 
accumulator (this is the last accumulator just before the quantizer) at n = 0 is 
^L.o = • Th^ output of the modulator y[n] in terms of the quantizer error 

sequence, input sequence and initial conditions is then given as below: 






V /..0 + 4O ]; 
/.-( 



« = 0 



+ x[n - 1 ] + ^ (- 1 / e[n - 



-1 



L 






(4.23) 



w = + 1, - •• 




ANALYSIS OF SINGLE-STAGE DELTA-SIGMA MODULATORS 



87 



where = y[«]- v^[«]; « = - , is the quantizer error sequence, and 



w ]_ ml 
k ]~\m-kYk\ 



(4.24) 



are the binomial coefficients, which were defined in Chapter 3 (see 
(3.116)). 

Proof: The proof is by induction on L . Observe that for L = 2 






v2,o +40]; 

v,_,-V2o+40] + 41]-2e[0]; 

x[n - 1] + e[n] - le[n -\] + e[n-2]; 



rt = 0 

rt = 2,3,--- 



and for L = 3 



y[n] = 



Vj,o + e[0]; 

v,,_i + - 2v, 0 + x[0] + e[l] - 3e[0] ; 

- Vj _| + Vj 0 + x[l] + e[2] - 3e[l] + 340] ; 
x[n - 1] + e[w] - 2e[n - 1] + 3e[w - 2] - e[n - 3] ; 



« = 0 
n = \ 
n = 2 
n = 3,4, 



(4.25) 



(4.26) 



Equations (4.25) and (4.26) are the same as (4.21) and (4.22), 
respectively. So (4.23) is true for L = 2 and L = 3. Now assume that (4.23) 
holds for any L>2. Using this we wish to prove that (4.23) is also valid 
L+1. 



L"" order AZ modulator 




Figure 4-7. {L + 1)'* -order single-stage AI modulator 

An (I + 1)"" -order modulator, which is formed by adding one more 
accumulation stage before an i"" -order modulator, is shown in Figure 4.7. 
The accumulator outputs within the Z,''' -order modulator are labeled as h^n] 
where i =l,2,...,L.The input to the Z,'^ -order modulator is x'[n], which is the 
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output of the accumulator preceding the L'^ -order modulator. At n = 0 we 
have 



y[0] = /.,,o+e[0] (4.27) 

By inspection, at n = 1 one can write (see Figure 4.7) 

L-i 

>'[l] = 40]->'[0](i + l)+r , +5^A„_, +^,„ +<1] (4.28) 

/=l 

where =x [-l] is the initial condition on the first accumulator in the 
(Z^ + 1)'* -order modulator. Using 40] = y[0]-^^o, the last equation can be 
given by 

L-\ 

= -ii,„ + x[0]+e[l]-(i + l>[0] (4.29) 

;=l 

The input to the I'* -order modulator x’ [n] is given in the form 

= + « = 0,l,2,--- (4.30) 

from which it can be written 

x’[n - 1] = x'[« - 2] + x[n - 1] - y[n - 1]; n = 1,2, ■ ■ • (4.3 1 ) 

By the induction hypothesis stated in (4.23), we have 



y[n] = 






+ r[rt-l] + ^(-l)*e[«-A:]|^^j; 



n = 0 



w = l,2,---,Z.-l 



n = L,L + i,--- 



(4.32) 



Specifically, for n =1,2,--, L-1 
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;<[«]= |;(- ir p )+(- 1)" ) 

From (4.33) it is immediate that 

" n-\ n 

'=» \ / \ / 

*=o V*J 



(4.33) 



(4.34) 



In order to obtain an expression for x’[n-2], we write (4.34) again as 
follows: 



r [h - 2] = )^[n-l]-^ 

/*«-! V / v” 



(4.35) 



Note that the last equation is valid for n = 2,3,---, L . Substituting (4.34) 
and (4.35) into (4.31) 

(4.36) 

+ ^(-l)'c(»-i][^l-||)(-l)‘e[«-l-i:]f^\ « = 2,3,...,Z, 



In order to proceed further we will make use of the following property 
[69], 

Property 4.1: For any positive integer numbers m and k , 



m\ ( m ^_r^ + l 

k k -\ r k 



(4.37) 



Proof: From the definition of the binomial coefficients stated in (4.24) 
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m\ 

{m-k)\k\ 

m{m - l)(w - 2) - ■ ■ (w - A + 1 ) 
_ 



Similarly, 



m\ 

^k -\ J (m-A + l)!(/t-l)! 

_ m{m - l)(w ~2)"-{m-k + 2) 

- (T-1)! 



Adding (4.38) and (4.39) gives 



m 

k-l 



m(m~])---(m~k + \) ^ m(m-l)---(m-k + 2)k 



k\ 



k\ 



which can be also be written as 

m \ m{m-\)---{m-k-\-2){m + \) 



k-\ 



k\ 



{m + 1 )! 
(m+I-A)!A! 



flj + l' 
V ^ . 



The terms in the right-hand-side of (4.36) can now be 
Property 4.1. However, first write one term in (4.36) as 



/-I ^ 
2 






/-I 

n-2 






By applying (4.37) in (4.41), we get 



/-I ) 
2 



( 4 . 38 ) 



(4.39) 



(4.40) 
simplified using 



( 4 . 41 ) 



( 4 . 42 ) 
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which can also be re-written as 



i-l / . • \ Z.-1 / • I \ 



V / \ / /««-! 

Another term in the right-hand-side of (4.36) is 




V(-l)”6.,« 



= (-!)■ Vo 
= (-l)”Vo 



'L-l^ 

'L-\ 

n 



(L-\ 

n-\ 



Using (4.37) in (4.44) results in 




A final term in the in the right-hand- side of (4.36) is 



Note that we have made a variable transformation in the 
summation asp = k + 1 . Further proceeding to write that 

g(- 1)* e[« - j- j 



( 4 . 43 ) 



( 4 . 44 ) 



( 4 . 45 ) 



( 4 . 46 ) 

second 



( 4 . 47 ) 



Using (4.37), we obtain 
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L + \ 



( 4 . 48 ) 



This can be summarized as: 






Substituting (4.43), (4.45) and (4.49) into (4.36) 

rz,' 



«-l n 



+ x[« - 1] + ^ (- l)* e[n - A] 



( 4 . 50 ) 



; n = 2A,---,L 



We now write y[n] for n = L,L+l,---, using again our induction 
hypothesis stated in (4.32), as follows 



L ^ / j ' 

y[n] = x'[n - 1] + ^ (- 1/ e[n - A] 

*=o , 



I; n = L,L + \,-- 



( 4 . 51 ) 



It follows that 



r[K - 1] = y[n]-^(-l)* e[« j; n = L,L + \,--- 



( 4 . 52 ) 



We again want to obtain an expression for x'[n-2 ] , and thus, write (4.52) 



as 



L 

r[«-2] = y[n-l]-21(-l)^4n-l-*] 



rr\ 






|: w — Z + 1, £ + 2, 



( 4 . 53 ) 



Substituting (4.52) and (4.53) into (4.31) gives 
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y[n] = i[n - 1] + ^(- 1)* e[« - j 



(4.54) 



I; w = Z + 1 ,£, + 2 , • 



The summation terms in the right-hand-side of (4.54) could be simplified 



to 






k 

\ f 



L ( I \ f I 



(4.55) 



Note that we again made a variable transformation in the second 
summation as p = k + i . Writing (4.55) again 



l/elw - ^][ ^ |- 1/ e[n-l- k] 






f. I' I \ ^ I 

*=o p=i \P * 



(4.56) 



= e[n] + (- ly-^' e[n~L-\] + '^{~\y e[n- ^:]|f ^ 
t=i IV , 






Using (4.37) in (4.56) 



^(- l)^e[« - A:]f^ 1- ^(- 1)* e[n - 1 - A] 



(L 



= e[n] + (- if"' e[« - Z, - 1] + ^(- 1)* e[n - k{^'l ^ 
*=i V > 



(4.57) 



which can also be written as 



^(- 1/ j- ^(- ^ X^' 



^L + \ 
k 



(4.58) 



By substituting (4.58) into (4.54), we get 
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t+i 

t=0 




w = I + l,L + 2,--- 



(4.59) 



We now define new labels for the accumulator outputs within the 
(I + 1)'* -order modulator as r,=v,_, and = v,^i_, where i = 1,2, •••,!- 1 , and 
finally ^ . We re-write (4.27), (4.29) and (4.50) using these newly 

labeled initial conditions 



y[0] = + e[0] 



(4.60) 



y[i] = Vj + 40] + 40- (i- + ikO] 

1=1 
i. 

= V|.^, + ^ ^/.-I - + 40 ] + 40 - (i + 1 )e[0] 

(=2 
L 



(4.61) 






(L 



‘'^A+I.O 



+ 4«-0+^(-iy4«-^] 



z, + i 
k 



I; « = 2 , 3 ,--,Z 



The last equation (4.62) could be simplified as 



4«]=T(-ir'vir 1 

+ 4« - 1] + ^ (- l)* e[n - ’’’ * 



; n = 2,X-’\l 



(4.62) 



(4.63) 



Using (4.60), (4.61), (4.63) and (4.59), the output of the (i + i)"' -order 
modulator is given by 
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K^i.o + ^[0]; 



X + ^[01 + e[\] ~{L + 1>[0] ; « = 1 

J|Hr.„(;:;)*Hrv„(;) 

+ x[«-i] + X(-iy4«-*l ; 

*=0 V * > 



n = 2,X---,L 



(4.64) 



« = I + l,L + 2,- 



The two terms in the above equation, namely for n = 1 and n = 2, 3, •••,/. 
can be combined into one term to obtain a final relationship 



ViHo+^0]; 



^ ^ ^ « = U-.L 



wit 7 4- 1 \ 

;c[«-l] + ^(-l)‘e[w-^] I; n 

*-o V ^ J 



= L + l,I + 2,- 



(4.65) 



which does comply with (4.23) and hence the proof is completed. 

4.4 No-Overload Stability Criterion 

In this section, we determine the criterion that ensures the no-overload 
operation of the quantizer within an I'* -order modulator. We start our 
analysis again with the second-order case, as it is the simplest form of the 
higher-order single-stage modulators, and then extend the results to the 
general case. 

In the previous section, we found that the output of the second-order 
modulator is given by (4.21). By definition v^in]^ y[n]-e[n] for all n. 
Hence, using (4.21) and this definition, the input to the quantizer (or the 
output of the second accumulator) is given by 

v2,o; « = 0 

^2 [n] = ■ Vi,_, - Vj 0 + x[0] - 2e[0 ] ; n = \ 

x[n-\]-24n-\] + e[n-2]', /i = 2,3,--- 



(4.66) 
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First assume that the initial conditions are appropriately set such that the 
quantizer operates in the no-overload region for n = 0 and n = 1. Remember 
that for an M-bit quantizer to work in the no-overload region, its input (this 
is the output of the second accumulator vjwl) should be confined in the 

range vjwje A1 for all n . Suppose that the maximum amplitude 

of the input is c (i.e., |a-[«]| < c for all n ). Further suppose that the quantizer 
has not been overloaded up to the time slot , where ^ 2 . This implies 
that |e[rt]| < A/2 for all « < . Using (4.66), we write the quantizer input at 

any Wq > 2 as 

v-2 [«o ] = -^[«o ' >1 - -11 + - 2] (4. 67) 



Using the well-known Cauchy-Schwarz inequality [22], 

|>^2K ]| ^ - ']| + - I]| +|eK - 2]| (4.68) 

and from (4.68) one can deduce that 

|v,K|<^' + 3A/2 (4.69) 

Since we want to keep the amplitude of the quantizer input within the no- 
overload operation region A, it can be concluded 

c = 2'"-'A-3A/2 (4.70) 

Notice that the above equation is only meaningful when M>2. This 
shows that for a second-order modulator, it is impossible to have a non- 
overload operation when a single-bit quantizer is utilized. When a single-bit 
quantizer is used in a second-order modulator, then the quantizer will 
certainly be overloaded. 

For the sake of minimum hardware complexity, we choose M = 2 (a two- 
bit quantizer) and in this case c = A/2. Now consider the second-order 
modulator shown in Figure 4.5 employing a two-bit quantization with output 
levels ±b . Assume that the input to the modulator is confined to be within 
half of the bin width A (i.e., |4/7]|<A/2 for all n ), where A = 2/)/3 . Further 
assume that the initial conditions on the first and second accumulator v, 
and ^2 0 , respectively, satisfy jvjol^^A and |v, - V 2 „| < A/2 , then the 

modulator will work in the no-overload region and the quantizer error 
sequence will be bounded to be within half of the bin width all times (i.e.. 
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|e[»]|<A/2 for all n). The following theorem will now give a general 

description for the no-overload stability criterion in an Z,'* -order modulator. 

Theorem 4.2: Consider the I'* -order modulator in Figure 4.6 that 
employs an L-bit quantizer with output levels ±b . Assume that |;c[«]| ^ A/2 ; 

for all n , where A = 2/>/(2^ -l) is the quantizer bin width. Further assume that 
the initial conditions satisfy the following criterion 

kJ<2^-'A (4.71a) 






i-U 











(4.71b) 



then the modulator is no-overload stable, and hence k[n]|<2^"‘A and 
|e[«]| < A/2 for all n . 

Proof: For the Z.'* -order modulator shown in Figure 4.6, we can write the 
quantizer input v^[«] by using Theorem 4.1 as stated in (4.23), and the 
definition for the quantizer error sequence, as before, 
n[«] = >'[«]-e[n]; n = 0,1,2,- 






L-\ 






+ - 1] + ^ (- ly e{n - 



^ «- 



w = l,2,---,I-l 



n = Z.,L + l,- 



(4.72) 



At n = 0 we have |v^(,|<2^'A, implying |e[0]|^A/2 (see also Figure 4.3 

and 4.4 and the discussion in Section 4.2). Assume that the quantizer has not 
been overloaded up to time slot , where \<nQ<L . This suggests that 
|<;[«]| < A/2 for all « < . At time , where I ^ < Z, , using (4.72) we write 
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L-\ 



Wn 

\ " / 



'■="0 V ® > 

+ x[n^ - 1 ] + ^ (- 1 )* ^[Wo - *]f ^ 

A =1 V , 

Again using the Cauchy-Schwarz inequality [22] results in 

'=% V J 

+ |4«o-l]| + ^|e[Ho-^]|[ 



(4.73) 



L-l 



(4.74) 






Since |4«]|<A/2 for all n, and |e[«]|<A/2 for all «<Wo, using the 

criterion for the initial conditions given in Theorem 4.2 (4.71a) and (4.71b), 
one can write 



vJ«,]|<2'-'A 




which is simply 
|v,K]|^2^-'A 



(4.75) 



(4.76) 



where . 

Now assume that the quantizer has not been overloaded up to time slot 
, where riQ>L. This suggests that \e[n]\ < A/2 for all w < Wo . At time Wq 

where %>L, using (4.72) we write 



L 

[«o ] = -^[«o - 1] + 1/ e[«o - k] 



t: 



(4.77) 



From (4.77), we obtain 



L 

|vi [«o ]| ^ - i]| + 

k=\ 




Since |x[/j]| < A/2 for all n , and |^[«]| < A/2 for all w , 



(4.78) 
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(4.79) 



We know that for any real numbers a and b , and for any non-negative 
integer number L we have 




(4.80) 



This is the Binomial expansion stated in (3.74). Setting a = ^ = 1, we can 
obtain 




(4.81) 



(4.82) 



Using (4.82) and (4.79) results in 

lvj«.]|<2‘-'i (4.83) 



for any Wq > /, . This completes the proof. 



4.5 Solution to the Non-linear Difference Equation 

The following theorem will now give the solution to the non-linear 
difference equation in (4.66) for a second-order system. 

Theorem 4.3: Consider again the second-order modulator shown in 
Figure 4.5 employing two-bit quantization with output levels ±b . Assume 
that the input to the modulator is confined to be within half of the bin width 
A (i.e., |jd[«j| < A/2 for all n ), where A = 26/3 . Further assume that the initial 

conditions satisfy |v2p|<2A and |v,^, - ^ A/2 , so that the modulator is 

overload free (i.e., |v2[w]|<2A and |e[n]<A/2| for all «). The normalized 
quantizer error sequence is then given by 
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Vw 



( 4 . 84 ) 



Proof: Since the modulator is overload free, we make use of (4.4) to 
write 




2 




« = 0 , 1 , 2 ,- 



Using (4.85) and (4.66), at n = 0 we write 




( 4 . 85 ) 



( 4 . 86 ) 



and at n = 1 

1 /v,,_|-V2o+x[0]-2£r[0] 

z 



Substituting (4.86) into (4.87), one obtains 




^ Vi+V2.O+40] ^ 



( 4 . 87 ) 



( 4 . 88 ) 



where we made use of the fact that 

± = {a±Fj ( 4 . 89 ) 



for any integer m , n , and real numbers a and b . 
Using (4.84), we write 



£[0] = :-- and e[\] = -- 



^^.-1 + ^ 2.0 + 40 ] 



This means that (4.84) satisfies the real solution for n = 0 and n= 1. Now 
assume that (4.84) is true for any n'^2 . Using (4.85) and (4.66), the 
normalized quantizer error sequence e[n] at n + 1 can be written as 






( 4 . 90 ) 
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By induction hypothesis we have 



g[«] ^ 1 / ^2.0 , 



A 2 






(4.91) 






Substituting (4.91) and (4.92) into (4.90) gives 






2 \.yfM4K 



*-o V / *-o V. / / 



which after some careful manipulations could be simplified as 



e[« + l] = --P ^ ^ 

2 \ rrr«-n v\k^ lY 



,[ 4^-1] , 1 

A 2 



4 M 4 2/ ' 



(4.93) 



(4.94) 



Note that we have again used (4.89). Further manipulations yield 

(4.95) 



which proves the theorem. 

Theorem 4.4: Consider the L"" -order modulator shown in Figure 4.6 
employing L-bit quantization with output levels ±b . Assume that the input to 
the modulator is confined to be within half of the bin width A (i.e., 
|;c[n]|SA/2 for all n), where As2ij/(2^-l)* Further assume that the initial 
conditions satisfy the criterion of Theorem 4.2, so that the modulator is 
overload free (i.e., |v^[n]| <2^''A and |e[«]<A/2| for all n ). The normalized 
quantizer error sequence is then given by 
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A 2 




;/7 = 0 , 1 . 2 , ••• ( 4 . 96 ) 



where the jV'^' -order summation operator ) is defined by 



‘^,v(4"])=^ ^ (4.97) 

i V=<) iS|=0 

Proof: As a matter of convenience, we start by introducing a notation for 
the A^'^' -order difference operator V^ (), such that 



;V 

V,v («[«]) ^ X 
*=0 



iV 



(4.98) 



Note that by definition, zero-order differentiation of any time sequence 
gives the sequence itself or V„(4/7j) = t7[rt]. 

As the modulator is guaranteed to work in the no-overload region, we 
can use (4.4), and write 




2 







(4.99) 



Observe that (4.96) satisfies the initial conditions. Now assume that 
(4.96) is true for any n>L. Using (4.99), (4.72) and the definition in (4.98) 
one can write the normalized quantizer error e[n] atn+1 as 



£[« + !] = -!- 
2 



/ 4^] I ^ i e[/^ + lJ 
\ A 'I A 



e[n + IJ 



(4.100) 



By our induction hypothesis stated in (4.96) 



'’/..O . (,«+! 



e[n+l] _ i / A 
A 2 ' 



+ 5," 



A 2 



'+ sr 






A 2 



1 -I !■ 6'!^ 






A 2 






I A 2 



(4.101) 
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Substituting (4.101) into (4.100), and using the linearity of the difference 
operator together with (4.89), we get 




We now give some useful properties regarding the difference operator. 

Property 4.2 [22]: L‘' -order (where L > 1 ) differentiation of a constant is 
zero. 

= 0 (4.103) 

Proof: Using (4.98), we have 

We know that for any integer number L > 1 , 

(4.104) 

Thus the proof is completed. 

Property 4.3 [22]: U*' -order differentiation of an -order summation is 
(n -order summation. Mathematically, 

V,(C(«[«]))=s;-i («[«]) (4.105) 

Proof: The proof is by induction on L . Observe that for L = 0 (4.105) 
holds. Now assume that (4.105) is true for any L . From the definition of the 
difference operator we have 



(«[«])- ^t-i («[«])- ^L-\ («[« - 1 ]) 



(4.106) 
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Now we write the Z,'* -order differentiation of an iV'* order summation 
using (4.106) 

V,., fe"' («[«]))= V, fc*' («[«]))- V, fc (a[n])) (4. 107) 

Using the induction hypothesis, (4.107) could be written as 

V,., (sr' («[«]))= 5") (a[«])- 5;^, („[«]) (4. 1 08) 



Using the definition of the summation operator in (4.97), 



*1=0 



( 4 . 109 ) 



^S-L 



M- I t -t.K, 

*/V-i=0*A'-i.-l=0 *1=0 



( 4 . 110 ) 



Subtracting (4.109) and (4.110) yields, 

*^-i=''*jV-t-l=0 *1=0 *A'-/.=0*Ar./,-I=0 *1=0 

- £ £ -£4m 

*.V-t=0*Ar_/.-l=0 *1=0 

Further simplifying, we get 
5;!)(oW)-5;-.k«])= £ *£ - 

*A/-i=«*A'-i-l‘=0 *1=0 

= £ *£' 

*A'-i-l=0*A'-t-2=0 *1=0 



( 4 . 111 ) 



Substituting (4.111) into (4.108) proves the property. Two direct 
consequences of this property, which might be very useful in evaluating 
(4.102), follow: 

v,.(sr («[«]))= 0 (4.112) 



where L> N . And, 
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( 4 . 113 ) 



where L = N . Equation (4. 112) states that -order differentiation of an 
N'*' -order summation is zero, when the order of the differentiation is greater 
than that of the summation. On the other hand, when the order of the 
differentiation is equal to that of the summation, (4.113) dictates that Z."’- 
order differentiation of an -order summation is the function itself inside 
the summation. Note that both of these rules are well known in calculus [22], 
here they are proved and stated as a matter of completeness. 

By making use of{4.112), (4.113), and (4.103), (4.102) becomes 



I / A A 
£[« + !] = --( 



1 4^0 _ x[n] _ j ^ ^ 




A 2 



\ 



I A 2 



Using (4.89) and simplifying we get 



£[« + !] = --( 



A 

,+- + 5 ^^ 



[a 2j ' ( A 2J 



( 4 . 114 ) 



( 4 . 115 ) 



which does agree with (4.96), and completes the proof. 

As in the previous chapter, we again confine our attention to only dc 
inputs x[n] = ATe [- A/2, A/2] for all n , and thus the normalized quantizer error 
can be given by 



e[«] = --a«]; « = 0,1,2,- 



( 4 . 116 ) 



where 






a«]=( 






+ S" 



+ + - |+ 5 "(^) 




( 4 . 117 ) 



with /3 = A'/A + l/2 as before. 
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4.6 Statistics of the Quantizer Error Sequence 



In this section, the long-term statistics of the normalized quantizer error 
sequence are studied. Notice that the solution of the difference equation in 
this case remarkably resembles that of the MASH modulators, and hence the 
theoretical framework developed in the previous chapter will be applied 
here, without loss of generality. As a matter of clarity and for the sake of 
brevity, we follow a similar notation to that of the Chapter 3, and only give a 
brief outline of the derivations. 

The nested summations in (4.1 17) could be simplified as 



A 



A 2 






1 \ 



' V, 






( 4 . 118 ) 



We have used Lemma 3.3 in (3.82). Recall that 5?i(«) represents a 
generic polynomial of degree L in n . We again perform our analysis for the 
special case of the rational dc input (i.e., A is a fraction of the quantizer bin 
width A ). It is further assumed that the initial condition on the first 
accumulator is an irrational number (i.e., Vj _, cannot be represented as a 

fraction of the quantizer bin width A ). With these specifications, the first 
and second moments of the sequence are given by: 






( 4 . 119 ) 






( 4 . 120 ) 



Note that (4.119) and (4.120) are only valid when L>2. Using these 
equations, one can determine the mean value A/(e) and the average power 
M(e^) of the normalized quantizer error sequence as 

M(e)=0, 



( 4 . 121 ) 
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Additionally, the cross-correlation between the input and the normalized 
error is 

r,,(/) = M{Xe)=M{e)X = 0; V/ (4.122) 

which implies that provided an irrational initial condition is imposed on 
the first accumulator in a higher-order single-stage NL modulator, the 
quantizer error sequence has a zero mean and an average power of A^/l2. 
Furthermore, it is uncorrelated with the dc input. We will now derive the 
autocorrelation of the normalized quantizer error in second-order and higher- 
order modulators. 



4.6.1 Second-order A£ Modulator 

For a second-order modulator we have 



C[«] = 




(4.123) 



Note that this differs from the second-order 1-1 MASH case by only the 
sign of the «/2 term. So the derivation is almost the same, and hence only 
a brief outline will be given here. 

Using again the Fourier Series representation for the modulo 1 operator 
stated in (3.92) and (3.93), one can write (4.123) as 



(/) = ihk) 



(4.124) 



where g() is the Fourier Series coefficients given in (3.93), and 4>,(/,/r) is 
the joint characteristic function of ^[«] , namely, 

* (a)= (4,125) 



An application of the Wely’s Theorem, which was stated in (3.47), to 
(4.125) would result in 

0 , 



<&,(;.*)= 



< + A: = 0 and /%/ = int eger 
otherwise 



(4.126) 




108 



ANALYSIS OF SINGLE-STAGE DELTA-SIGMA MODULATORS 



Hence, the autocorrelation of the quantizer error in a second-order single- 
stage modulator is 



12 ’ 

1 _ 



/ = 0 



where this time y = + ^ + + ^) 



(4.127) 



and /j is defined in (3.107). 



4.6.2 Higher-order AS Modulator 

Because the solutions of the non-linear difference equations are exactly 
the same, one can follow the same mathematical methodology as in Chapter 
3 to derive the autocorrelation function of the quantizer error sequence. The 
derivations are quite straightforward. For a higher-order single-stage system 
having a rational dc input and irrational initial condition on the first 
accumulator, we have 

f[«] = (3,,(«) + S,(«)) (4.128) 

The autocorrelation of ^[«] is 

= (4.129) 



where 



1 _ 

0,(/,^)= lim — ^ 






(4.130) 



Applying Wely’s Theorem to (4.130), 
/ fl, 



(4.131) 



It should be also highlighted that Wely’s Theorem can only be used here 
when the order of the modulator is greater than two (i.e., L>2). 

Using (4.131) it immediately follows that 
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'■.(/)= H’ ' = “ (4.132) 

[o, 

4.7 Fundamental Result 

Our theoretical results have shown that, once the initial condition on the 
first accumulator is set to an irrational value (or it is not a fraction of the 
quantizer bin width), then the other initial conditions on the remaining 
accumulators will not have any effect on the long-term behavior of the 
quantizer error sequence. It is, therefore, very desirable to update the no- 
overload stability criterion in Theorem 4.2 by setting only the first initial 
condition to be a non-zero value. This would lead to a simpler expression 
and a better understanding of the effect of the initial condition. The 
following theorem will now give this updated no-overload stability criterion 
along with all our key findings. 

Fundamental Theorem: Consider the -order modulator in Figure 4.6, 
which employs an L-bit uniform quantizer with output levels ±b. Assume 
that \x[n]\<A/2'y for all n, where A~2b/i^^ ~\) is the quantizer bin width. 
Further assume that the first accumulator initial condition is set to an 
irrational number satisfying |v,_,| <2^ 'A-(L + l)A/2 , whereas all other initial 

conditions are set to zero. Then the modulator is no-overload stable, and 
hence |v^[n]| < 2^"'A and |^[«]|^A/2 for all n . Moreover, the quantizer error 

sequence is consistent with the signal-independent uniform distribution 
white noise approximation, when L > 3 . 

4.8 Simulation Results 

In order to validate our theoretical developments, in this section, we 
present autocorrelation simulation results for various modulators. A time- 
difference equation model was developed in the C programming language 
and the outputs of the simulations were passed to Matlab™ for further 
processing and display. 

Figure 4.8 depicts the simulated autocorrelation function in the second- 
order modulator along with the analytically derived results, for two different 
input and initial condition combinations. The exact matching between our 
theoretical results and simulations strongly and clearly support the validity 
of equations (4.121) and (4.127). 




no 



ANALYSIS OF SINGLE-STAGE DELTA-SIGMA MODULATORS 



mean value=1 .411 54373e-S, average power^O. 08333440 
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(b) 

Figure 4-8. Autocorrelation of the quantizer error sequence in the second-order modulator: (a) 
X/A = 0.05 and v, _, = 0.01929710159476, and (b) A'/A = -0.1 and 

v,,_, =0.11762271073900 
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Figure 4-9. Autocorrelation of the quantizer error sequence in the third-order modulator: (a) 
X/A = 0.1 and v,__, = 0.46466656988786, and (b) = and 

v,,,i =0.00563502894910 

The simulated mean value and average power of the quantizer error 
sequence is also listed in the plots. As predicted by our theory, the mean 
value of the quantizer error sequence approaches zero and the average power 
nearly equals to aV 12 . For the second-order case, the autocorrelation does 
not have zero values for non-zero lag indexes. Therefore, it cannot be 
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regarded as white noise even when an irrational initial condition is used in 
the first accumulator. This completely agrees with our findings. 

Similar results for the third-order system is shown in Figure 4.9. It is 
clear from these figures that the autocorrelation of the quantizer error 
approximates an impulse, which complies with equation (4.132). Moreover, 
the mean and average power are zero and NjM , respectively. 

The quantizer error from a third or higher-order modulators (of course 
when an irrational initial condition is used) does exhibit identical first and 
second order moments, which are the same as those for a uniform 
distribution white noise. 

4.9 Output Spectrum 

Output noise spectrum (the dc input is removed) of a second-order 
modulator is shown in Figure 4.10(a) and (b), for the zero and irrational 
initial conditions, respectively. 

As for all of the spectral density estimation results presented in this 
research monograph, these spectra were obtained using the Welch’s method 
[70] of averaging modified peridograms, as explained in Chapter 3 Section 
3.5.As illustrated in Figure 4.10(a), when the initial condition is zero, the 
output quantization noise exhibits pure oscillations and limit cycles, and 
there is very poor noise shaping. 

An irrational initial condition on the first accumulator, however, 
appreciably reduces the tones, and thus provides better noise shaping, as 
shown in Figure 4.10(b). Nevertheless, this improvement does not totally 
eliminate the tones throughout the whole spectrum. 

An irrational initial condition on the first accumulator in the third-order 
modulator, however, is very effective in completely eliminating the tones as 
shown in Figure 4.11. Third, fourth and fifth-order modulators with an 
irrational initial condition on the first accumulator are simulated for various 
dc input amplitudes and no tones are observed throughout the whole 
spectrum (see Figure 4.12). 
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Figure 4-10. Output quantization noise spectrum for second-order modulator where 
X/A = 0.05 , = 0 : (a) v, = 0 , (b) v, _, = 0.077046171 19143 
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Figure 4~IJ. Output quantization noise spectrum for third-order modulator where 
X/A = 0.05, V 2 ,_, =0, V3o-0:{a) v,.,^0,(b) v,_, =0.34676719059531 





ANALYSIS OF SINGLE-STAGE DELTA-SIGMA MODULATORS 



115 




02 

digital frequency 



3 

|--200 

O 

05 



(b) 



116 



ANALYSIS OF SINGLE-STAGE DELTA-SIGMA MODULATORS 




(c) 



Figure 4-12. Output quantization noise spectrum for modulatore: (a) third-order, (b) fourth- 

order, (c) fifth-order 



4.10 Error-Feedback Topology 

Because the governing difference equations are the same, our theoretical 
derivations are also valid for some other AS topologies. For instance, Figure 
4.13 shows a third-order error-feedback topology. This modulator is 
simulated for a normalized dc input amplitude A'/A = 0.05 with zero and non- 
zero irrational initial conditions. 






— [ZD^ 






thrce-bii 

quantizer 









initial condition 



Figure 4-13. Third-order error-feedback topology AI modulator 




ANALYSIS OF SINGLE-STAGE DELTA-SIGMA MODULATORS 



117 



Figure 4.14 demonstrates the output quantization noise spectra for this 
modulator and the effectiveness of the initial condition is obvious. 




(a) 




(b) 

Figure 4-14. OuQ)ut quantization noise spectrum for third-order error-feedback 
modulator where X/A = 0.05 : (a) ic = 0 , (b) ic = 0.52394539423623 
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4.11 Conclusion 

This chapter has proved that higher-order AZ modulators are capable of 
producing tone free output spectra (for dc input up to maximum depth of 
modulation) when an irrational initial condition is used. It has been long 
believed that higher-order modulators are less prone to generating limit cycle 
oscillations and tones. Here we have presented a rigorous mathematical 
analysis that supports and proves this hypothesis. To the author’s best 
knowledge, this is the first time this has been quantitatively and rigorously 
proved. 

We started our analysis with the second-order single-stage AZ modulator, 
which consists of two cascaded integrators (the first one is non-delaying and 
the other is delaying). The analysis was then generalized to any order case. 
A non-linear difference equation was derived that relates the output of an 
L’’ -order modulator to the input and initial conditions on the integrators 
(accumulators). From this non-linear difference equation, we found the 
necessary conditions that must be met to operate the modulator in the no- 
overload region of the quantizer. By solving the non-linear difference 
equation, we have obtained first and second-order moments of the quantizer 
error sequence in an D' -order modulator. The long-term statistics of the 
quantizer error sequence in a higher-order single-stage AZ modulator were 
found to be the same as for the higher-order MASH type modulator. 

The no-overload stability criteria states that in a single-stage modulator 
the number of bits in the quantizer should be the same as the order of the 
modulator. Furthermore, the maximum input amplitude cannot exceed half 
of the bin width. In practice, however, second-order single-stage AZ 
modulators with a pure differentiation NTF as shown in Figure 4.5, which 
employs a single-bit quantizer, have been designed and used in commercial 
products [7]. Although, the quantizer significantly overloads, the modulator 
still works reasonably, because in this case it is Bounded-Input Bounded- 
Output (BIBO) stable. The BIBO stability is looser than the no-overload 
stability, thus it may allow some useful modulator designs (with a single-bit 
quantizer) to pass the stability test. However, the second-order single-stage 
modulator consisting of two cascaded integrators is a special case. The same 
statement is not at all true for third or higher-order modulators. It was shown 
in [34] that higher-order single-stage modulators, which employ single-bit 
quantization, with pure differentiation NTF are not BIBO nor, of course, no- 
overload stable. The exact mathematical analysis of higher-order single- 
stage AZ modulators with a single-bit quantizer is not possible, at least with 
the method used here. To the author’s best knowledge, there is no reported 
exact analysis result in the state-of-the-art literature that concerns the 
overloaded quantizer error. 
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5. FRACTIONAL-N PLL FREQUENCY 
SYNTHESIZERS 

This chapter focuses on the design, modeling and simulation of 
Fractional-N PLL synthesizers as well as the pipelined implementation of 
MASH type AX modulators to be used in such synthesizers. We start by 
providing the motivation for the development of a simulation model for PLL 
frequency synthesizers. Next, a linearized frequency-domain model of the 
PLL frequency synthesizer is derived, serving as a very useful tool to 
understand the basic operational principles. Following this, the main 
limitations of the traditional synthesizers are discussed and the Fractional-N 
concept is reviewed to reveal how fine frequency resolution can be obtained 
without compromising and corrupting the phase noise performance. Later, a 
simulation model is developed for a complete Fractional-N PLL synthesizer. 
This chapter concludes by presenting design considerations and 
implementation aspects of pipelined MASH AX modulators suitable for 
Fractional-N PLL frequency synthesis applications. Specifically, in an effort 
to reduce the hardware complexity and power consumption the alignment 
registers, which are normally required in pipelined adders, are eliminated by 
taking advantage of the static modulator input. 
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5,1 Introduction 

PLL based frequency synthesis techniques are the most popular 
approaches to providing accurately defined frequencies and desired channel 
selection in wireless transceivers [60], [61]. The spectral purity and timing 
specifications of various wireless standards, for instance Global System for 
Mobile communication (GSM) [73] or Digital Enhanced Cordless 
Telecommunication (DECT) [74], [75], put stringent requirements on the 
synthesizer performance criteria such as phase noise and settling time. To 
mitigate the problems encountered with the classical PLL synthesizers, a 
technique called “Fractional-N synthesis” has recently received a lot of 
attention for providing step sizes much smaller than the reference frequency 
[64]. Unlike traditional PLL frequency synthesizers, in the Fractional-N 
approach the divider modulus in the feedback loop is allowed to vary 
between predetermined integer values so that in the time average sense an 
accurate fractional division ratio is obtained. This technique allows very 
narrow channel spacing when compared to the reference frequency and 
relatively wide loop bandwidths while improving the phase noise 
performance of the PLL [61]. In the Fractional-N method, the fractional 
division is achieved by varying the modulus of a frequency divider upon the 
output of a digital A£ modulator. 

Both the PLL and AS modulator are non-linear systems and the 
mathematical descriptions, which rely on the small signal model of the PLL 
components and white noise quantizer error assumption for the modulator, 
are not enough to predict the actual performances (phase noise, spurious 
level and settling time). PLLs are inherently analog systems and are usually 
implemented by using discrete components or on-chip mixed-signal VLSI 
techniques [60]-[64]. There is, however, a lack of information in the 
literature concerning the simulation of Fractional-N PLL synthesizers. The 
difficulty involved in the simulation of PLLs comes from the fact that the 
time constants of different components within the system vary over a very 
wide range. To further elaborate on this issue, consider the AZ modulator 
alone. Like pseudo-random generators, many samples (on the order of 
millions) are required from a AZ modulator output to accurately predict its 
noise shaping performance [7]. Now assume that the division in the feedback 
loop within the PLL is set to 40. For this condition, the clock frequency of 
the modulator will be 1/40 of the PLL output frequency. Assuming that the 
system is simulated for 1 Million AZ output samples and taking into account 
the Nyquist frequency for the PLL output, a simple calculation reveals that 
the number of samples that must be acquired at the PLL output should be at 
least 80 Million. If the circuits that constitute the Fractional-N PLL was to 
be simulated at the transistor-level with SPICE [10], it would take an 
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exorbitant amount of time. This is certainly not an efficient way forward as 
an initial design tool. 

It is therefore obvious that a behavioral level simulation model is needed 
to accurately define the PLL parameters before the real implementation, and 
hence save considerable turnaround time in the design process. A second 
benefit that can be obtained from such a simulation model is that the 
spurious suppression performance of a AZ modulator with different initial 
conditions can be easily assessed. For these purposes, we have developed a 
difference equation simulation model on a mixed MATLAB”^^ [76] and C 
[77] platform, where all the components of the PLL have been considered to 
be ideal, except the AS modulator for which a fixed-point exact model was 
used. 

5,2 Analysis of PLLs 

A PLL is basically a feedback network, where the phase of the output 
signal is locked onto that of the input signal [78]. A typical PLL frequency 
synthesizer is shown in Figure 5.1, consisting of a Phase Detector (PD), a 
low-pass loop filter, a Voltage-Controlled Oscillator (VCO) and a frequency 
divider. The input and output signals have frequencies of and (in 

units of Hz) and phases of and (in units of radians), respectively. 

The main function of the PLL is to track the small differences in phases 
between its input and feedback signal derived from the divider. This is 
accomplished by measuring the phase difference in a PD, which generates a 
voltage v^(f) that is proportional to the phase difference. This voltage is then 
filtered out by a low-pass filter and applied, as a correction voltage Vj.(r), to 
the VCO to adjust the output oscillation frequency in such a way that the 
phase error is reduced. 




Figure 5-1. Block diagram of a traditional PLL frequency synthesizer 

Through negative feedback in conjunction with frequency division, the 
PLL results in = Nf^j - , where N is the division ratio of the divider. More 
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specifically, when the loop is locked the average frequency of the VCO is 
approximately equal to N times the input frequency (i.e. for each input cycle 
there are N output cycles). The frequency divider in the feedback loop is an 
edge triggered digital device and thus its division ratio is restricted to be an 
integer. Therefore the output frequency of the PLL can only be represented 
by integer multiples of the input reference frequency. That is, can be 
varied in increments of as N is changed in an appropriate manner. 

5.2.1 Small-signal Model 

Although the real mathematical description of the PLL is quite 
cumbersome, it can be accurately modeled as a linear system when the loop 
is in the lock condition [61], [78]. Specifically, when the PLL is locked the 
PD acts as a linear element and its output v^{t) can be modeled as follows: 

vAO = RA<t>,.rU)-<!>^,iO/N) (5.1) 

where is called the PD constant and measured in units of V/rad, and 
the reference phase is radians. The PD compares the phase 

of the reference signal with the phase of the divided-down output, and 
produces an output voltage proportional to this phase difference. The output 
of the PD is composed of a dc component and a superimposed ac 
component. To suppress the undesirable high-frequency ac components, 
vAt) is then processed by a low-pass loop filter, which has a transfer 
function, say, F{s). The VCO has an oscillation frequency , which is 
controlled by the output of the loop filter v^(/). The instantaneous angular 
frequency (in radians/sec) at the VCO output can be written as: 

^oAt) = MouA) = Wo + RMt) (5.2) 

where and are the free running frequency of the VCO and VCO 
constant measured in Hz and rad*HzA^, respectively. Since phase is the 
integral of the angular frequency [13], one can write the VCO output phase 
as: 






(5.3) 
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Assuming that all the frequency and phase variables in the system are 
measured with respect to the free running frequency of the VCO, , and 
taking the Laplace Transform of (5.3) and (5.1), we get [78] 

= (5.4) 

S 

and 

VM = (5.5) 

Having established these small-signal properties, we can now produce a 
linear model for the PLL. A signal-flow diagram based on the small-signal 
properties of the PLL components is shown in Figure 5.2. Such a model can 
be used as an initial design tool to roughly set the PLL parameters. Fine 
adjustments must be obtained by means of simulations. This model relates 
the changes in reference input phase to the output phase, and is valid as long 
as these changes are slow and small about an operating point [78]. It is also 
worth noting that there are two different frequency variables in this 
development. For instance the input as well as the output of the PLL are 
periodic voltages with temporal frequencies of and , respectively. On 

the other hand, in (5.5) we have seen that the output of the PLL is also a 
function of the complex frequency variable s (Laplace variable). This can 
lead to a confusion and must be handled carefully. In fact, the complex 
frequency variable s here represents the rate of change in output frequency 
(or phase) about the free running oscillation frequency (or phase) of the 
VCO. 

Phase Detector 



Figure 5-2. Linear model of a PLL frequency synthesizer 

Carefully examining this figure and working around the loop yields the 
following input-output transfer functions: 
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H(s) = 






0) 



E,K,F(s) 



s + KjK,F(s)/N 

The phase error transfer function is given by: 



(5.6) 



E(s)^ 









s^K^^f\s)iN 



(5.7) 



5.2.2 Second-order Systems 

The transfer function equations derived in the previous subsection are 
general in that they apply to almost every PLL. What differs in the dynamic 
responses between different PLL designs is the choice of the low-pass filter 
transfer function F(s) (in the ideal case). In general, a very simple first-order 
filter is commensurate as the low-pass loop filter in PLL designs. For 
example Figure 5.3 shows a simple lag-lead low-pass filter with transfer 
function: 



1 + JT, 



with time constants t, = (/?, + )C and = R 2 C . 




(5.8) 



Figure 5-3. A passive lag-lead low-pass filter 
Inserting the transfer function of (5.8) into (5.6) and (5.7), one can obtain 
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H{s) = 





24 - — 


■\-Nwl 












and 



E{s)^ 



N 2 

s + w 



+ 2(^W„5 + w„ 



where 






V 



(5.9) 



(5.10) 



(5.11) 



e * 



T, +- 









(5.12) 



Because the highest power of s in the denominator of these transfer 
functions is two, the PLL is classified as a type-1 second-order system. 
Second-order systems received a lot of attention in the literature, because of 
the simplicity of the low-pass loop filter and because the input-output 
transfer function can be represented in a standard form. By this way many 
useful results from control theory can be directly used in the analysis of 
PLLs [78]. In agreement with servo terminology from control theory, w„ 
(which is a measure of the PLL settling time) and | are usually called the 
natural frequency and the damping factor, respectively. 




Figure 5-4. A feedback control system 
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The type of the PLL refers to the number of poles of the open-loop 
transfer function located at the origin. To further elaborate on this issue, 
consider a basic linear feedback system as shown in Figure 5.4, where G{s) 
and B{s) are called the feed-forward and feedback transfer functions, 
respectively. Comparing the block diagrams in Figure 5.2 and Figure 5.4, the 
PLL can be considered as a feedback system with 



G{s) = 



K,K,F{s) 

s 



(5.13) 






(5.14) 



Using the transfer function for the passive lag-lead low-pass filter given 
in (5.8), the open-loop transfer function 0(5) ofthe PLL is given by 



0{s) = G{s)B{s)^ 



K,KX\^sxY) 

i'A^(l + 5T,) 



(5.15) 



It is clear from (5.15) that the open-loop transfer function of the PLL has 
one pole located at zero frequency, and hence the name type-1. The main 
drawback of type-1 systems is that the steady-state phase error does not 
approach zero for step changes in reference frequency (corresponding to a 
channel switch in frequency synthesis applications). To see this, let us 
consider a step change of A/^y in the reference frequency. The Laplace 

transform of the reference phase is now given by 0^y(5) = 2;iA/’„y/5'^ (recall 

that phase is the integral of the instantaneous frequency with respect to 
time). Also note that the Laplace transform of the phase error is defined as 
d>^(5) = 0^^(5)-<I>^„,(5)/A^. Using the definition for the phase error transfer 

function stated in (5.7), we get 



(D,(^)=£(5)0^y(5) 



(5.16) 



which can also be written as: 






s{s + KjK,FU)/N) 



(5.17) 
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The steady-state phase error, which is the value of the phase error after 
all the transients have died out, can be found from the final value theorem 
[78]: 



lim0g(/) = limsOg(5) 



( 5 . 18 ) 



Substituting (5.17) into (5.18) gives 



lim0,(O = 






( 5 , 19 ) 



Unless the dc gain of the loop filter is very large, the steady-state phase 
error will not be zero. 

To have zero phase error response to step changes in reference frequency 
the low-pass filter in the PLL should ideally have infinite dc gain. The 
realization of very very large (infinite) dc gain, however, requires an active 
filter as shown in Figure 5.5. Assuming that the Op-amp is ideal the transfer 
function of this active filter is: 



F{s) = - 



l+^T; 

JT, 



( 5 . 20 ) 



where t, =/?,C and Tj =RiC . The advantage of the active filter is that it 
can be designed to provide near ideal integration. Substituting (5.20) into 
(5.6) and (5.7), we get 



H{s) = 






( 5 . 21 ) 



and 









( 5 . 22 ) 



where this time 
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(5-24) 

With this active low-pass filter, the open-loop transfer function is given 
by: 



OU)=^~ 



KjKA^ + sxA 
s'^Nx^ 



(5.25) 



It is now clear from (5.21) and (5.25) that the active filter in Figure 5.5 
makes the system a type-2 second-order PLL. While type-2 PLLs benefit 
from the reduced steady-state phase error response, the presence of an Op- 
amp in the active filter introduces significant challenges in the 
implementation [60], [61], and [78]. Specifically, the Op-amp is very noisy 
and may severely corrupt the phase noise performance of the PLL, and 
therefore it is undesirable in high-performance systems. 




Figure 5-5. A first-order active low-pass filter 
5.2.3 Charge-pump PLL 

A simple way to realize a type-2 PLL without using an active filter is the 
so-called charge-pump technique. Shown in Figure 5.6 is a charge-pump 
PLL, which includes a Phase-Frequency Detector (PFD) and a 
corresponding charge-pump loop filter. Among other PDs such as four- 
quadrant multiplier, EX-OR gate or J-K flip-flop, the PFD has some 
remarkable advantages such as offering a virtually unlimited pull-in range‘s, 
which guarantees the PLL acquisition [78]. Additionally, the PFD also has 
the ability to measure the frequency difference between its inputs, and thus 



The pull-in range of a PLL is the frequency range within which the system will always 
become locked (although the process is much slower when compared to the normal 
operating range). 
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enables the PLL to attain lock very quickly [61], provided does not 
exceed the normal operating range. 




Figure 5-6. Charge-pump PLL frequency synthesizer 

As its name implies, the output of the PFD is not only dependent on the 
phase difference, but also on the frequency error while the PLL is acquiring 
lock. Figure 5.7 shows one possible implementation of the PFD, which 
basically consists of two edge-triggered D-type flip-flops and an AND gate. 
The PFD generates two outputs Up and Dn according to its input voltages 
Vin and Vdiv, which are digital signals, e.g. square waves. Specifically, Vin is 
a periodic square wave oscillating at a frequency of and Vdiv is the 

output voltage of the frequency divider in the feedback path. The output 
signals. Up and Dn, from the PFD are used to control the VCO output 
frequency by applying positive or negative charge to a charge-pump loop 
filter as shown in Figure 5.8. 

The PFD either increases, decreases or retains (unchanged) the charge 
stored in the loop filter. When Si is closed (and S 2 is opened), the charge- 
pump current flows into the loop filter, increasing the control voltage 

Vp(0 to the VCO. Conversely, when S 2 is closed (this time Si is opened), 7^ 

flows out of the loop filter, decreasing the control voltage of the VCO. When 
both Si and S 2 are opened no current flows in the loop filter and the control 
voltage v^(r) remains unchanged. The magnitude of the charge-pump current 
is an important parameter, as it determines the PPD constant , which 

affects the dynamics of the PLL. The average current flowing into the loop 
filter can be given as [61], [78]: 







(5.26) 



Noting that, by definition, the average output current is the PFD constant, 
Kj , multiplied by the phase error one can easily conclude that [61] 
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(5.27) 

which is measured in A/rad. 




Figure 5-7. Phase-Frequency Detector 

A very simple low-pass filter can be used in the charge-pump loop filter, 
which is made of a series connection of a capacitor C, and a resistor R (see 

Figure 5.8). The resistor in the loop filter is used to realize a zero in the loop 
filter transfer function (which is not located at the origin), so that the 
damping factor ^ of the PLL can be set to a desirable value independently 
from the loop natural frequency w„ . In response to the PFD output, the pump 
current 1^ is driven into (or out of) the loop filter, which responds with an 
instantaneous voltage jump in the control voltage Av,, = Rip . At the end of the 
charging interval, the current switches off and another voltage jump with 
negative sign cancels off the first abrupt change. Although the net effect is 
zero, this causes the VCO to be frequency modulated and is responsible for 
the generation of reference spurs (recall that the comparison in the PFD is 
made approximately at the reference frequency) [61], [78]. To suppress these 
undesirable spurs, a capacitor C can be connected in parallel with the 
resistor. The capacitor has very little effect on the dynamics of the PLL, 
if it is chosen small compared to C, [13]. We will investigate this issue in 
further detail in Section 5.4.1 when we derive the small-signal model for the 
charge-pump PLL. 
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Figure 5-8. A Charge-pump loop filter 

5.3 The Fractional-N Concept 

In frequency synthesizer applications, the main goal is to derive a 
number of periodic signals whose frequency can be externally programmed. 
An important criterion, in addition to high frequency agility and low cost 
implementation, is the signal purity of the synthesized outputs in terms of 
phase noise (jitter) and spurious content [60]. One drawback of the 
traditional PLL frequency synthesizers described above, is that the output 
frequency resolution is the same as the input reference frequency, thus 
imposing some compromises on the phase jitter performance, frequency 
resolution and transient response of the synthesizer [61]. The phase noise of 
the PLL output is multiplied by N within the loop bandwidth, whereas 
outside the loop it follows that of the VCO [64]. On the other hand, the 
switching or settling time of the synthesizer is determined by the inverse of 
the loop bandwidth [64]. Therefore, it is always desirable to use a small 
division ratio N with a large loop bandwidth [60], [61]. In these conditions, 
the desired frequencies can only be achieved with a relatively high input 
frequency, which of course corrupts the output frequency resolution or the 
step size. Thus, in order to provide small step sized moves between adjacent 
channels, a very low reference frequency is required. However, the use of a 
very low reference frequency results in very high division ratio, which 
corrupts the phase noise performance of the PLL [64]. Additionally, the PLL 
loop bandwidth must also be decreased in this condition resulting in long 
settling time [64]. 

The Fractional-N concept has been radically proposed to produce a step 
size finer than the input reference frequency [62]-[64]. The name of the 
technique stems from the fact that the division ratio in the feedback loop is 
not constrained to an integer number only, but additionally, a fractional 
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number is also allowed. There is no frequency divider that can divide by a 
fractional number, but by making the division ratio in the feedback path to 
alternate between predetermined integer numbers in a controlled and 
repetitive fashion, a fractional number is obtained on the average sense. The 
advantages of this new system are obvious. First, for the same frequency 
resolution, a lower division ratio can be used. This results in significant 
reduction in the overall phase noise performance of the PLL. Secondly, as 
the reference frequency is much higher compared to the frequency 
resolution, wider loop bandwidth can be used, which in turn reduces the 
settling time of the PLL. 

5.3.1 First Generation Fractional-N PLL Synthesizers 

The first generation Fractional-N loops, as shown in Figure 5.9, included 
an accumulator whose carry bit was used to vary the division ratio of a dual 
modulus divider between N and N+l. The duty-cycle for the accumulator 
carry bit is always K/F(K and F are integers), where K is the dc input 
(binary word) to the accumulator, and F is the accumulator size. In other 
words, the accumulator carry flag set to “1” K times in every F cycles. 




Figure 5-9. First generation Fractional-N PLL synthesizer 

Referring to Figure 5.9, if K is represented with m bits, then F is given 
by 2* . For F= 16 and K = 3, the fractionality x = K/F is given by 3/16 and 
the dc input to the accumulator is represented by (0011)2. In this case the 
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accumulator will increment to 3, 6, 9, 12, 15, 2 (carry flag set to “1”). 5, 8, 
11, 14, 1 (carry flag set to “1”), 4, 7, 10, 13, 0 (carry flag set to “1”) and so 
on. Note that there are exactly 3 accumulator carry outputs for each 16 
cycles. Each time when the carry flag is set to “1”, the dual modulus divides 
by TV + 1 instead of TV . So the fractionality is obtained by dividing TV + 1 in 
K number of cycles and TV in the remaining F - K cycles [65]. Therefore, 
the output oscillates at a fractional number times the input or 
foul = where X = / F is the fractionality factor. 

However, the first generation Fractional-N synthesis technique suffers 
from an abrupt change in the phase error, thus the process causes spurious 
signals that have to be attenuated by a lower loop bandwidth [62], [63]. 
While the divider is programmed to divide by TV , the phase error begins to 
advance. After a certain time, a carry is activated and the dual modulus 
divider divides by TV + 1 representing a removal of In phase difference. This 
periodic modification results in a saw-tooth phase error, which is responsible 
for the undesirable spurious generation [64]. 

The phase perturbation at the phase detector is entirely predictable and an 
analog compensation signal can be derived from the accumulator content to 
suppress the spurious components [60], [61]. This analog spurious 

suppression technique, however, suffers from implementation imperfections 
because of a D/A converter required to convert the accumulator content into 
voltage. The D/A converter puts more stringent requirements on the 
implementation while increasing the cost and power consumption of the 
overall synthesizer. Some other forms of suppression techniques were also 
developed in the art such as those based on phase interpolation and Wheatley 
random jittering [61]. Both of the methods have their own inherent 
disadvantages. For example, the phase interpolation technique suffers from 
the phase jitter problem during the interpolation process, whereas in the 
Wheatley random jittering approach the resulting phase noise is nearly 
white. Consequently, this kind of suppression techniques are not suitable for 
wireless transceiver applications where low-cost and low-power 
implementations are important while still requiring low-phase-noise and 
spurious free performance [64]. 

As an alternative approach, oversampling noise shaping techniques have 
recently received a lot of attention for spurious suppression in Fractional-N 
synthesizers [62]-[64]. The basic idea of the noise shaping technique is to 
shape the spectrum of the quantization noise such that its power within the 
useful signal band becomes very much smaller. Similarly, in the context of 
the Fractional-N synthesis, the phase error spectrum is shaped, so that the 
phase noise is pushed further away from the carrier frequency, where the 
existing loop filter can easily attenuate them. 
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5.3.2 Higher-order AL Modulation for Modulus Control 

Figure 5.10 illustrates a Fractional-N PLL synthesizer that exploits a 
higher-order AD modulator to control the division ratio of a multi-modulus 
divider in the feedback path. The idea of using a AD modulator to control the 
division ratio in a Fractional-N PLL comes from the fact that an accumulator 
is an equivalent digital implementation of the first-order AD modulator. To 
see this, let us consider a first-order AD modulator with error-feedback 
topology as shown in Figure 5.11(a). This structure is often used in D/A 
applications since the loop can be perfectly implemented in the digital 
domain [11]. The signal-flow graph for this modulator in its digital 
implementations is illustrated in Figure 5.11(b), where the m -bit input signal 
K[n] is summed with the m -bit register content to produce the (m+l)-bit 
quantizer input signal v[n] . The single-bit quantization process is 
accomplished by simply taking the most significant bit (MSB) of v[n] . The 
residual m -bit signal, which represents the negative of the quantization error 
signal, is then stored in an m -bit register to be summed with the input signal 
at the next clock cycle. Figure 5.11(c) shows an accumulator, which can be 
used as the digital implementation of the first-order AD modulator having the 
error-feedback topology of Figure 5.11(a). The accumulator overflow and 
the accumulation result correspond to the single-bit quantizer output and the 
negative of the quantization error at any time, respectively, as suggested by 
the signal-flow graph in Figure 5.11(b). So an accumulator has the same 
first-order noise shaping property as the first-order AD modulator with error- 
feedback topology. 




Figure 5-10. NL controlled Fractional-N PLL synthesizer 
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Figure 5- II. First-order AL modulator with error-feedback topology: (a) block diagram, (b) 
signal-flow graph, (c) its equivalent digital implementation 
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Figure 5-12. A fourth-order all-digital MASH modulator: (a) block diagram, (b) its digital 

implementation 
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Note that in Figure 5.1 1(b) and (c) we denote the input and output K[ri\ 
and b[ri\ as opposed to x[n\ and y[n\ in Figure 5.1 1(a). This convention has 
been used to indicate that the input and output in Figure 5.11(b) and (c) are 
binary words denoting the respective signals. 

In Fractional-N frequency synthesis applications, the divider modulus 
control using a first-order AL modulator produces significant unwanted 
spurious components since lower-order modulators are prone to generating 
limit cycle oscillations and tones when driven by a constant input. As we 
have proven in Chapter 3 and Chapter 4, higher-order AZ modulators (both 
single-stage and multi-stage) are capable of producing tone free outputs 
provided that the first accumulator is started with an irrational initial 
condition. 

Higher-order single-stage modulator, however, require that the amplitude 
of the input signal is limited to half of the quantizer bin width for no- 
overload stability, which in turn reduces the available tuning range in a 
Fractional-N frequency synthesis application. 

Besides providing 100% depth of modulation, a MASH modulator also 
benefits from the advantage of lending itself to simple implementation with 
a minimum of hardware [79]. MASH type modulators have, therefore, 
enjoyed a wide-spread usage in such applications [62], [64], and [80]. A 
fourth-order MASH AX modulator, which is formed by cascading four first- 
order AX modulators configured with the error-feedback topology, is 
presented in Figure 5.12(a). 

Using the fact mentioned earlier on, that a digital accumulator is the 
compact realization of a first-order AX modulator, the digital implementation 
of the fourth-order MASH of Figure 5.12(a) can be simply accomplished by 
cascading four digital accumulators as shown in Figure 5.12(b) [71], [79]. 

The main disadvantage of MASH type modulators is that the 
quantization noise power at the output is related to that of the single-bit 
quantizer in the last stage. Although, the output is multi-bit in MASH 
modulators, this does not result in reduced quantization noise power, unlike 
multi-bit single-stage modulators. 

It is also worth to mention that the fourth-order MASH architecture of 
Figure 5.12(a) generates a 4-bit wide output and therefore a multi-modulus 
divider with 16 programmable division ratios will be required to scan all of 
the desired tuning range. 
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5.4 Design and Simulation of Fractional-N PLL 
Frequency Synthesizers 

5.4.1 Linear Model of the Charge-Pump PLL 

Since the charge-pump loop filter is a time-varying system, it should be 
approximated by an equivalent linear-time-invariant (LTI) system. Our 
analysis will be based on an equivalent LTI system, which has the same 
average current flowing into or out of the loop filter as shown in Figure 5.13. 
Referring to this figure, the voltage-to-average current transfer function for 
the low-pass loop filter can be given as follows (ignoring ): 






(5.28) 



Assume that the loop begins with a phase error = 
then the average current charging or discharging the loop filter is given by 
. The average change in the control voltage can be 

found using (5.28) as 

= (5.29) 



Noting that the output phase is related to control voltage with 
f one can obtain the closed loop transfer function for the 
charge-pump PLL as: 



s + 2^„s + w,; 



(5.30) 



where 



w„ 




and 

^=~wRC, 
^ 2 " ' 



(5.31) 



(5.32) 



Similarly, the phase error transfer function can be found as 
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E{s) = 






(5.33) 



Note that these are the same transfer functions as for the second-order 
PLL with an active loop filter as analyzed in Section 5.2.2. The open-loop 
transfer function for the charge-pump PLL is given by 



0 ( 5 ) = 



\ + sRC, 
N 



(5.34) 



from which it can be deduced that the system is a type-2 second-order 
system. 



/ (t) 




Fi^re 5-13. Loop filter with equivalent average current 

The Bode plot for the closed loop and open-loop transfer functions are 
shown in Figure 5.14(a) and (b), respectively, for different values of 
damping factor | . We see that the second-order PLL is a low-pass filter for 
its input phase signal. The damping factor has an important influence on the 
dynamics of the PLL. As shown in Figure 5.14(a), lowering ^ from its 
optimal value of 0.707 results in overshoots in the response. The degree of 
overshoot is dependent on how small ^ is. Making ^ larger than 0.707 
makes high-frequency suppression worse. For an optimally-flat response ^ 
is usually set to 0.707 so that the closed loop transfer function has a second- 
order Butterworth filter response [78]. 

We choose to ignore in the linear analysis so that the resulting PLL 
would not be a third-order system. Third-order loops cannot be put into a 
standard form, as is the case for the second-order systems. If Cj is chosen to 
be around 1/10 of C, , then the system can be well approximated as a second- 
order system [13], which greatly simplifies the linear analysis. 
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(a) 




(b) 



Figure 5-J4. Bode plot for second-order charge-pump PLL with N =45 for different values 

of^:(a) E(s) 

Nevertheless, the effect of Cj will be included in the simulation to show 
that it does not severely degrade the PLL stability. The inclusion of , in 
fact, introduces a third pole in the PLL transfer function, which is normally 
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at a frequency much higher than the loop resonant frequency. To see this, we 
can write the charge-pump loop filter transfer function including Cj as 
follows: 



F(5) = 



1 

5 C, 



R 

+ - 

sRC^ + 1 



5T; +1 

SC,(5T, + 1) 



(5.35) 



where x^=RC 2 and r, = /?{C, + ) . Using (5.35), the PLL transfer 

function H(s) can be found as 



H{s) = - 



N{sx, + 1 ) 



NC, 



-T, +S 



NC, 



K,K. 



— + JT2 + I 



The open-loop transfer function is now 
+ 1 ) 



(5.36) 



(5.37) 



The transfer function in (5.36) has a third degree denominator and the 
open-loop transfer function in (5.37) has two poles located at the origin. So 
in this case the charge-pump PLL is a third-order type-2 system. Note that 
(5.36) and (5.37) reduce to (5.30) and (5.34), respectively, when ^=0. 
Provided that C 2 /C, is sufficiently small (say between 1/8 and 1/10) [13], 
the addition of a third pole in (5.36) only affects the high frequency filtering 
properties of the PLL. Thus, the low frequency properties of the third-order 
system should nearly be the same as the second-order system. 

5.4.2 Design Issues 

In this subsection, we present the design procedure for a charge-pump 
PLL that can be used in a Fractional-N frequency synthesis application. As 
an example, we have undertaken a design for the GSM900 base station 
receiver systems. 

In GSM900 base station systems there are 175 channels, separated by 
200 kHz from each other. The total bandwidth occupied by this GSM900 
signal is 35 MHz, which must be scanned by a synthesizer. The base station 
has a frequency range of 880 MHz to 915 MHz and 925 MHz to 960 MHz in 
receive and transmit modes, respectively. The reference frequency is set to 
20 MHz, and the settling time for the synthesizer is 10 ps (i.e., the 
synthesizer should switch from one frequency to another within 10 ps). The 
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specifications for GSM900 base station system are summarized in Table 5.1 
[81]. 



Table 5-1. GSM900 base station specifications [81] 



Mode 


Frequency 
range (MHz) 


Channel 

spacing 


Reference 

frequency 


Settling time 


Rx 


880-915 


200 kHz 


20 MHz 


10 


Tx 


925-960 


200 kHz 


20 MHz 


1 0 



Table 5.2 presents the corresponding channel number as well as the value 
of the integer division ratio and the fractional factor X for each carrier 
frequency (output frequency of the synthesizer) in Table 5.1. 

Table 5-2. Channel numbering in GSM900 base station systems 



Mode 


Output 

frequency 

(MHz) 


Channel 

number 


Integer 
division ratio 


Fractionality 
factor (X) 


Rx 


880+0.2n 


0< n <99 


44 


0.0 In 






100<n<174 


45 


0.01*(n-100) 


Tx 


925(0.2n 


0< n <99 


46 


0.01*n-0.25 






100<n<174 


47 


0.01 *(n- 










I00)H).25 



A step-by-step procedure for Fractional-N charge-pump PLL synthesizer 
design is shown in Figure 5.15. Using this, we will now design a frequency 
synthesizer for GSM900 base station receiver systems. 

Step 1: We start by stating the system specifications. For our case study, 
which is the GSM900 base station receiver, the output frequency range is 
880-915 MHz =880 MHz and =915 MHz) in increments of 200 

kHz, where the reference frequency is = 20 MHz. The settling time of the 
Fractional-N PLL should be The charge-pump current is 

usually set according to implementation issues such as power consumption 
and speed [61]. For instance, it is very common to use a charge-pump 
current 1^ = \QfiA [78]. Another system parameter that must be considered is 

the damping factor ^ . We know that the optimal value for the damping 
factor is ^ = 0.707 [78]. The extra capacitor (Cj) in the charge-pump loop 
filter will, however, result in reduced damping factor. As a general rule of 
thumb, the damping factor is set to be 20% higher than the optimal value, to 
compensate for this effect [13]. So a good choice for the damping factor 
would be ^ = 1 . We also need to know the voltage levels and that 
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represent logic “1” and “0”, respectively. For instance, we choose = 5 V 
and t/^j„ = -5 V to represent logic states, which are practical. 



Step 1 



Step 2 



Step 3 



Step 4 



Step 5 




Figure 5-15. Flowchart for Fractional-N PLL frequency synthesizer design 

Step2: In this step, the division ratio in the feedback path of the PLL 
must be determined. In every Fractional-N applications, the ratio of output 
frequency to the reference frequency is variable. Thus, N changes within a 
range ^ . When the division ratio N is variable, the loop natural 
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frequency w^, and damping factor | vary in proportion with as 

dictated by (5.31) and (5.32). In the variable division ratio case, one has to 
define a mean for N , so that w„ and ^ become optimum for this mean 
value. As w„ and ^ depend on N with an inverse square root relationship, 
we make use of the geometric mean of and , which is given by 






( 5 . 38 ) 



As long as the ratio is moderately large (say 10), the variations 

in loop dynamics can be acceptable [78]. If is larger than 10, then 

it is often required to define more than one frequency range, and switch the 
parameters accordingly. 

In our particular application (where y„=880 MHz, ./;„=915 
MHz, = 20 MHz and frequency resolution 200 kHz), = 44 and 
=45.75, hence =44.86. 

Step 3: This step involves the calculation of the PFD constant and the 
VCO constant K^, . As the charge-pump current value is known beforehand 
can be calculated using (5.27). For = the PFD constant is 

= 1.5915x10“^ A/rad. In this step, the characteristic of the VCO should also 
be determined. Let the maximum and minimum output frequencies of the 
VCO be and respectively. Assume that usable range for the 

control voltage v^(/) (that is the input voltage of the VCO) is given by the 
supply voltages and . This means that the VCO was designed to 
generate output frequencies and when its control 

voltage is and , respectively. The (ideal) characteristic of the VCO 
is plotted in Figure 5.16. The gain of the VCO can now be calculated from 
the slope of this characteristic curve as: 



•max ,mm ) 

^max ^mir 



( 5 . 39 ) 



For = 5 V and =-5 V, K„ =2.2x10^ rad*HzA^. 
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Figure 5-16. Input voltage-to-output frequency transfer characteristic of the VCO 

Step 4: Given the settling time of the synthesizer, the loop natural 
frequency can be calculated from [78] 




(5.40) 



After w„ is determined, it is now possible to calculate the value of the 
first capacitor C^ according to (5.31). For 7’^=l0/is we have w„ = 628 
rad*kHz and C, = 2 pF. 

Step 5: The last step is to calculate the second capacitor of the loop filter 
by choosing Cj^C, /lO, which leads to C2=0.2 pF. As a final step, since 

= 1 the value of the resistor R can be found using (5.32) as if = 1610 kU. 
5.4.3 Computer Simulation Model 

Computer simulation of a PLL is a non-trivial task. Although the code 
that must be written is of moderate size, the time-varying nature of the 
charge-pump loop filter makes the simulation very challenging. In this 
subsection, we, therefore, will focus on the most confusing components, 
which are the charge-pump loop filter along with the PFD and VCO models. 
Difference equation simulation of the digital AL modulators is fairly simple 
and straightforward [7], and hence details are omitted here. 

5.4.3.1 State Model of the PFD 

Referring to Figure 5.7, the PFD can be in one of the following four 
states: 

State 1 : Up=0 and Dn=0 
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State 2:Up=0 and Dn=l 
State 3:Up=l and Dn=0 
State 4:Up=l and Dn=l 

However, observe that the last state is practically not possible (i.e., 
ideally it can never happen), by virtue of the AND gate in the feedback path. 
Whenever, the outputs of both flip-flops are logic high, a logic high level 
appears at their reset input clearing the flip-flops, and hence the PFD 
immediately returns to state 1. Consequently, the device can be modeled as a 
tri-state device as illustrated in the state diagram of Figure 5.17. 




Figure 5-17. State diagram model of the PFD 

Some example responses for the PFD outputs are given in Figure 5.18, 
which provide some insights into the operation of the PFD. Consider the 
case where reference frequency is greater than the frequency of the 

divider output. In this case, the PFD generates Up signals with logic high 
state for some time interval^'^, whereas Dn remains always zero. Conversely, 
when is higher than , Dn signal is repeatedly asserted high for a time 

period and Up remains always zero. When , the PFD generates 

pulses at either Up or Dn signals with a width equal to the phase difference 
between its inputs (see Figure 5.18(b)). So, in summary, the PFD produces 
positive or negative pulses to equalize the phase or frequency difference 
between its two inputs. Clearly, Up is used to increase and Dn is used to 
decrease the VCO output frequency. As a final note it is also worth to 
mention that because of the propagation delay time in the AND gate (Figure 
5.7) the PFD may present dead zones, when both of its inputs are triggered. 
This causes the feedback loop of the PLL to be broken (for a short period of 
time) and leads to phasejitters or spurs [78]. We have, however, ignored this 
implementation imperfection in our PFD simulation model. 



16 



Note that in this case, the duty-cycle of the Up signal is determined by the frequency 
difference between the reference signal and the divider output. 
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(b) 



Figure 5-18. Examples for PFD output signals: (a) frequency difference, (b) phase difference 

S.4.3.2 Difference Equation Model of the Loop Filter 

In difference equation simulation of continuous-time systems a fixed 
time step should be chosen, which is the sampling rate of the equivalent 
discrete-time system. This sampling rate should be commensurate to satisfy 
the Nyquist rate for the signal with the highest bandwidth in the system. In 
Fractional-N PLL simulation, however, the simulation step time must be 
chosen very small (say around one-tenth to one-hundredth of the period of 
the PLL output frequency) for an acceptable precise result. By choosing such 
a small T , , it is guaranteed that the changes in voltages and currents are very 
small during the fixed time step. As long as the charge and current changes 
are small, the average current is given by the change in charge divided by 
the time step. So the voltage-to-current transfer function describing the 
charge-pump loop filter can be transformed into an equivalent discrete-time 
system with a voltage-to-charge difference transfer function. 

The reason for transforming the voltage-to-current transfer function into 
a voltage-to-charge difference transfer function is that it is much easier to 
calculate the charge difference during two consecutive time steps rather than 
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the average current. This is the same approach used to simulate switched- 
capacitor systems [13]. 

We define the following transfer functions that describe the charge-pump 
loop filter in the continuous-time and discrete-time domains. 

V (^) 

• F{s) = continuous-time voltage-to-average current transfer 

function 

V (j) 

• continuous-time voltage-to-average charge transfer 

function 

V (z) 

• M(z) - discrete-time voltage-to-average charge transfer function 

y (z) 

• P(z) = discrete-time voltage-to-charge difference transfer 

A^(2) 

function 

Since charge is the integral of the instantaneous current [13], we have 

= (5.41) 

s 



From (5.41) we write 



K(s) = 



Q(s) 



= sF(s) 



(5.42) 



2 l-z~' 

Using the bilinear transformation — T (where E is the sampling 

interval of the discrete-time system) [12], one can write 

M(z)^K(s)\ (5.43) 

1^,-1 



As we are using a sampling interval which is much smaller than (10 to 
100 times smaller) than the Nyquist interval, the commonly encountered 
distortion problem associated with the bilinear transformation need not be 
compensated for. During each simulation time step, we are only interested in 
changes in charge rather than the total charge stored in the system. Noting 
that the charge difference is equal to present charge from which the previous 
value is subtracted, we write the charge difference as. 



A0(z) = (l-z-')0(r) 



(5.44) 
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Using (5.44), we can find the voltage-charge difference transfer function 



as 






V(s) M(z) 



A0(z) \-z- 



(5.45) 



Substituting (5.43) into (5.45), one can obtain 



, 2 l-z"' J 2 1-z"' 

P(z) = r^l 



T l + r‘ T, l + z"' 1-z 



1 



(5.46) 



The charge-pump loop filter can be represented as an equivalent voltage- 
average current transfer function, which is given below (this time including 
C,): 






I^is) jC, 5^C2 + 1 



(5.47) 



2 1-z”' 

By replacing all s variables in (5.47) with = — r» 



2 1-r-') r.(l + z-') ^ 



^r,(i+z'‘) 



T; 1 + z"' I 2C,(l-z"') 2 ^C 2 (l-z"') + r,(l + z-‘) 



Substituting (5.48) into (5.46) results in the following equation: 



(5.48) 









Further manipulation on (5.49) yields, 

j /g(C, + C2) + ^ + z-^^-/e(C,+C2)j 



F(z) 



1 - ?■' T ( T 

/?C,Cj -H ^ C, + z”' I C, - RC,C, 



(5.49) 



(5.50) 



As a matter of convenience, let us define 
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^(C,+C,) + ^ 
D 



(5.51a) 



T. 



m, 



^-/e(C, +Q) 



D 



(5.5!b) 



k = 



c,r,. 

D 



(5.51c) 



where Z) ^ flCiC, +-^C, . With these definitions, 






I m, +z 'wi 

1-7 ' 1+V'(^-I) 



(5.52) 



Hence, the discrete-time voltage-charge difference transfer function can 
be implemented as shown in Figure 5.19. 







-► v.(«] 



Figure 5-19. Discrete-time equivalent model of the charge-pump low-pass filter 

For each time step 7; , the charge difference can be calculated as below 
according to the PFD outputs Up and Dn. 

• if Up=l and Dn=0 => A^[«l = I J, 

• if Up=0 andDn=l Aq[n] = -I^T^ 

• ifUp=0 and Dn=0 => A^[«] = 0 

It is now possible to calculate the control voltage v^. [«] to the VCO 
according to the signal-flow graph in Figure 5.19. 
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S.4.3.3 Difference Equation Model of the VCO 

A VCO produces an output voltage, say , oscillating at a frequency, 

which is a function of its input voltage. Remember that the total phase of the 

/ 

VCO is given by = + JvAtVt , where /(, is the free-running 

frequency. The square wave output voltage of the VCO is then defined by 
= sgn[sin(0,^(O)] (5.53) 

where sgn[-] represents the signum function. 

In computer simulations, however, one has to use a Numerically- 
Controlled Oscillator (NCO) instead of a VCO. Since the discrete-time loop 
filter produces an output voltage v^[n], which is only known at certain time 
instants nT ^ , a NCO can only generate an output voltage sampled at this 
fixed rate. Let us denote the sampled output of the NCO as v^^,[n] . As the 
time step r, is chosen to be very small, we can assume that the output 
voltage of the discrete-time loop filter vjw] does not change during the time 
interval {n - l)r, <t <nT,, Then, the change in total phase of the NCO during 
this time interval is: 

= (5.54) 

This means that if the total phase of the NCO is known at a sampling 
instant (« - IIT, , then it is possible to find the [w] by using 



0„.[«] = 0«>-l] + 2;tfo + ^vVcW7; (5.55) 

Given the value of the total phase, one can directly calculate 
from the following rules: 

• if 2;iA^0^,[«]<(2A + l)?r => = t/,,, 

• if 2;r(^ + l)^0„„,M<2ff(A + 2) => [«] = !/„,„ 

where k is any integer. This comes from the fact, the sign (signum) of a 
sinusoid is positive when its argument is in [o,;r), and negative when it is in 
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5.4.4 Overall Fractional-N PLL Simulation Results 

Figure 5.20 shows a complete simulation model for a Fractional-N PLL 
synthesizer, which includes a 24-bit fixed-point model for the fourth-order 
MASH modulator, as depicted in Figure 5.12(b), as well as a variable 
counter acting as the multi-modulus divider model. To demonstrate the 
functional correctness of this complete model, we have performed various 
simulations. For example, Figure 5.21 shows the spectrum of the PLL output 
for =890 MHz (channel number n = 50), where the integer division ratio 
was set to =44 and the fractionality factor was X = 0.5 . The binary input 
to the modulator was A' = (l00000000000000000000000)j , which corresponds to 
0.5 in a 24-bit representation. Figure 5.21(b) was obtained again with the 
same parameters but this time an initial condition of “1” LSB was set on the 
first accumulator in the MASH modulator. The effect of this irrational initial 
condition in eliminating the spurious components is clearly obvious. The 
only visible spurs are those of the reference frequency (located at a 20 MHz 
distance from the synthesized frequency), whose power is -75 dB with the 
respect to the synthesized signal. In this case, the irrational initial condition 
provided an overall spurious free phase noise floor, which is 10 dB better 
than the uncompensated system. This reduction in the phase noise floor 
comes from the fact that the energy contained in the discrete spurs are 
randomly distributed to the higher frequencies, where the low-pass loop 
filter can attenuate them. It should also be noted that the spurious free 
operation of the synthesizer achieved using the “1” LSB initial condition is 
very important. This is because if any discrete spur is not converted into a 
random phase noise, it can be combined with a strong signal from an 
adjacent channel allocated for another user, resulting in a significant amount 
of interference between message signals [60], [61]. 




Figure 5-20. Complete Fractional-N PLL simulation model 
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(b) 

Figure 5-21. Fractional-N PLL output spectrum ( = 20 MHz, = 44 , X = 0.5 and 
lout - MHz): (a) zero initial condition, (b) “1” LSB initial condition 

Figure 5.22 was obtained for an output frequency of /^, =905 MHz 
(channel number n = 125 ), where this time = 45 and X = 0.25 . The binary 
input to the MASH modulator was A' = (010000000000000000000000)j. Again 
the spurs are eliminated by using the “1” LSB initial condition. Note that in 
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Figure 5.21(b) and 5.22(b), the phase noise of the synthesizer is shaped by 
the high-pass NTF of the MASH modulator multiplied with the closed loop 
transfer function H(s) of the PLL. Therefore, the phase noise increases up to 
a certain point, and then it starts to decrease from where the attenuation 
provided by the closed loop transfer function of the PLL begins to dominate 
the amplification of the NTF of the modulator. 

uu 




(a) 



BO r 




a. 



- 60 f 



. 60 ' ' 1 

860 890 900 910 920 930 
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Figure 5-22. Fractional-N PLL output spectrum ( = 20 MHz, iVj„, = 45 , X = 0.25 and 
= 905 MHz): (a) zero initial condition, (b) “I” LSB initial condition 





FRACTIONAL-N PLL FREQUENCY SYNTHESIZERS 



155 



ai& 

910 

5 905 

s 

5“ 900 
e 

^895 

5 

$ 890 - - - - ----- - 

0 

9 

§885 

1 880 

8 875 
> 

870 

865 

0 10 20 30 40 50 60 70 80 90 

(im« (micro uc) 



(a) 



930 

925 

5 920 

S 

S’ 915 

c 

0) 

a- 910 
£ 

$ 905 - • - 

0 

1 900 

I 895 
O 

u 890 
885 
880 

0 10 20 30 40 50 60 70 80 90 

time (micro sec) 



(b) 



Figure 5-2i. Transient response with “1” LSB initial condition in the MASH AS modulator: 
(a) for Figure 5.21(b), (b) for Figure 5.22(b) 

Figure 5.24(a) and (b) were obtained for an output frequency of 
= 900.2 MHz (channel number n = 101) using zero and “1” LSB initial 
conditions, respectively. The binary input to the modulator was 
if -(010000101000111 101011 101)j . When the binary input to the MASH 
modulator experiences activity near or at the LSB position, the dc value 
resembles an irrational number, and the MASH modulator produces tone 
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free outputs even for the case of zero initial condition. In this case, the “1” 
LSB initial condition is not necessary, but it should be incorporated in any 
Fractional-N design to ensure that the PLL output spectrum is spurious free 
for all input cases. 
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Figure 5-24. Fractional-N PLL output spectrum ( = 20 MHz, = 45 , A” = 0.0 1 and 
= 900.2 MHz): (a) zero initial condition, (b) “1” LSB initial condition 
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Figure 5.25 presents a zoomed-in version of Figure 5.24(b) to show that 
the PLL exactly locks onto the desired frequency for “1” LSB initial 
condition. In other words, the use of the “1” LSB initial condition does not 
result in any noticeable frequency deviations. 




897 898 899 900 901 902 903 904 

Frequency (MHz) 



Figure 5-25. Zoomed-in version of Figure 5.24(b) 

To show that the PLL locks within the specified settling time Tj = IO/l? , 
we investigated the transient response of the loop too. In Figure 5.23(a) and 
(b) we present the VCO instantaneous frequency for the cases of Figure 
5.21(b) and Figure 5.22(b), respectively. As shown in these figures, the 
synthesizer attains lock in less than 10 micro-seconds meeting and 
surpassing our design specification. It should be also noted that the shaped 
phase noise (because of the high-pass MASH NTF and low-pass PLL 
transfer function) is clearly seen in Figures 5.21(b), 5.22(b) and 5.24. The 
quantization noise in higher-order MASH modulators increases very rapidly 
at higher frequencies, and the second-order loop filter can only provide 
limited suppression. If it is desired to further suppress the out-of-band phase 
noise in the PLL output, then a higher-order (for instance a third-order) 
charge-pump loop filter should be used, which in turn makes the PLL a 
fourth-order system [13]. 

5.5 Pipelined Implementation of MASH AZ Modulators 

In this section, a pipelined implementation of a fourth-order MASH AL 
modulator is presented in detail. The reason we choose a fourth-order 
MASH modulator is that apart from providing more suppression of the in- 
band phase noise, a fourth-order MASH also benefits from being less 
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sensitive to the digital approximation of the irrational initial condition when 
fixed-point arithmetic is used. While a third-order structure would be 
simpler, it, on the other hand, suffers from structured spectral components 
present in the output spectrum. For a third-order modulator, increasing the 
dynamic range of the accumulators improves the tracking accuracy, but 
lessens the effect of the “1” LSB initial condition in whitening the 
quantization error sequence (see Chapter 3). This trade-off becomes less 
sensitive for a fourth or higher-order modulator, where the number of bits of 
the accumulators at which the structured components begin to contaminate 
the output spectrum, is beyond the practical limits of the fixed-point 
implementations (i.e., more than 45 and 60 bits for fourth and fifth-order 
modulators, respectively) [82]. 

In the actual implementation of the fourth-order MASH, the limited 
precision of the accumulators results in the loss of accuracy. Therefore, 
before proceeding to the implementation it was necessary to determine the 
dynamic range of the accumulators to furnish the desired accuracy for the 
Fractional-N synthesis. Several fixed-point simulations of the fourth-order 
MASH were performed through Simulink MATLAB™ [83] employing our 
proprietary fixed-point libraries developed for Simulink, and it was observed 
that accumulators with 24-bit dynamic range provided substantially good 
tracking on the desired mean value. Afterwards, a structural bit-level model 
was developed in Simulink MATLAB™, which employs 24-bit 
accumulation. The bit-true model was thoroughly simulated to validate its 
conformance to the desired specifications. 

Figure 5.26 illustrates the circuit realization of the fourth-order MASH of 
Figure 5.12(b). Each accumulator employs a pipelined 24-bit adder and a 24- 
bit register. It can be easily seen that there is a long delay chain between the 
accumulator overflow outputs, and therefore they appear at different time 
instances. To provide synchronization among the carry overflows, they are 
captured in 1-bit registers before being forwarded to the error cancellation 
network. The error cancellation network performs the specific function of 
canceling the quantization noise from the first three stages and produces a 4- 
bit wide output. 

A 24-bit adder could be easily realized by using 24 full-adder (FA) logic 
circuits connected in cascade. However, carry propagation between each FA 
may significantly reduce the operation speed of the overall design. For 
speed-up purposes, we have decided to use six cascaded 4-bit Carry-Look- 
Ahead (CLA) adders in the realization of 24-bit adders, as indicated in 
Figure 5.27. Typically, a CLA adder circuit consists of three sub-systems, 
which are called the propagate generator, the carry generator and the sum 
generator. The logic equation for each sub-system is given below [84]: 
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Pi = -"f/ + yt 
Si={p,*g,)^c, 



(5.56) 



where , _v, , c, and 5 . represent the C bit of the first input, second 
input, carry out and sum out, respectively. Furthermore, • , + and ® 
represent the logic AND, OR and EX-OR gates, respectively, g. is the 
inverse of the g, . 




Figure 5-26. Circuit realization of l-l-l-l MASH 

The carry propagation between each CLA stage results in the critical 
delay path within the 24-bit adder, and thus the maximum achievable 
operating frequency in a given process is limited by the overall computation 
speed of the six 4-bit CLA adders. A 24-bit CLA would not suffer from the 
carry chain problem and may provide faster operation speeds. However, the 
hardware complexity increases exponentially with the number of bits at the 
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CLA adder input. A simple approach to break the long carry chain in adders 
is to employ the pipelining technique so that the carry information is only 
forwarded by one stage [82]. 



c. 




Figure 5-27. Pipelined 24-bit adder formed out of six cascaded 4-bil CLA circuits 
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As shown in Figure 5.27, a 1-bit register is incorporated between each 
CLA stage to pipeline the 24-bit adder. 



pipe shift align shin 




Figure 5-28. A pipelined accumulator topology as proposed in [85] 

It should also be noted that in normal pipelined adders, additional 
registers are required to provide time alignment between the input, delayed 
carry information and the output. The bits of both inputs of the adder should 
be appropriately delayed to synchronize them with the true carry signals. 
Conversely, the output bits are realigned in time by likewise employing the 
appropriate delays. Previously, it was reported in [85], that an accumulator is 
also pipelined in a similar way but without the need of time alignment 
registers for the second input, since it is taken from the delayed output of the 
accumulator. 

A three-bit example of such a pipelined accumulator is shown in Figure 
5.28. Clearly, this topology operates on general time varying inputs. Since 
the AS modulator input is constant for the entire time, the input and output 
alignment registers can be completely eliminated in Fractional-N frequency 
synthesis applications. This is because the magnitude of the input to the 
modulator does not change with time, and hence there is no need at all to 
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store the same input bits and their corresponding output bits in the registers 
for synchronization. Consequently, a pipelined accumulator as shown in 
Figure 5.27 without including the time alignment registers at the input and 
output of the adder, has been successfully utilized in the design of our 
MASH modulator, resulting in considerable savings in area and power. 

Figure 5.29 depicts the logic diagram of the error cancellation network 
used in the realization of the fourth-order MASH of Figure 5.12. The error 
cancellation network basically performs the following function [86]: 

b[n] = CAn]-^C,[n~\] + ^CM~2]-CAn-Y\^C,[n] 

-2Cj[n-l] + C,[fj-2] + C,[t7]-C^[r}-{] + q[n] ^ ‘ ^ 

where Q , Cj , Q and C^ are the carry outs of the first, second, third and 
fourth accumulators, respectively. Basically, there are three stages in the 
error cancellation network. The first stage performs the following operation: 

A[n] = CAn]~^C,[n-\] + 2C,[n-2]-C,[n~3] (5.58) 

This is accomplished by using a 4-bit two-input carry ripple adder and a 
simple logic as shown in Figure 5.29. It is clear from (5.58) that A[n] can 
only vary between -4 and +4 and that is why it is necessary to choose a 4-bit 
adder in the first stage. In the second stage we desire to perform the 
operation A[n] + B[n], where 

B[n] = C,[n] - 2C, [n - 1] + Cd« - 2] + [n] - C,[« - 1 ] (5 .59) 

The result obtained in the previous stage (i.e., A[nj) is directly forwarded 
to the second stage by taking the output summation bits of the first stage as 
the input bits of a 4-bit two-input carry ripple adder in this second stage. The 
third stage is designed to perform one more addition with the C, to produce 
the final output. 

The logic incorporated in the error cancellation network takes four inputs 
and outputs two signals according to the following relation: 

Out\ = In\ • Jn2 • /«4 + In \ • In2 • /n3 • InA + ln \ • /w3 • //i4 

+ In \ • ln2 • /«3 + ln\ • ln2 • Ini • /«4 (5.60) 

Out2 — In\ • Ini + In2 • Ini • In4 + In\ • In2 • In4 + Ini • In4 
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Figure 5-29. Logic design of error cancellation network 

The output of the MASH is a signed multi-bit (4-bit) number, and two’s 
complement binary number representation was chosen since it allows simple 
addition and subtraction. However, the 4-bit two’s complement system 
covers the numbers between -8 and +7, while the 1-1- 1-1 MASH produces 
output levels between -7 and +8. To eliminate one additional bit coding at 
the output, we have coded the output level +8 as -8. The coding table of the 
MASH output indicating the output bits and their corresponding level is 
shown in Table 5.3. 

A VHDL implementation of the architecture of Figure 5.26 has been 
undertaken, and the final design was targeted to an FPGA device provided 
by Altera™ [87] where a Flex lOK (EPF10K10LC84-3) device family was 
chosen to fit the design. The reported maximum achievable operating 
frequency was 34 MHz (this was 13 MHz before the pipelining was 
applied), which provides the desired clocking frequencies for the Fractional- 
N synthesis applications. 
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Table 5-3. Coding table for the MASH output 



Output Level 


Outbit4 


Outbit3 


OutbitZ 


Outbitl 


-7 


1 


0 


0 


1 


-6 


1 


0 


1 


0 


-5 


1 


0 


1 


1 


-4 


1 


1 


1 


0 


-3 


1 


1 


0 


1 


-2 


1 


1 


1 


0 


-1 


1 


1 


1 


1 


0 


0 


0 


0 


0 


1 


0 


0 


0 


1 


2 


0 


0 


1 


0 


3 


0 


0 


1 


1 


4 


0 


1 


0 


0 


5 


0 


1 


0 


1 


6 


0 


1 


1 


0 


7 


0 


1 


1 


I 


8 


1 


0 


0 


0 



The MASH modulator was simulated for all possible values of static 
input word K, and no tones in the whole spectrum were observed. For 
instance, the modulator output noise spectra (dc input removed) for a binary 
input = (010000000000000000000000)2, which corresponds to 0.25 as a 
fraction, is shown in Figure 5 . 30 , where a clock frequency of 20 MHz was 
used. These spectra were obtained by (Hanning) windowed 2 **-point PTT 
data that was normalized to the single-bit quantization step size, which in 
our case was one. Effectively, both of the structural bit-level Simulink 
simulation and FPGA measurement results in Figure 5.30 confirm the 80 dB 
per decade increase in the spectrum, validating the fourth-order noise 

1 Q 

shaping. Furthermore, for 2 samples of modulator output the fraction was 
represented to an accuracy of 99 . 994 %. 
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Figure S~i0. Output quantization noise spectrum of fourth-order MASH: (a) structural bit- 
level Simulink simulation, (b) FPGA measurement 

5.6 Conclusion 

The basic concepts of PLLs are reviewed [78]. We have seen that the 
PLL can be model as a linear system when it is in the lock condition. A 
small-signal model was derived for a PLL from which the closed loop and 
open-loop transfer functions have been derived. The effect of the choice of 
the loop filter on the PLL dynamics has also been investigated. An active 
low-pass loop filter with a very large dc gain is required to achieve a 
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relatively low static phase error [78]. The Op-amp within the active loop 
filter, however, introduces implementation challenges. A simple way for 
achieving low static phase error is the charge-pump technique [61], [78]. 

An overview of the Fractional-N synthesis technique was given. First 
generation implementations suffer from an abrupt change in the phase error, 
thus producing significant amounts of undesirable side-band spurs [60], [65]. 
Several spur reduction techniques for Fractional-N PLL synthesizers were 
also discussed. Phase error compensation by using a D/A converter cannot 
provide a low-cost high-performance solution [60], [62], and [63]. In 
random jittering techniques, the division ratio of the dual modulus divider is 
changed in a pseudo-random fashion. This technique, however, increases the 
phase noise floor of the synthesizer, because the injected phase jitter is 
nearly white [61]. The phase interpolation technique introduces frequency 
and phase jitter in the frequency divider [61]. As an ultimate solution, 
higher-order AS modulator controlled Fractional-N techniques were 
discussed. We have demonstrated that an accumulator is a compact 
realization of an all-digital first-order AS modulator with the error-feedback 
topology [79]. By using this, it was deduced that the first generation 
technique was actually a AS controlled Fractional-N synthesizer, where the 
order of the modulator was one. Better phase noise and spur performances 
can be achieved by using a higher-order modulator, instead of a first-order 
one [62], [63]. 

In this chapter, we have also developed a behavioral simulation model 
for the complete AS controlled Fractional-N PLL frequency synthesizers. 
The simulation model is especially very useful for accurately setting the PLL 
parameters prior to the real implementation. It is also very useful in 
assessing the spurious suppression performance of different AS modulator 
topologies in Fractional-N synthesis applications. For instance, by using this 
simulator, the ultimate effect of the “1” LSB initial condition was clearly 
shown. 

We have also presented the design and FPGA implementation study of a 
fourth-order all-digital MASH AS modulator, which can be used in 
Fractional-N PLL frequency synthesis applications. To achieve the desired 
operation frequency range (20 MHz or higher) while providing low-power 
dissipation and small area, the pipelining technique was utilized in the 
design, where the input and output alignment registers were eliminated by 
taking advantage of the constant nature of modulator input. The computer 
simulation and FPGA experimental results agreed very closely with one 
another, validating the functionality of the developed bit-true simulation 
environment in Simulink MATLAB™ [83] and the physical design for the 
implemented MASH modulator. 
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6. TIME-INTERLEAVED DELTA-SIGMA 
MODULATORS 

In this chapter, a new design methodology utilizing the concept of time- 
interleaved oversampling delta-sigma conversion is developed and explored 
to obtain efficient architectures. In this approach, the time-domain 
expressions of a standard modulator are rearranged according to the desired 
channel count to produce a modular structure with reduced hardware 
requirements. It is shown that the proposed approach results in an 
architecture, which is functionally equivalent to that of the existing method 
based on the block digital filtering concept, but with a reduced hardware 
complexity. The theoretical results are verified by means of behavioral 
simulations. 

We start our development with the modulators formed out of the 
cascaded integrators with distributed feedback topology. We then extend our 
method to cover a more advanced sub-class of modulators containing 
cascaded integrators with weighted feed-forward summation and cascaded 
integrators with distributed feedback as well as feed-forward branch 
topologies. A new time-interleaving concept based on zero-insertion 
interpolation is also proposed, which eliminates the high sampling rate 
multiplexer at the input stage, resulting in a further significant simplification 
in hardware complexity. In this scheme, the input signal is sampled at the 
operation frequency of the channels and applied only to the first channel 
whereas all other channels are fed with zeros all the time. The low-pass filter 
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at the output of the modulator serves two purposes; (i) it rejects the spectral 
replicas of the input signal, and (ii) it attenuates the out-of-band quantization 
noise. 

6.1 Introduction 

As noted before, oversampled AZ modulators have become very popular 
for the design of high-performance A/D and D/A converters in the last 
decade as they comprise very few components compared to conventional 
Nyquist rate converters [5]-[8]. In addition to their reduced hardware 
complexity and inherent high linearity, they also benefit from being robust 
against implementation imperfections such as circuit element inaccuracies 
[14]. A major shortcoming of AZ modulators, however, is the limitation in 
their applications to relatively narrow bandwidth signals due to the 
oversampling requirement [5]-[7]. 

Several parallel AZ approaches have been proposed to date offering the 
potential of extension to high-resolution oversampling noise shaping 
conversion for signals with higher bandwidths. The parallel approaches can 
be classified into three main categories. These are; Frequency Division 
Multiplexing (FDM), Code Division Multiplexing (CDM) and Time 
Division Multiplexing (TDM). In the FDM approach, a bank of AZ 
modulators having different band reject NTF operate on small sub-bands of 
the original input signal [88], [89]. A set of digital band-pass filters then 
attenuate the out-of-band noise at each channel allowing reconstruction of 
the frequency decomposed input signal. This scheme however, may not be 
practical in terms of hardware complexity, since band-pass AZ modulators 
and band-pass filters are required for each channel having different center 
frequencies [88]. Another parallel AZ conversion technique utilizes 
Hadamard modulation to effect a CDM lookalike scheme [90], [91]. In this 
approach, referred to as ITAZ converter, the input signal is multiplied with a 
Hadamard sequence at each channel and quantized using standard low-pass 
AZ modulators. The output of each modulator is filtered using low-pass 
filters, multiplied with the delayed version of the Hadamard sequence and 
finally added together for re-composition [90]. A significant advantage of 
this method is that it includes identical standard low-pass AZ modulation at 
each channel. Another feature of the TIAZ converter comes from the fact that 
the Hadamard modulation sequences consist of ±l’s corresponding to only 
sign changes at the input, hence no multipliers are required. 

Among parallel structures, the TDM scheme or Time-Interleaving (TI) 
provides the simplest form of parallelism [92]. In this approach, the input is 
sampled at a high sampling rate and then distributed through a multiplexer 
over different channels, each operating at a slower rate. For Nyquist rate 
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converters, the TI concept was proposed in the early 80’ s to increase the 
effective sampling rate [93]. However, straightforward application of this 
approach results in only a small increase in SNR performance of the AS 
modulators [94]. Due to the recursive nature of these converters it was 
difficult to formulate the equivalent TI structure for a given modulator till 
the mid 90’s. Fortunately, in 1993 Poorfard et. al. [94], [95] has shown that 
TI AX modulators can be successfully implemented by deploying the block 
digital filtering technique. This approach makes use of M mutually cross- 
coupled AX modulators each operating at a sampling rate of resulting in 
an effective sampling rate of . 

In this chapter, a new design methodology utilizing the time-interleaving 
concept in noise shaping oversampled converters is explored and described. 
In this approach, the time-domain internal node expressions of standard 
modulators are rearranged according to the desired channel counts, to 
produce efficient hardware architectures with modular structures and 
reduced hardware element counts. The TI counterparts of single-stage 
modulators containing Cascaded Integrators with a Distributed Feedback 
(CI-DF) topology are presented along with comparisons to those obtained 
using block digital filtering. It was shown that our novel approach resulted in 
functionally equivalent structures as to those derived from the block digital 
filtering start-point, but with a considerable reduction in the hardware 
complexity. We then extend this approach to a sub-class of modulators, 
which provide more suitable structures for higher-order oversampled noise 
shaping conversion. In fact, the primary focus of this chapter is the 
investigation of a new TI conversion technique based on the zero-insertion 
interpolation technique. In this new approach, the input is sampled at the 
operating frequency of the channels, and applied only to the first channel 
while all other channels are fed with zeros at all time. The readily available 
low-pass filter at the output of the AX modulator serves two purposes; (i) it 
rejects the spectral replicas of the input signal arising form the up-sampling, 
and (ii) it attenuates the out-of-band quantization noise. 

In the next section, the block digital filtering approach as applied to the 
AX modulators is briefly discussed to provide the further details of this 
technique. Section 6.3 covers the proposed approach along with the 
derivation of time-interleaved versions of first, second, and third order 
single-loop modulators for a variety of interleaving numbers (channels). 
Section 6.4 demonstrates some ideal behavioral simulation examples. A 
generic architecture for an l!'' -order M-channel time-interleaved modulator 
is presented in Section 6.5 including a quantified comparison for the 
required hardware of both our presented and the block digital filtering 
approaches. Some practical issues regarding the TI modulator 
implementation such as finite gain of Op-amp, dc offset, mismatch effect 
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and sampling jitter problems are also briefly discussed in a comparative 
manner. In Section 6.6, TI versions of higher-order AS modulator topologies 
are derived. Specifically, Cascaded Integrators with weighted Feed-Forward 
Summation (CI-FFS) and Cascaded Integrators with Distributed Feedback as 
well as Feed-forward Branch (CI-DFFB) topologies are studied and their M- 
channel TI counterparts are presented along with theoretical analysis as well 
as simulation results. The effect of mismatches between the inter-channel 
scaling coefficients is also investigated in this section. In Section 6.7, our 
novel Time-Interleaved conversion technique based on Zero-insertion 
Interpolation is exposed. Verification of this approach is accomplished 
through examples on single-stage modulators with CTDF topology for a 
variety of modulator order ( L ) and interleaving number ( M ). Section 6.8 
discuses some implementation issues regarding the novel Zero-insertion 
Interpolation Time-Interleaving (ZI-TI) concept in a comparative manner to 
the regular TI approach, with the emphasis put on the effects of clock jitter 
and branch mismatches. Finally, we provide some concluding remarks in 
Section 6.9. 

6.2 Block Digital Filtering Approach 

Recently, Poorfard et. al. developed a technique based on block digital 
filtering to realize time-interleaving in oversampling converters [94], [95]. 
To set the scene, in this section we briefly discuss the basics of this 
approach. 

Consider the basic block diagram of a A£ modulator as shown in Figure 
6.1. The equivalent block filtering structure for a given single-input single- 
output (SISO) feed forward loop filter H(z) is presented in Figure 6.2 [94], 
where, H{z) is an MxM transfer function matrix that can be written in the 
form of: 



£■ 0 ( 2 ) 


^l(z) 


E 2 O) ■ 


■ ■ 




^ 0 ( 2 ) 


E,(z) . 


■ ■ Ei,_Az) 






Eo(^) ■ 


• ■ E/^_Az) 


z-'E^(z) 




z-'EAz) . 


■ ■ Eo(z) 



Here, Hij represents the contribution of the input to the output, and 
El's (/ = 0,1,...,M- 1) are poly-phase components of H{z) and can be 
represented mathematically by the following equation [96]: 
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M-l 

/-o 



( 6 . 2 ) 



Since each row is essentially a circularly shifted version of the previous 
row except for the elements below the diagonal entries, H{z) is said to be 
pseudo-circulant [96]. This condition is necessary and sufficient for a block 
digital filter to represent a SISO linear-time-invariant transfer function [95]. 



H(z) (integrator) 







^y\n\ 



Figure 6-1. Basic AI A/D structure 

In the simplest case, the feed forward loop filter is given by the transfer 
function /»(z) = z“'/(l-z"'), resulting in the first order modulator, which is 
also shown in Figure 6.1. The poly-phase decomposition of H(z) for 
channel count M = 2 can be found as: 



l-z"' 



2 -'(z"‘+l) 

l-z-^ 



E,{z^) + z-'E,{2^) 



(6.3) 



and thus 

and £,(z^)=—^- (6.4) 

1-z 1-z 

Hence, 



H{2)^ 




1 1 

2^' 1-Z*' 



(6.5) 
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Figure 6-2. Equivalent time-interleaved structure of a basic modulator obtained using block 

digital filtering 

The corresponding time-interleaved structure for this block digital filter 
is given in Figure 6.3(a). Although the block filter reveals that three 
additional delay elements are needed for realization, one simplification, we 
have made here is to take the first cross-connection from the delayed version 
of u\n] instead of taking it directly from w,[«] (see Figure 6.3(a)). 
Therefore, with the exception of the input multiplexer and output 
demultiplexer, in addition to two explicit integrators, two quantizers, two 
two-input subtractors and two DAC’s, it also needs at least two delay 
elements, two two-input adders and two cross-connections between the 
channels. 

By performing linear analysis on the structure shown in Figure 6.3(a), 
one can easily derive the following input-output transfer function [95]: 

Y{z) - z-^X{z) + (1 - z ' )E{z) (6.6) 

Note that this is the same transfer function as that for the standard A£ 
modulator except with an additional delay in the input-output path, which 
has no effect on the noise shaping characteristic [5]. 
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(a) 




(b) 



Figure 6-3. TI first-order AS modulators derived using the block digital filtering approach: (a) 

two-channel, (b) four-channel 






174 



TIME-INTERLEAVED DELTA-SIGMA MODULATORS 



The hardware requirements for this approach increase very rapidly as the 
channel count M increases. Let us examine again the same feed forward 
loop, but this time for M = 4 with its poly-phase decomposition: 



1-z 



z ‘(z ^ + l + z' 






= E,(z‘) + z-'E,(z') + z-^E,(z') + z-'E,{z‘) ( 6 . 7 ) 



which results in 






,EAz') = ~^,E,(z') = 



-4 » ^^2 



1-Z 



1-Z 



(6.8) 



Hence, the block digital filter of the first order AS modulator for M = 4 is 
given by. 



H{z) = 



I-z' 



z 

z"' 2 

Z“* Z 

z"' z 



1 

1 

z ' 
z"' 



(6.9) 



Figure 6.3(b) depicts the four-channel time-interleaved (M = 4) first- 
order modulator which is derived from the block digital filter given in 
equation (6.9). Again for reducing the number of delay elements the cross- 
connections from the lower to higher channels are taken from the readily 
delayed versions of «. ’s ( i = 1,2,3 ) where they are needed. In this case, with 
the exception of the input multiplexer and output demultiplexer, four explicit 
integrators, four quantizers, four two-input subtractors and four DAC’s, as 
well as at least four delay elements, twelve two-input adders and twelve 
cross-connections between channels are needed. 

The input-output transfer function can be obtained as: 



//(z) = z-^3r(z) + (l-z-')^(z) 



( 6 . 10 ) 



The block digital filtering approach can be easily applied to higher-order 
modulators ( L > 1 ) by taking into account the appropriate SISO feed forward 
loop transfer function [94]. 
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6.3 Time-Interleaved AL Modulators with Reduced 
Complexity 



In this section, alternative to the block digital filtering approach, we will 
show that the time-interleaved counterpart of a AX modulator can be derived 
directly from its time domain behavior. We present all details with examples 
for our approach to achieve efficient architectures for time-interleaved 
oversampling converters. Our approach includes the following steps: (1) 
Express the time domain node equations of the AX modulator for at least M 
consecutive time slots, (2) Re-label the nodes within the modulator in 
accordance with the desired channel count M and modify the time-domain 
equations for each time slot, (3) Combine the resultant time-domain 
equations into a single equation corresponding to one time slot, and (4) 
Derive the desired time-interleaved architecture from the time-domain 
behavior obtained in the previous step. The resulting time-interleaved AX 
converter architectures from our novel approach are modular and result in 
reduced hardware element counts [97]. 

For the sake of clarity, here we perform an example based on the 
simplest case which is the derivation of a two-channel time-interleaved 
(M = 2) version of the first-order AX modulator, where all steps are shown 
for a complete description. 

Consider again the first-order modulator shown in Figure 6.1. Assume 
that the DAC in the feedback loop is ideal, so that y[n] and its analog 
representation can be used interchangeably. Therefore, referring to Figure 
6.1 at the (2n)'* time slot one can write 

v[2«] = u[2n - 1] + v[2« - 1] (6.11a) 

>'[2«] = G(v[2«]) (6.11b) 

u[2n] = x[2n]-y[2n] (6.11c) 

where Q(-) represents the quantization function. Similarly, for the next 



time slot: 

v[2n + 1] = «[2n] + v[2n] (6. 1 2a) 

y[2n + 1] = Q{y[2n + 1]) (6. 1 2b) 

u{2n + 1] = x[2n + 1] - y[2n + 1] (6.1 2c) 

Assume that the input x[n] is obtained from an input multiplexer 
operating at twice the rate of the internal clock frequency of the channels. 
Thus, we shall label the inputs as: 

x,[«] = jc[2«] (6.13a) 

X; [rt] = x[2n - 1] (6. 1 3b) 

Similarly, the nodes of the modulator can be labeled as: 

«,[«] = «[2k] 



(6.14a) 
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Wjlw] - u[2n - 1] 


(6.14b) 


Vi[«]= v[2n] 


(6.14c) 


vAn] = v[2n-\\ 


(6.14d) 


y,[n] = y[2n] 


(6.14e) 


ydn]=^y[2n-l] 


(6.l4f) 


Substituting these newly labeled nodes into equation (6.11) results in the 


following modified time domain equation: 


+ 

II 


(6.15a) 


yM = Q(vy[n]) 


(6.15b) 




(6.15c) 


By further substituting the expressions given in (6.13) and (6.14) into 


(6.12), one may write 


v4rt + l] = i/,[ff]+v,[/?] 


(6.16a) 


>',[h + 1]-0(V2[« + ]J) 


(6.16b) 


U2[n + 1] = A'd/? + IJ- >’2[« + 1] 


(6.16c) 


Eventually, the time-domain behavior of the two-channel time- 


interleaved first-order modulator can be 


found by combining (6.15) and 


(6.16) into one equation as given in (6.17) 


corresponding to only the n'" time 


slot: 


V|[/?]=wd«] + vd«] 


(6.17a) 


vd«] = W|[«-l] + V| [w-1] 


(6.17b) 


y,[«] = 0('^i[«]) 


(6.17c) 




(6.17d) 




(6.17e) 


Wd«] = -«2l«]-.V2[«] 


(6.170 


The corresponding time-interleaved architecture can easily be formed 



using (6.17) as depicted in Figure 6.4(a), where its counterpart derived using 
block digital filtering was already shown in Figure 6.3(a). The input 
multiplexer and output demultiplexer are the same in both structures, 
however the internal blocks are different. In comparing the two structures it 
is convenient to quantify only the hardware required by each parallel 
channel. By carefully examining Figure 6.3(b), it can be found that our 
proposed architecture only needs two quantizers, two two-input subtractors, 
two DAC’s, one delay element, two two-input adders, and four cross- 
connections between channels. 

The proposed converter produces the same transfer function as that of the 
structure in Figure 6.3(a) that was derived using the block digital filtering 
approach. To show this, we once again examine our proposed structure in 
the time domain as follows: the quantizer inputs can be written as: 
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V|[«] = Vj[n] + w 2 [/j] (6.18a) 

V; [n] = V, [« - 1] + M| [n - 1] (6. 1 8b) 

From these we obtain the following results: 

[”] = ^1 [« “ 1 ] + "i [« ” 1 ] + «2 [”] (6 • 1 9 ^) 

V 2 [»] = Vj [k - 1] + M; [«-!]+ u^ [n - 1] (6. 1 9b) 

Similarly, the quantizer inputs of the structure in Figure 6.3(a) can be 
written as: 

w,[rt] = H,[/j-l] + «,[«] (6.20a) 

H’2[n] = M|[n-l] + M2[n-l] (6.20b) 

V, [n] = V| [w - 1] + w, [n] (6.20c) 

Vj [»] = Vj [« - !] + W 2 [«] (6.20d) 

By substituting (6.20a) into (6.20c) and (6.20b) into (6.20d), one can 
obtain: 

V| [«] = V, [n - 1] + U 2 [«] + «[[«- 1] (6.2 1 a) 

[n] = Vj [n - 1] + Uj [n - 1] + w, [n - 1] (6.2 1 b) 



Since the quantizer inputs are the same for both structures, they produce 
the same outputs for the same inputs. The input-output transfer function for 
the proposed structure shown in Figure 6.3(b) can also be easily derived 



from the following equations: 

= / = l ,2 (6.22a) 

U,(z) = X,{z)-YAz), 1 = 1,2 (6.22b) 

K,(z) = t/,(z)4-F,(z), K,(z) = z-'[C/,(z) + f',(z)] (6.22c) 

Y,.(z) = VAz) + E,(z), / = 1.2 (6.22d) 

r(z) = y,(z^)z-'+y,(z') (6.22e) 

E(z) = (z' )z-' + £3 (z^ ) (6.22f) 

Using these expressions the following transfer function can be obtained 
which is the same as the equation (6.6): 

7(z) = z-'Z(z) + (1 - z“‘ )E{z) (6.23) 



In the first application presented so far, the proposed method decreases 
the number of analog delay elements at the expense of a small increase in the 
number of cross-connections between channels, when compared to the block 
digital filtering approach. The increase in the required hardware by the block 
digital filtering approach is proportional to the square of the channel count 
M, whereas the hardware complexity of the proposed architecture is linearly 
related to the channel count [98]. This can be shown by utilizing the block 
diagram of a four-channel time-interleaved first-order AZ modulator 
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depicted in Figure 6.4(b) which is derived using the same procedure outlined 
above. The proposed scheme does solely necessitate four quantizers, four 
two-input subtractors, four DAC’s, one delay element, four two-input adders 
and eight cross-connections between channels. Greater savings through our 
simplifications are now more evident as the interleaving number M 
increases. 

The transfer function of the proposed four-channel time interleaved first- 
order modulator can be found from the following set of equations: 



X,(z) =— i = \,2 M (6.24a) 

M 

V,{z)^X-Xz)-y,{z), i = (6.24b) 

^^.(z) = ^/,,.(z) + ^^.,,( 2 ), / = F„(z) = z-'[t/,(z) + F,(z)] (6.24c) 

y'(z) = f^,{z)+£;,.(z), i = l,2,...,A/ (6.24d) 

r(z) = ]£r,(z*')z'-*' (6.24e) 

(=1 
M 

£(z) = ££,(z“)z'-" (6.24f) 

(=1 



By performing linear analysis on equation (6.24) and after some algebra, 
one may obtain the following transfer function, which is the same as 
equation (6.10): 

H{z) = z^X{z) + (1 - z^' )E{z) (6,25) 

To demonstrate that our proposed method is also applicable for any 
higher-order AE modulator as well, the two-channel time-interleaved 
versions of a second-order modulator as shown in Figure 6.5(a) have been 
explored. 

The resultant architectures of both the block digital filtering and our 
proposed approaches are presented in Figure 6.5(b) and (c), respectively. 
Referring to Figure 6.5(b), in addition to four explicit integrators, two 
quantizers, four two-input subtractors, two DAC’s, the block digital filtering 
approach further needs two delay elements, four two-input adders and four 
cross-connections between channels. Conversely, our proposed architecture 
needs two quantizers, four two-input substractors, two DAC’s, two delay 
elements, four two-input adders and six cross-connections between channels 
(see Figure 6.5(c)). 

Single-stage modulators of order higher than two ( L > 2 ) have some 
stability problems when a single-bit quantization is used [5], [7] and [37]. 
However, such modulators can be made conditionally stable by using the 
appropriate scaling factors in the feedback loops. A third-order AE 
modulator is shown in Figure 6.6(a). 





Figure 6-4. TI first-order AS modulators derived using our proposed approach; (a) two- 

channel, (b) four-channel 
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We have demonstrated that this modulator ensures stable operation up to 
the input signal amplitude of 0.7V, where V is the maximum output level of 
the quantizer [98]. In other words, the DAC in the feedback loop supplies V 
or -V volts if its input is “1” or “-1”, respectively. Using our proposed 
approach, the two-channel time-interleaved version of this-third-order 
AZ modulator has been derived which is demonstrated in Figure 6.6(b). 

It should also be noted that the proposed architectures could be derived in 
the frequency domain utilizing the general block digital filtering theory too 
[97]. To further elaborate on this point, now consider a SISO integrator 
= with its input and output called U(z) and V(z ) , respectively. 

The equivalent block digital filter for this integrator for M = 2 with cross- 
coupling term k is given by [95]; 



H{z)^ 



1 

1 - 2 -' 





( 6 . 26 ) 



Assume that the input and output of the block digital filter is 
U{z) = [Uy{z) UAz)f and V{z) = [V^{z) Viiz)]^ , respectively, thus we may write 
V(z) = H(z)U(z) . Let us assume that H(z) can be decomposed into two 
matrices such as H{z) = (I - B(z)y' A{z) . Hence, V{z) = A(z)U(z) + B{z)V{z) 
where ^( 2 ) and B(z) are called the feed forward and feedback matrices, 
respectively. For the block digital filter implementation of (6.26), A(z) and 
B{z) are as follows: 



A(z)^ 



k-\ 




z~' -k 

7 ^ 




i 

0 



( 6 . 27 ) 



Ideally ^ = 1, and hence A{z) reduces to B(z). If M = 4, then A(z) and 
B(z) are (assuming again ^ = 1 ): 



0 10 0 
0 0 10 
0 0 0 1 
2 "' 0 0 0 



A(z) = B(z) = 



( 6 . 28 ) 
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x|n) — 



quantizer 





Figure 6-5. Second-order AI modulators: (a) standard, and two-channel TI derived using (b) 
block digital filtering, (c) proposed method 



Note that in this special case (i.e., when k = l), A(z) and B(z) are both 
sparse matrices and this is why the complexity of the proposed TI 
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modulators grows linearly with M , rather than M' . However, the proposed 
method may not necessarily result in simpler structures when the cross- 
coupling term is reduced from its nominal value of unity. Reducing the 
cross-coupling term is effective in alleviating the undesirable effects of some 
circuit imperfections in analog implementations [95]. As (6.27) indicates 
when the cross-coupling term is different than unity, A(z) is no longer a 
sparse matrix, and the complexity of the implementation for the proposed 
method is as large as the block digital filtering approach. 

6.4 Simulations 

Let us consider again the first-order modulator shown in Figure 6.1 and 
its four-channel time-interleaved (M = 4) counterparts as depicted in Figure 
6.3(b) and Figure 6.4(b). To validate our theoretical results, simulations were 
carried out for these structures. Throughout the simulations, all modulators 
were excited by a single tone sinusoidal signal at a frequency of 8 KHz with 
amplitude of 0.5. It should also be noted that, in all cases single bit 
quantization has been employed with output levels of ±1 . 

The internal clock frequency of the TI structures and the operation 
frequency of the standard modulator was set to y' = 2.048 MHz. Clearly, for 

TI modulators the input and output multiplexers are operating at the rate of 
four times the internal clock frequency corresponding to an effective 
sampling rate of =8.192 MHz. In all simulations, the in-band SNR was 

computed within the ±16 KHz band (i.e., /^ = 16 KHz, OSR=64 and 
OSRej/=256) frequency range. Therefore, to provide the same resolution per 
bin in the frequency domain, the length of each periodogram was set to 32K 
for the standard modulator, whereas periodograms of 128K length were 
necessary for the four-channel time-interleaved modulators. To prevent limit 
cycle oscillations, the same amount of dithering was applied to each 
quantizer [5]-[7j. The dithering signal was uniformly distributed over an 
interval of [-1/2, 1/2] and directly added to the input of the quantizer as 
shown in Figure 2.10. 

The output spectra for the standard and time-interleaved by four first- 
order modulators are shown in Figure 6.7. These spectra were generated 
using the Welch’s spectrum estimation method, where a Hanning window 
has been applied throughout. 

The tall peak at 8 kHz is of course the desired sinusoid, where the 
remaining components are the shaped quantization noise. Such power 
spectra are often used to numerically calculate the in-band SNR of AL 
modulators. The results in Figure 6.7 reveal that the first-order standard 
modulator has an in-band SNR of 44.6 dB, whereas as anticipated the time- 
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interleaved modulators produces the same output, hence the same SNR 
which is approximately 18 dB better than that of standard modulator 
operating at the same clock frequency. 

It is well known in the art that for a first-order modulator every doubling 
of oversampling ratio (OSR) corresponds to a 9 dB improvement in SNR 
performance [5]. A as a result, the four-channel time-interleaved first-order 
modulators as shown in Figure 6.3(b) and 6.4(b) operate at a frequency of 
/^, but produce approximately the same SNR performance as the standard 
structure of Figure 6.1 clocked at 4/, . 





Figure 6-6. Third-order AZ modulators: (a) standard, (b) two-channel TI derived from the 

proposed method 
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Figure 6-7. Output spectra for standard and four-channel first-order AS modulators 

6.5 Implementation Issues 

An M-channel TI Z,'* -order AL modulator theoretically produces an in- 
band SNR improvement over its standard counterpart for every doubling of 
the interleaving number A/given by the following equation: 

ASNR = (6L + 3) (6.29) 

In D/A conversion, the AZ modulator is implemented digitally, and 
therefore there is no circuit imperfections which might adversely affect the 
performance of the modulator [7], [11]. The proposed TI modulators can be 
implemented perfectly in the digital domain for such applications. However, 
in A/D conversion applications the modulator must be realized in the analog 
domain where some circuit imperfections would degrade the SNR 
performance compared to the expected values given by (6.29). Some of the 
practical limitations considered in this chapter are; finite Op-amp gain, Op- 
amp dc offsets, mismatch effects and sampling clock jitter. These are briefly 
discussed in the following subsection. 

6.5.1 Finite Op-amp Gain 

The finite gain of Op-amp’s (A) generally does not pose any significant 
problem in TI modulators derived from the block digital filtering theory as 
the effective Op-amp gain is Mtimes the nominal gain (i.e., Agj[f^MA) [95]. 
As a result, the TI modulators are more robust against leaky integrators. 
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Unfortunately, the above statement is not correct for TI modulators derived 
using our proposed method. This is because of the cross-terms taken from 
the adders’ outputs. Hence, the effect of the finite gain of the Op-amps is 
similar to the mismatch effect between branches, as aliasing occurs due to 
the discrepancy in the coefficients of the B{z) matrix from their ideal values. 
Figure 6.8(a), (b) and (c) show the simulation results for the standard 
second-order modulator of Figure 6.5(a) and its two-channel TI versions of 
Figure 6.5(b) and (c), respectively, for different values of A. In these 
simulations, the standard modulator was clocked at twice the rate of the 
internal clock frequency of the TT modulators so that the noise floor is the 
same in all structures if A was assumed to be infinite. The simulation results 
reveal that the block digital filtering TI modulator with A=50 dB has nearly 
the same noise floor performance as that of the standard modulator when 
A=56 dB confirming the fact that Aeff=2A. On the other hand, in our 
proposed TI modulator the noise floor is marginally increased when 
compared to the standard modulator for the same finite values of A due to 
the aforementioned aliasing effect. As a result, our proposed structures are 
slightly less robust against the finite Op-amp gain non-ideality and may 
require the use of Op-amps with higher dc gains. 

6.5.2 Op-amp dc Offset 

Op-amp dc offsets normally cause a significant problem in TI modulators 
derived using the block digital filtering approach since the integrator outputs 
may become saturated causing the noise shaping property to be degraded. 
Fortunately, there is a method to alleviate the effect of dc offset to some 
extent in this kind of TI modulators [95]. By reducing the cross-coupling 
coefficient k from its nominal value of unity (1), the degree of overloading 
in quantizer inputs can be substantially reduced. Hence, the block digital 
filtering TI modulators can be made robust to this non-ideality, provided that 
k is carefully chosen. 

On the other hand, in our proposed TI structures the dc offsets of Op- 
amps do not cause any problems as the quantizer input at the first and second 
channel is fed back to the second and first channel, respectively, thereby 
essentially canceling the dc offsets in the respective branches. Our proposed 
TI modulators, therefore, have the capability of inherent overload protection 
at the quantizer inputs in a similar manner to what happens in standard 
modulators. The two-channel TI first-order modulators of Figure 6.3(a) and 
6.4(a) were simulated to validate the above statements assuming ±10 mV dc 
offset and the results of this test are shown in Figure 6.9(a) and (b), 
respectively. Actually, the quantizer inputs for the TI modulator of Figure 
6.3(a) are completely saturated for k = i and therefore the modulator is 
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unstable for this dc offset value. When the cross-coupling term is reduced to 
k = 0.97, the modulator works quite satisfactorily (Figure 6.9(a)). As shown 
in Fig. 6.9(b), our proposed TI modulator, however, produces 
indistinguishable outputs under the ideal and dc offset conditions. 
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Figure 6-8. Output spectra for second-order modulators built out of leaky integrators: (a) 
standard, and two-channel TI derived using (b) block digital filtering, (c) proposed method 

6.5.3 Mismatch Effects 

One of the critical problems in the implementation of TI modulators is 
the mismatch effect between branches. In [99], it was shown that due to 
mismatches, the block digital filtering implementation becomes time varying 
and hence the spectrum component around if^ {i = i,2,...,M- 1) would first 
get attenuated and then folded back into the lower frequencies. The cross- 
coupling term k , which was introduced for handling the dc offset effects, is 
also very effective in correcting the aliased noise floor of TI modulators 
derived employing the block digital filtering approach [95]. The two-channel 
TI first-order modulators were simulated with a 0.1% mismatch ratio. 

The results given in Figure 6.10(a) and (b) show that in both of the 
structures the noise floor is flattened as aliasing occurs. The block digital 
filtering TI modulator, however, produces nearly the same results as the 
ideal case when k = 0.97. Reducing k from its nominal value introduces a 
notch at //, in the noise transfer function so that the noise energy is 
minimized at locations from where the aliased noise is expected to fold back 
into the base band. As discussed in the previous section, the /:-factor method 
was not employed in our proposed structures. If the k-factor method (more 
flexibility) were used in our approach it would result in additional hardware 
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(increase in complexity), which defeats our ultimate goal (simplified 
hardware). 
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Figure 6-9. Output spectra for two-channel first-order modulators with ±10mV Op-amp dc 
offset: (a) block digital filtering, (b) proposed method 
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Figure 6-10. Output spectra for two-channel first-order modulators with 0.1% mismatch: (a) 
block digital filtering, (b) proposed method 
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6.5.4 Sampling Clock Jitter 

In general, any TI A/D converter built out of either Nyquist rate 
converters or AZ modulators is quite sensitive to the input sampling clock 
jitter [100]. This is because of the time uncertainty of the stored value in the 
S/H devices preceding the A/D converters. To overcome this problem, a two 
rank S/H structure can be used, in which the input signal is held constant by 
the first S/H device and then distributed to each channel by M S/H devices 
each operating with 2nl M phase difference [101]. 

6.5.5 Critical Delay Problem 

At a first glance, the TI modulators derived using the block digital 
filtering approach seem to have M times the critical delay path when 
compared to the standard structure [95]. The same statement seems also to 
be true for our proposed TI modulators. In fact, both of the TI modulators 
need the usage of the clock staggering strategy in the analog implementation 
or pipelining techniques in the digital implementation to fully benefit from 
the parallel processing propriety offered by the approaches. Otherwise, the 
maximum achievable clock frequency in TI modulators would be M times 
slower than that of the standard modulator rendering the approaches 
infeasible. 

A first-order AS modulator with error-feedback topology is shown in 
Figure 6.11(a), which contains no explicit integrators. Its two-channel TI 
version derived using the proposed method is also given in panel (b) of this 
figure. The critical delay path is defined as the signal path, which includes 
no delay elements. The block diagram of Figure 6.11(a) reveals that in the 
standard modulator the critical delay path consists of one multi-bit adder and 
one single-bit subtractor. In digital implementations, however, the single-bit 
quantization and single-bit subtractor are embedded in the overflow 
mechanism of the multi-bit adder (see Figure 5.15), and hence there is no 
physical quantizer and single-bit subtractor required in the implementation. 
In reality, the critical delay would be determined by the computation speed 
of one multi-bit adder (accumulator). In other words, the calculation of u[ri\ 
should be accomplished before the next clock cycle. In the case of the time- 
interleaved modulator, the critical delay is determined by the total 
computation delay of two multi-bit adders. Obviously, the required 
computation speed is doubled in the time-interleaved case, which contradicts 
with the aim of the parallelism. A simple solution to this is to utilize 
pipelined adders in the implementation. 
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Figure 6-11. First-order A1 modulators with error feedback topology: (a) standard, (b) two- 
channel TI derived using our proposed approach 

In pipelined operations, additional single-bit registers are inserted in the 
carry chain of the adders to reduce the critical delay path. Since each carry 
information is stored in a register, one adder can start its calculation as soon 
as the LSB of the result from the other adder is ready In this case, the 
critical delay path is determined by the time each full-adder unit in the 
adders needs to complete its calculation. Assuming that both of the standard 
and two-channel TI modulators are implemented using pipelined adders, we 
may conclude that the critical delay path would be the same in both 
architectures. Therefore, with appropriate exploitation of pipelining 
technique the TI modulator can fully benefit from its parallel operation to 
increase the effective oversampling ratio. 

This argument is based on assuming that adders are composed of full-adder units, i.e. carry- 
ripple adders. However, the same reasoning can be extended to the carry-save adders or 
carry-lookahead logic adders. 
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6.6 Hardware Comparison 

The approach presented in this chapter can be used to derive the time- 
interleaved version of any standard modulator. The technique can be 
systematically applied to obtain a generic time-interleaved structure. An U*' - 
order modulator and its M-channel time-interleaved counterpart are shown in 
Figure 6.12(a) and (b), respectively. 

The proposed TI architecture is repetitive, thus has a modular structure. 
Referring to Figure 6.12(b), the internal nodes of the modulator are formed 
by using appropriate cross-connections form a higher channel, except at the 
last channel where the cross-connections should be taken from the first 



channel via delay elements. 

Specifically, by investigating Figure 6.12(b) one may write 
V,jn]-U,jn] + Vi„^,,[n] for ( / = - 1 and m = l,2 M-\) (6.30a) 

[«1 foi" ( w - 1 , 2 ,.,., M -\) (6.30b) 

However, at the A/* channel 

mIwI+v,,!!"-!] for(/ = o,l....,l-l ) (6.31a) 

[«] ^ [«-!]+ I'/.j [« - 1] (6.3 1 b) 



The structure involves a number of delay elements which is equal to the 
order of the modulator L, where (L -1) of them are embedded in the cross- 
connections taken from the first channel to last channel and the remaining is 
incorporated on the last stage of the last channel. 

By using the block digital filtering approach [94], [95], one can also 
derive the M-channel time-interleaved version of the U -order modulator 
shown in Figure 6.12(a). The resulting modulator is very cumbersome and 
therefore not shown here, but its hardware requirements are quantified for 
comparison. 





(b) 

Figure 6-I2. F’' -order AS modulator structures formed from the CI-DF topology: (a) 
standard, (b) A/-channel proposed TI 

Table 6.1 summarizes the hardware complexity for M-channel time- 
interleaved structures of both approaches. The number of quantizers, two- 
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input subtracters and DAC’s are the same for both of structures and equal to 
M, LM and M, respectively. The block digital filtering approach necessitates 
LM explicit integrators, - M) cross-connections between channels, M 
delay elements and L{M^ - M) two-input adders. In contrast, our proposed 
structure uses M{L -i-l) cross-connections between channels, L delay 
elements and LM two-input adders. 

Table 6-1. Comparison of the block digital filtering and our proposed approaches in terms of 

hardware requirements 





Block digital filtering 


Proposed method 


Input multiplexer 


needed 


needed 


Output demultiplexer 


needed 


needed 


U of quantizer 


M 


M 


# of two-input subtractor 


LM 


LM 


# of DAC’s 


M 


M 


# of explicit integrators 


LM 


none 


# of cross-connections 


L{M^ - M) 


(/. + 1)A/ 


# ofdelay elements 


M 


L 


U of two-input adders 


L{M'^ -M) 


LM 



6.7 Higher-order Time-Interleaved AS Modulators 

Higher-order AZ modulators offer improved SNR performance for a 
given OSR by pushing more of the noise power to frequencies outside the 
signal band. They are also less prone to generating idle tones and limit 
cycles. A major obstacle, however, is the stability issue, especially when a 
single-bit quantizer is used. Oversampling noise shaping functions of order 
higher than two consisting of pure differentiation result in inherently 
unstable systems [7] (with a single-bit quantizer). Consequently, higher- 
order single-stage modulators are realized by making use of some feedback 
scaling factors, which reduces the dynamic range of the modulator. Previous 
sections dealt with modulators formed out of the Cascaded Integrators with 
Distributed Feedback (CI-DF) topology. 

As an alternative to the CI-DF structure shown in Figure 6.12(a), some 
stable higher-order single-stage topologies have also been proposed in the 
literature which employ high-pass Butterworth or inverse Chebychev NTF 
[36], [37]. In this section, we extend our TI approach to a more general class 
of single-stage modulators by covering two different topologies as suggested 
by Adams in [38]: Cascaded Integrators with weighted Feed-Forward 
Summations (CI-FFS) and Cascaded Integrators with Distributed Feedback 
as well as Feed-forward Branch (CTDFFB). 
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6.7.1 Cascaded Integrators with Feed-Forward Summation 
Topology 

To provide further details of our approach, here we perform an example 
on the derivation of the two-channel TI counterpart of a second-order single- 
stage AS modulator with CI-FFS topology. 

Derivation of the TI structure: 

Step (1): Consider the second-order modulator formed from CI-FFS 
topology as shown in Figure 6.13(a). Assuming that the DAC in the 
feedback loop is ideal (i.e., its transfer function is 1), we may write the time- 
domain node equations of this modulator at the (2«)'^ and (2« + 1)"' time slots 



as follows: 

V, [In] = u[2n - 1 ] + v, [2n - 1 ] (6.32a) 

v2[2rt] = »j[[2/7-l] + V2[2rt-l] (6.32b) 

r[2n] = b{\>\2n]->r bjViiln] (6.32c) 

y[2n] = Q{r[2n]) (6.32d) 

u[2n] = x[2n] - y[2n] (6.32c) 

and, 

V| [2n + 1] = u[2n] + v, [2n] (6.33a) 

v 2[2« + 1] = vT2n] + V2[2n] (6.33b) 

r[2« + 1] = b^v^[2/l + l] + ^jjV 2 [ 2 n + 1] (6.33c) 

y[2« + l] = Q(r[2n + l]) (6.33d) 

u[2n + 1] = x[2n -f- 1] - y[2n + 1] (6.33e) 



where Q(-) represents the quantization function. 

Step (2): Assume that the input x[n] is being distributed over two 
channels through an input multiplexer operating at twice the rate of the 
internal clock frequency of each channel. Thus, we shall label the inputs of 
each channel as follows: 

X|[n] = x[2n], = (6.34) 



Similarly, the internal nodes of the modulator can be re-labeled as: 



v, ,[rt] = vi[2n], v,,Jfl] = Vi[ 2 «-I] (6.35a) 

V2.M=V2[2n], V3j[«] = V3[2n-l] (6.35b) 

^i[«] = d2«], rAn] = r[ 2 n-\] (6.35c) 

= w[2ff] , «2[«] = «[2n-l] (6.35d) 

yM = y[2n], yM = y[2n-\] (6.35e) 
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Substituting these newly labeled input and internal node expressions 
given in equations (6.34) and (6.35), into equation (6.32) results in the 



following modified time-domain equation: 

\i[«] = « 2 [«] + ^i. 2 [”] (6.36a) 

^ 2,1 [«] = v ,,2 [«] + vj 2 [rt] (6.36b) 

= [«]-!- (6.36c) 

y\n] = Q{rM) (6.36d) 

«,[«] = ;c,[«]-yi[«] (6.36e) 

By further substituting the same expressions into equation (6.33), one 
may write 

v,_ 2 [n + l] = W|[«] + V|_|[«] (6.37a) 

" 2,2 [« + !] = v,_| [«] + Vj , [n] (6.3 7b) 

r^[n + 1 ] = 6 ,v, Jrt + 1 ] + 62 V 3 2 [« + 1 ] (6.37c) 

+ = + (6.37d) 

« 2 [« + 1 ] = JCjlrt + l]-^ 2 [w + 1 ] (6.37e) 

Step (3): The time-domain behavior of the two-channel TI counterpart 
can be found by combining equations (6.36) and (6.37) into one equation 
corresponding to only the time slot as given in equation (6.38). 

V|,i[«] = « 2 W-*-^i, 2 [«] (6.38a) 

V| 2 [n] = m,[h - 1] + V, |[n - 1] (6.38b) 

V 2 ,|[«] = ^l, 2 [«] + *’ 2 , 2 [«] (6.38c) 

V2.2 [«] = V,.| [»-!] + Vj,, [n - 1] (6.3 8d) 

[«] = [»] + h '' 2.1 [«] (6.38e) 

f‘2 [«] = i» 2 V 2.2 [«] (6.38f) 

yM = Qinm (6-38g) 

y2{n] = QU2[n]) (6.38h) 

u^[n] = x^[n]-y^[n] (6.38i) 

«;[«] = Xj[n]-y 2 [«] (6.38j) 



Step (4): The desired TI modulator can be easily formed from equation 
(6.38) as illustrated in Figure 6.13(b). It is worth noting that the outlined 
procedure has resulted in a structure where there is no explicitly visible 
integration block. The accumulation process is, however, accomplished 
through the cross-connections between channels and the add-and-delay 
circuitry embedded on the second channel. Consequently and as expected, 
the TI modulator produces the same transfer function as that of the standard 
modulator. 
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(b) 

Figure 6-JS. Second-order A£ modul^or structures formed from the CI-FFS topology: (a) 
standard, (b) two-channel proposed TI 

Assuming a linear model for the quantizer (i.e., additive white noise 
approximation) and performing linear analysis on the structure, the input- 
output transfer function of the TI second-order CI-FFS topology modulator 
given in Figure 6.13(b) can be found as: 



Y(z) = z-'X(z)S(z) + E{z)N{z) 



( 6 . 39 ) 
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where N{z) and 5(z) are the NTF and STF, respectively, and given by: 



N(z) = 



1 



! + 



z-1 



+ 



U-\f 



S{z) = \-N{z) 



(6.40) 



Hence one can identify from (6.39) and (6.40), the XI modulator shown 
in Figure 6.13(b) produces the same transfer function as the standard 
modulator of Figure 6.13(a), with an additional delay incurred on the input- 
output path due to the increased parallelism [38], [102]. 

By systematically applying the method outlined here for different 
modulator order (L) and interleaving number (M), one can come up with 
generalized XI structures with an arbitrary interleaving number M as it 
comprises repetitive structures. 

Xhe I'* -order CI-FFS topology modulator is illustrated in Figure 6.14(a) 
where its M-channel XI counterpart is depicted in Figure 6.14(b). 

6.7.2 Cascaded Integrators with Distributed Feedback as well as 
Feed-forward Branch Topology 

For the CI-FFS topology when the NXF or SXF is specified, the other is 
fixed automatically as dictated by equation (6.40). Xhis property of the CI- 
FFS topology is restrictive, as one is unable to set NXF and SXF 
independently. In order to allow a certain degree of freedom in setting NXF 
and SXF, we also studied the CI-DFFB topology. 




(a) 
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The L'* -order modulator formed with CI-DPTB topology and its M- 
channel TI counterpart are shown in Figure 6.15(a) and (b), respectively. 
Again by performing linear analysis, the transfer function of the M-channel 
TI modulator ofFigure 6.15(b) can be formulated as: 



Y{z) = z^-'X{z)S{z) + E{z)N(z) 
where, 






(^- 1 )^ 



(6.41) 



(6.42a) 



S{2) = 



b^ + Z>2(z“0 + -" + ' 

(z-1)^ +«3, +a 2 (z-\) + ... + ai^(z-])^~^ 



(6.42b) 



This is the same transfer function as that of the standard Z,'* -order CI- 
DFFB modulator ofFigure 6.15(a) with the additional delays appearing on 
the input-output path [38], [102]. Note that this linear analysis was aided 
with the use of the Matlab symbolic equation handling tool [103]. 

6.7.3 Simulations 

The theoretical results are verified through simulations performed on the 
four-channel TI versions of fourth and third-order modulators formed from 
the CI-FFS and CI-DFFB topologies, respectively, with their respective 
coefficients listed in Table 6.2. 

Table 6-2. List of coefflciente for higher-order modulators 



Topology 


Coefficients 


fourth-order CI-FFS 


b 1=0.806538 




b2=0.307230 




b3=0.063500 




64=0.005620 


third-order Cl-DFFB 


a 1=61=0.044083 




a2=b2=0.287163 




a3-b3-0.800431 



These coefficients were synthesized by specifying a high-pass 
Butterworth NTF with a pass-band gain of nearly 3 dB as explained in [38]. 
Throughout the simulations, a single tone sinusoidal excitation signal at a 
frequency of 8 KHz with an amplitude of 0. 1 was applied to all modulators. 
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with the modulator clock frequency set to 1.024 MHz (i.e., /, = 1.024 MHz). 
For the four-channel interleaved modulators, the effective sampling rate 
is four times the modulator clock frequency (i.e., =4.096 MHz). The 

in-band SNR was again computed over a ±16 KHz frequency range (i.e., 
/ft = 16 KHz) corresponding to OSR=32 and OSRejf=\2S for the standard and 
four-channel TI modulators, respectively. 
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Figure 6-J6. Output spectra for higher-order modulators (standard vs. four-channel 
proposed TI): (a) fourth-order CI-FFS, (b) third order CI-DFFB 
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Figure 6.16(a) illustrates the output spectra of the fourth-order standard 
modulator having the CI-FFS topology with the result for its four-channel TI 
counterpart overlaid on the same plot. These spectra were generated using 
the Fast Fourier Transformation (FFT) with a Harnning window applied 
throughout where 16 50% overlapping periodogram averages were used for 
spectral estimation. 

The simulation results in Figure 6.16(a) reveal that the standard fourth- 
order CI-FFS modulator achieves an in-band SNR performance of 50.13 dB 
at an OSR=32. Its four-channel TI counterpart produces an in-band SNR, 
which is approximately 54 dB better than that of the standard modulator 
operating at the same clock frequency. For a fourth-order modulator, the 
SNR performance increases by 27 dB for every octave increment in the OSR 
[5]. In conclusion, despite the fact that the four-channel TI structure was 
clocked at a frequency of for each of its channels, its SNR performance 
came very close to that of its standard counterpart being clocked at 4/^ . 

The same numerical experiment was repeated for the third-order 
modulator with the CTDFFB topology. The simulation results presented in 
Figure 6.16(b) demonstrate that the four-channel TI third-order modulator 
generates approximately 42 dB better in-band SNR performance in 
comparison to its standard counterpart, which corresponds to the same SNR 
enhancement for a third-order modulator when the sampling rate was 
quadrupled. 

It is well kwon in the art that higher-order AS modulators are 
advantageous especially when a high OSR is used [7], [38]. The simulation 
results presented here indicate that there is only approximately 3.5 dB 
difference in the SNR performances of the fourth-order CI-FFS and third- 
order CI-DFFB modulators at OSR=32. However, when the effective 
sampling rate is increased by four through interleaving (i.e., 05i?g^l28), the 
benefit of higher-order noise shaping becomes more noticeable as the fourth- 
order CI-FFS modulator provides approximately 16 dB better in-band SNR 
performance than the third-order CTDFFB modulator. 

6.7.4 Coefficient Mismatches 

A major problem that arises in the analog implementation of TI 
modulators, is the issue of mismatches between the inter-channel scaling 
coefficients. Generally speaking, any TT modulator is a multi-rate system 
and its successful operation relies on perfect aliasing cancellation [99]. In 
other words, the aliased components generated due to the down-samplers 
will be completely canceled in the final output if the system is ideal. 
However, because of the analog circuit element inaccuracies the scaling 
coefficients will always be mismatched among various channels and aliasing 
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occurs, since the overall structure becomes time varying. We have 
performed here similar simulations to those reported in the previous 
subsection, in order to investigate the effect of scaling coefficient 
mismatches in the CI-FFS and CI-DFFB topologies. 

In the CI-FFS topology, the quantizer input at each channel is obtained 
by multiplying the respective adder outputs with scaling coefficients and 
summing them. Therefore a small deviation in the scaling factors among the 
channels does not pose any problem, as the sign of the quantizer inputs 
remains intact. In actual fact, the four-channel TI fourth-order CI-FFS 
modulator was simulated with up to 10% scaling coefficient mismatches and 
the measured in-band SNR was as expected within 1 dB of the ideal 
performance. Hence, we conclude that this topology is robust against this 
kind of mismatches when a single-bit quantizer is used. However, the above 
statement is not at all true for the CI-DFFB topology because of the 
feedback scaling factors. The four-channel TI third-order CI-DFFB 
modulator was simulated with 0.1% mismatches in scaling coefficients and 
the result of Figure 6.18(a) shows that the noise floor is raised and flattened 
due to aliasing and the in-band SNR performance is deteriorated by 
approximately 22 dB compared to the ideal case. 

Fortunately, there is a way to reduce the amount of aliased noise in the 
baseband to a great extent. Poorfard et. al. [99] has shown that any 
mismatches among the channels causes the frequency portions around 
InijM {i = 1,2,...,M-1) to be folded back to the baseband, and in order to 
alleviate this effect further suggested to modify the noise shaping function in 
such a way that the noise energy is minimized around ImjM . This can be 
easily accomplished by introducing notches in the NTF of the modulator at 
iTtijM . Replacing all the three integrators within the third-order CI-DFFB 

modulator with the transfer function of h{z) = (z"^ +/b"‘)/(l - z~^) and applying 
the proposed method results in the TI structure shown in Figure 6.17, where 
the block filter H{z) is given by: 



z-' 


k 


k 


k 


Az"' 


z“' 


k 


k 


kz~' 


kz~' 


z ' 


k 


kz^' 


kz~' 


kz-' 


z~' 



( 6 . 43 ) 



Our method, however, results in the same structure as the block digital 
filtering approach proposed by Poorfard et. al, indicating the fact that the 
method may not necessarily give simpler TI architectures. In section 6.3, it 
was shown that our method corresponds to a sparse matrix decomposition of 
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the block digital filter and such decomposition for any arbitrary block digital 
filter may not necessarily exist. 




Figure 6-1 7. Four-channel TI third-order block digital filtering CI-DFFB modulator 

The four-channel TI block digital filtering CI-DFFB modulator with 
0.1% mismatches in the scaling coefficients was simulated by gradually 
reducing k from its nominal value of unity until the aliased noise is not 
noticeable. Figure 6.18(a) also shows the result for this modulator when 
k = 0.78 . The in-band SNR performance is improved by 11 dB compared to 
our proposed structure with mismatches. However, while for high OSR the 
^-factor method is very effective in correcting the aliased noise, for low 
values of OSR the in-band SNR performance may become worse. Reducing 
the cross-coupling term k introduces two notches in the noise transfer 
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function at half-Nyquist and Nyquist frequencies of the modulator. These 
notches increase the overall noise floor of the modulator, as the total power 
is constant at the output [95]. For instance, Figure 6.18(b) depicts the results 
for the four-channel TI third-order CI-DFFB modulators where OSR=S and 
OSRef=32. Since the in-band quantization noise power is dominant in 
comparison to the aliased noise, the performance of the proposed TI 
modulator under these conditions is now better than that of the block digital 
filtering approach with k = 0.78 , in the presence of aliasing. 

6.8 Zero-Insertion Interpolation Time-Interleaving 

So far we have introduced a method for the architectural desig of TI AZ 
modulators and extended our work by showing that it can be successfully 
applied to different topologies. The resulting TI modulators do not employ 
any explicit integration blocks. All the accumulation processes are achieved 
at the A/* channel through the use of L add-and-delay elements. All other 
channels consist of only summation blocks, and therefore simplify 
processing requirements. Rather than utilizing M cross-coupled AZ 
modulator in parallel, our method suggests a way to share the already 
existing resources, and successfully reduces the complexity of the whole 
implementation. 

Although the resultant TI modulators are already efficient in terms of 
hardware complexity, the elimination of the high sampling rate input 
multiplexer would be advantageous. In an effort to reduce the hardware 
complexity involved in the implementation of the high sampling rate input 
multiplexer, we have explored in this section a new TI concept based on our 
novel application of the zero-insertion interpolation technique. The 
underlying idea here, is that the input signal is sampled at the operating 
frequency of the individual channels A (of course, provided that /, respects 
the Nyquist criterion), and applied only to the first channel whereas the 
remaining channels are being fed with zeros at all time. 

The output spectrum consists of the shaped quantization noise and the 
original input signal along with its spectral replicas arising from the up- 
sampling. The low-pass filter at the output of the AZ modulator suppresses 
the high frequency quantization noise as well as the spectral replicas of the 
up-sampled input. 

The way we distribute the input signal for the conventional TI and Zero- 
insertion Interpolation Time-Interleaved (ZI-TI) approaches are shown 
schematically in Figure 6.19(a) and (b), respectively. In the regular TI 
approach, the input signal is sampled at the effective sampling rate of the TI 
modulator and every sample is applied in a commuted manner to 

each of the channels operating at A • The outputs from each of the M 
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channels operating at /, are recombined as seen in Figure 6.19(a) and the 
ultimate modulator output >?[n] streams through at an effective rate of Mfs. 
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Figure 6-18. Output spectra for four-channel Tl third-order CI-DFFB modulators with 0.1% 
scaling coefficient mismatch (block digital filtering with k = 0.78 vs. proposed): (a) 
05/e=32, OSR^f^n^, (b) OSy?=8, OSR,Jr^2 




208 



TIME-INTERLEAVED DELTA-SIGMA MODULATORS 





(b) 



Figure 6-19. Input signal distribution for interleaving approaches; (a) regular TI. (b) ZI-Tl 

In the ZI-TI concept, the high sampling rate input multiplexer is 
completely avoided as the input signal needs only be sampled at the 
operating frequency of each individual channel / . A gain factor of M is 

incorporated, which scales the input signal, to compensate the input signal 
power loss at the output of the ZI-TI modulator due to the up-sampling, and 
will be explained in greater detail in the following paragraphs. 

The signal processing scheme for the ZI-TI approach is illustrated in 
Figure 6.20. As can be seen from this figure, the major difference between 
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the spectral characteristics ofZI-TI and regular TI scheme is the presence of 
the input signal replicas spaced at intervals, which were generated as a 
result of the zero-insertion up-sampling. However, these replicas are well 
below the modulator shaped quantization noise characteristic and lie outside 
the base-band region of the modulator. As a result, the decimating low-pass 
filtering (or smoothing low-pass filtering in D/A applications) which is an 
essential element of the system, would conveniently attenuate and filter out 
these noise hurried replicas as well as the out-of-band quantization noise. 

A standard first-order modulator and its two-channel TI version, which 
was derived using the method outlined in Section 6.3, are shown in Figure 
6.1 and 6.4(a), respectively. The transfer function of the TI modulator of 
Figure 6.4(a) is given by: 

Y{z) = z-^X{z) + (1 - z ‘‘ )E{z) (6.44) 

The two-channel ZI-TI version of the first-order modulator is depicted in 
Figure 6.21. Let us assume that there is a unity gain factor at the input. With 
this assumption, the transfer function for the ZI-TI modulator can be easily 
derived as: 

y(z) = L-^x(z)+o-z-')EU) (6.45) 

The analytically derived expressions given in (6.45) shows that the input 
signal component at the output of the ZI-TI modulator is halved, which 
results in a 6 dB SNR degradation in comparison to the equivalent TI 
structure. 




Figure 6-20. Signal processing scheme for ZI-TI approach 




210 



TIME-INTERLEAVED DELTA-SIGMA MODULATORS 



Generally, for an M-channel ZI-TI the SNR degradation is given by 
lOlog(M^) . To eliminate this power loss due to zero-insertion up-sampling of 
the input signal, a gain factor of Mis included at the input stage of the ZI-TI 
modulator. At a first glance, this gain factor may be deemed to be causing 
overloads at the quantizer inputs. However, a close investigation reveals that 
the internal nodes of the TI and ZI-TI modulators have nearly identical 
statistical characteristics such as mean, variance and maximum values. 
Therefore the ZI-TI approach with input gain factor would conveniently 
utilize the same Op-amp characteristic in terms of signal swings as the 
regular TI structures [102]. 

Regarding the ZI-TI concept, the work we have reported so far has been 
backed by theoretical analysis for the first-order modulator. Due to the 
considerable complexity involved in the algebra of higher-order TI 
modulators with an interleaving number M greater than 2, we have opted to 
deploy the same input gain factor principle for the M path ZI-TI modulators 
as was suggested by the analytical first order results, without resorting to 
cumbersome algebra. 




Figure 6-21. First-order two-channel proposed ZI-TI AL modulator 

The structures illustrated in Figure 6.22(a) and (b) are those for the two- 
channel ZI-TI second-order modulators derived out of the proposed method 
and block digital filter, respectively. The TI versions of these modulators can 
be easily reconstructed from their corresponding ZI-TI structures by 
replacing the input gain factor with a high speed input multiplexer, and 
therefore not shown here. Simulations of these structures were carried out 
with the simulation parameters as was used for the four-channel higher-order 
modulators reported in Section 6.7.3. As lower-order modulators are prone 
to generate idle tones, identical dither was added to each of the quantizer for 
a fair comparison, where the dithering signal was uniformly distributed over 
[-l/2,+l/2j. The output spectra for the standard, two-channel TI and two- 
channel ZI-TI second-order AL modulators are demonstrated in Figure 
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6.23(a) and (b). As expected and suggested by the theory both two-channel 
TI and ZI-TI modulators resulted in a net SNR improvement of 
approximately 15 dB in comparison to their standard counterparts. It should 
also be noted that the interleaved modulators derived out of the proposed 
method and block digital filtering are functionally equivalent when k = 1 and 
therefore they produce the same output. 




(a) 




(b) 

Figure 6-22. Second-order two-channel ZI-TI AZ modulators: (a) proposed, (b) block digital 

filtering 

This simulation experiment was repeated for different input amplitudes 
and each time the in-band SNR was measured, where the results are 
displayed in Figure 6.24(a). In order to prove the ZI-TI approach also works 
for low values of OSR, this simulation experiment was repeated for low 



212 



TIME-INTERLEAVED DELTA-SIGMA MODULATORS 



values of OSR (i.e., OSR=^ and 8, OSRe^% and 16), where this time the 
standard modulator was clocked at the effective sampling rate of the 
interleaved modulators. The results presented in Figure 6.24(b) confirm our 
assertion that the ZI-TI structure is also suitable for low OSR applications 
(i.e., wide bandwidth low-resolution). 
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Figure 6-23. Output spectra for second-order AT modulators: (a) standard vs. two-channel TI, 
(b) standard vs. two-channel ZI-TI 
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Input Amplitude (dB) 




Figure 6-24. Simulated in-band SNR performances of second-order modulators (standard, 2- 
channel TI and 2-channel Zl-TI): (a) OSR=32, 05/?e/=64. (b) 05/?=4, OSR,jr% and 

OSR,j=l6 
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Figure 6-25. Third-order four-channel proposed ZI-TI AX modulator 
120 




Input Amplitude (dB) 



Figure 6-26. Simulated in-band SNR performances of third-order modulators formed from the 
CI-DF topology (standard, and four-channel proposed ZI-TI) 
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In order to further confirm the proposed ZI-TI technique’s validity, we 
have also undertaken design and simulation studies on a third-order AS 
modulator with CI-DF topology as shown in Figure 6.6(a), where its four- 
channel ZI-TI counterpart is given in Figure 6.25. This third order modulator 
has a 70% depth of modulation (i.e., it is stable up to 0.7 full scale) [98]. 
Both the standard and four-channel ZI-TI modulators were simulated with 
different input amplitudes up to 0.7 (-3.01 dB full scale) and the in-band 
SNR plot against input amplitude is presented in Figure 6.26. On the 
average, the four-channel ZI-TI third-order modulator provides 
approximately 42 dB SNR improvement as expected, validating the 
functionality of the proposed ZI-TI concept for higher interleaving number 
M as well. 



6.9 Further Practical Issues 



6.9.1 Sampling Clock Jitter Effects 

A mathematical analysis to estimate the total error power due to the clock 
jitter on the input for the case of a single tone sinusoidal excitation was 
previously accomplished in [92] for the Nyquist rate TI converters. This 
analysis can be easily extended to TI AZ modulators to derive a formula for 
the total jitter error power, with the simple inclusion of the oversampling 
factor. In an M-channel TI modulator the total jitter error power at the output 
can be expressed as [92], [104]: 



E 



i.n ~ 



20 ^ 



(6.46) 



where A, and Oy, represent the magnitude and frequency of the input 
sine-wave, and standard deviation of the clock jitter error, respectively. 

To illustrate the SNR degradation due to sampling clock jitter, the 
second-order modulator shown in Figure 6.5(a), its two-channel TI versions 
in Figure 6.5(b) and (c), and its two-channel ZI-TI versions in Figure 6.22(a) 
and (b) were simulated under ideal and clock jitter conditions. Here, the 
sampling rate of the standard modulator (and hence the operation speed of 
the channels in the TI structure) was set to 2.048 MHz (i.e., A = 2.048 MHz) 
corresponding to OSR=6A and OSRejf=\2% for the standard and two-channel 
TI modulators, respectively. The input here was again an 8 KHz single tone 
sinusoidal signal. The second-order modulators were simulated for different 
input amplitudes and the result of this experiment is displayed in Figure 
6.27(a). 
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(a) 




(b) 



Figure 6-27. Simulated in-band SNR performances of second-order modulators under ideal 
and sampling clock jitter error conditions: (a) standard vs. two-channel Ti, (b) standard vs. 

two-channel ZI-TI 

We assumed that the sampling clockjitter is uniformly distributed over [- 
ATI2,+ATI2] (i.e., <7]~^UI\1) with NT~TjA where r, is the sampling 
period. The simulation results revealed that the standard modulator is almost 
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insensitive to this type of clock jitter, whereas its two-channel TI counterpart 
exhibits significant SNR degradation especially for large input amplitudes 
(e.g. A>0.1). The maximum in-band SNR is decreased by approximately 6.5 
dB when compared to the ideal case as seen from Figure 6.27(a). 

The sensitivity to sampling clock jitter error is surprisingly increased for 
the ZI-TI approach. Even though, there is only one component that 
contributes to jitter error for the ZI-TI modulator, the gain factor of M at the 
input stage increases the jitter power by a factor . Hence, the total jitter 
error power at the output of M-channel ZI-TI modulator is: 



p 






20SR. 



eff 



(6.47) 



The two-channel ZI-TI second-order modulators of Figure 6.22 were also 
simulated under the ideal and jitter error conditions in the same way as was 
done for the standard and two-channel TI modulators, where the SNR plot is 
given in Figure 6.27(b). These results indicate that the two-channel ZI-TI 
second-order modulator is marginally more sensitive than the TI modulator 
as expected and suggested by equations (6.46) and (6.47). The maximum in- 
band SNR for the jittered ZI-TI structure deteriorates by approximately 8 dB 
when compared to the ideal case. 

Table 6.3. summarizes the peak in-band SNR measurements for second- 
order TI and ZI-TI modulators under ideal and clock jitter conditions for 
different values of input signal bandwidth. Upon this analysis, we have 
established that the deterioration in peak in-band SNR is reduced with 
decreased OSR. This is because of the fact that the in-band quantization 
noise power is substantially increased in comparison to the increase of the 
jitter error power. 

Table 6-2. Measured peak in-band SNR performances of ideal and jittered second-order AS 
modulators for different values of input signal bandwidth ( /j„ = 8 kHz, = 2.048 MHz, 

=4.096 MHz and a] =r//l92 ) 





BW=32 kHz 
OSR=32 and OSR,.a-64 


nw=64 kHz 
OSR=16aiid OSR,,fl-32 


ideal 

jittered 


TI 

67.46 dB 
66.54 dB 


ZI-TI 
67.20 dB 
65.45 dB 


TI ZI-TI 

52.92 dB 52.82 dB 

52.63 dB 52.54 dB 



6.9.2 Branch Mismatch Effects 

A second major problem in the analog implementation of the interleaved 
modulators is the mismatches among various channels, which creates 



218 



TIME-INTERLEAVED DELTA-SIGMA MODULATORS 



aliasing distortions in the final output and significantly reduces the in-band 
SNR performance. In [95], branch mismatches were investigated thoroughly 
and the ^-factor method was introduced to alleviate the effects of this 
problem. The two-channel second-order block digital filter TI (Figure 
6.5(b)) and ZI-TI (Figure 22(b)) modulators were simulated with 0.1% 
mismatch ratio among the channels. Figure 6.28(a) depicts the output 
spectrum of the block digital filter TI modulator for OSR=6A and 
OSRej}=\2^> As expected, under the mismatch condition the noise floor is 
flattened and in-band SNR reduces by approximately 8 dB compared to the 
ideal case. The results for the block digital filter ZI-TI modulator are shown 
in Figure 6.28(b). In both of the regular TI and ZI-TI block digital filter 
structures, the A:-factor method is very effective in correcting the aliased 
noise floor. The in-band SNR performance of the two-channel proposed TI 
and ZI-TI (Figure 6.5(c)and Figure 6.22(a)) modulators with 0.1% 
mismatches is reduced to nearly the same value as that of the block digital 
filter modulators. Unfortunately, the A:-factor method is not applicable in our 
proposed modulators with their reduced complexity as discussed in Section 
6.3. Hence, we have concluded that, while the proposed modulators have 
simpler hardware, they are less robust against branch mismatches when 
compared to the cited block digital filtering idea. 

6.10 Conclusion 

In this chapter, a novel approach for the design of efficient architectures 
for TI AS modulators has been presented. Our proposed architectures result 
in modular structures and provide considerable reduction in the hardware 
requirements. The approach proposed here does not bring any restriction on 
the type of the modulator. However, it may not necessarily give simpler 
structures than the already existing approach based on the block digital 
filtering, for all cases (i.e., when the cross-coupling term is reduced from its 
nominal value of unity). Theoretical analysis and simulation results have 
shown that our proposed TI modulators produce the same transfer functions 
as those of the block digital filtering approach. 

We have also investigated some practical issues such as branch 
mismatches and clock jitter problems. Upon the investigation of these non- 
idealities, we concluded that our proposed TI modulators are marginally 
sensitive to the finite gain of the Op-amp as well as to the mismatches in 
scaling coefficients, therefore may require the design and use of better 
analog building blocks. In the digital implementations, however, there is no 
circuit imperfections that adversely affect the functionality of the modulator, 
thus in this case the proposed may be very advantageous because of their 
reduced hardware requirements. From an implementation perspective, both 
of the proposed and the block digital filtering approaches result in a critical 
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delay M times longer (where M is the interleaving number) than that of the 
non-interleaved counterpart. This problem can be overcome by using the 
pipelining technique in the digital implementations and the clock staggering 
strategy in the analog implementations [95]. 
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Figure 6-28. Output spectra for two-channel interleaved second-order AZ modulators with 
0. 1 % branch mismatch: (a) block digital filter TI, (b) block digital filter ZI-TI, (c) proposed 

TI and ZI-TI 

This chapter has also introduced a novel parallel conversion method 
called as the ZI-TI technique. This technique provides further simplifications 
in the implementation by eliminating the input multiplexer block involved in 
the regular TI conversion. We note that the ZI-TI technique is equally 
applicable to both our proposed method and the block digital filtering 
approach. It was shown that the ^-factor method could be successfully 
utilized in ZI-TI modulators derived out of the block digital filtering 
approach, in the same way as it is used in the regular TI approach. Hence, 
the ZI-TI technique preserves the robustness feature of the block digital 
approach in the analog implementations. 
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